REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704-0188 

The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of 
information,  including  suggestions  for  reducing  the  burden,  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188), 
1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any 
penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently  valid  OMB  control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 

1.  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE 

07-31-2008  Biannual  Technical  Report 

1  3.  DATES  COVERED  (From  -  To) 

1/1/2008  -  6/30/2008 

4.  TITLE  AND  SUBTITLE 

Bi-annual  (1/1/2008-6/30/2008)  Performance/Technical  Report 
for  ONR  Award  under  Grant  N000 14-07- 1-0395 

MARSnet:  Mission-aware  Autonomous  Radar  Sensor  Network  for 

Future  Combat  Systems 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

N00014  -  07  — 1  -0395 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

Liang,  Qilian 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

University  of  Texas  at  Arlington 

Office  of  Sponsored  Projects 

PO  Box  19145 

Arlington,  TX  76019 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Office  of  Naval  Research 

875  N.  Randolph  St. 

One  Liberty  Center 

Arlington,  VA  22203-1995 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

ONR 

11.  SPONSOR/MONITOR’S  REPORT 

NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  Public  Release;  Distribution  is  Unlimited. 

13.  SUPPLEMENTARY  NOTES 


14.  ABSTRACT 

During  the  period  of  1/1/2008  —  6/30/2008,  we  performed  the  following  studies  on  radar  sensor  network:  1 .  favor  weak  fuzzy  logic 
systems:  theory  and  applications  to  threat  assessment  in  disparate  sensor  networks;  2.  cognitive  radar  sensor  networks:  a 
human-inspired  information  integration  and  fuzzy  logic  system-based  approach;  3.  wireless  channel  modeling  in  foliage 
environment  using  UWB  radar  versus  narrowband  radar;  4.  tri-phase  coded  waveforms:  design  and  applications  to  radar  systems. 


15.  SUBJECT  TERMS 

Radar  Sensor  Network,  UWB  Radar,  fuzzy  logic  system,  threat  assessment,  tri-phase  coded  waveforms,  UWB  channel. 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

18.  NUMBER 

19a.  NAME  OF  RESPONSIBLE  PERSON 

a.  REPORT  1  b.  ABSTRACT 

c.  THIS  PAGE 

ABSTRACT 

OF 

PAGES 

Qilian  Liang 

U  U 

U 

UU 

134 

19b.  TELEPHONE  NUMBER  (Include  area  code) 

817-272-1339 

Standard  Form  298  (Rev.  8/98) 

Prescribed  by  ANSI  Std.  Z39.18 


20080807  1 79 


Bi-annual  (1/1/2008-6/30/2008)  Performance/Technical  Report  for 
ONR  Award  under  Grant  N000 14-07- 1-0395 
MARSnet:  Mission-aware  Autonomous  Radar  Sensor  Network  for 

Future  Combat  Systems 

Qilian  Liang 

Department  of  Electrical  Engineering 
University  of  Texas  at  Arlington 
Arlington,  TX  76019-0016  USA 
Phone:  817-272-1339,  Fax:  817-272-2253 
E-mail:  liang@uta.edu 


Abstract 

During  the  period  of  1/1/2008  -  6/30/2008,  we  performed  the  following  studies  on  radar 
sensor  network: 

1.  favor  weak  fuzzy  logic  systems:  theory  and  applications  to  threat  assessment  in  disparate 
sensor  networks; 

2.  cognitive  radar  sensor  networks:  a  human-inspired  information  integration  and  fuzzy  logic 
system- based  approach; 

3.  wireless  channel  modeling  in  foliage  environment  using  UWB  radar  versus  narrowband 
radar; 

4.  tri-phase  coded  waveforms:  design  and  applications  to  radar  systems. 

1  Favor  Weak  Fuzzy  Logic  Systems:  Theory  and  Applications  to 
Threat  Assessment 

In  current  and  future  DoD  operational  environments,  such  as  the  Global  War  on  Terrorism  (GWOT) 
and  the  Maritime  Domain  Awareness  (MDA),  there  is  a  need  to  exploit  and  integrate  heterogeneous 
sensor  networks  based  information  from  multiple  sensors  across  a  variety  of  modalities,  media  types 
including  radar,  video  imagery  (static  and  dynamic),  audio  (speech),  acoustic  stimuli,  IR,  text, 
as  well  as  signals  from  novel  chemical  and  biological  sensors.  Unfortunately,  current  DoD  tactical 
sensor  systems  are  not  networked  to  each  other  or  to  small  expeditionary  warfighters  units  operating 
on  the  battlefield.  A  divide  exists  between  receiving  timely  and  useful  tactical  sensor  data  and 
the  delivery  of  actionable  intelligence  to  the  Marine  Expeditionary  Force  (MEF)  and  subordinate 
expeditionary  force  units.  Moreover,  software  applications  do  not  exist  to  automatically  translate 
and  fuse  data  into  tactical  situational  understanding,  indications  and  warning  (I&W).  Algorithms 
do  not  exist  to  fuse  tactical  data  (level  2  or  3  fusion  for  a  GWOT  asymmetric  environment)  in  the 
expeditionary  forces  tactical  environment  using  tactically  available  computing  capabilities.  This 
need  is  also  motivated  by  the  fact  that  humans  display  a  remarkable  capability  to  perform  situation 
understanding  despite  noisy  sensory  signals  and  conflicting  inputs.  Humans  are  adept  at  network 
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visualization,  and  at  understanding  subtle  implications  among  the  network  connections.  To  date, 
however,  human’s  innate  ability  to  process  and  integrate  information  from  disparate,  network-based 
sources  for  situational  understanding  has  not  translated  well  to  automated  systems. 

Motivated  by  the  above  challenges,  in  [1],  we  applied  human  information  integration  mechanisms 
to  threat  assessment  in  heterogeneous  sensor  networks.  Humans  use  multiple  sources  of  sensory 
information  to  estimate  environmental  properties  and  has  innate  ability  to  integrate  information 
from  heterogeneous  data  sources.  How  the  multi-sensory  and  multimodal  information  are  integrated 
in  human  brain?  There  is  consensus  that  it  depends  on  the  prefrontal  cortex  (PFC).  The  PFC  has 
top-down  control  (favor  weak)  and  rule-based  mechanisms,  and  we  propose  to  incorporate  the  favor 
weak  mechanism  into  rule-based  fuzzy  logic  systems  (FLS)  via  using  upper  and  lower  membership 
functions.  The  inference  engine  of  favor  weak  fuzzzy  logic  system  is  proposed  under  three  different 
categories  based  on  fuzzifiers.  We  observe  that  the  favor  weak  FLS  is  a  special  type-1  FLS  which 
is  embeded  in  an  interval  type-2  FLS,  so  its  much  simpler  in  computing  than  an  interval  type-2 
FLS.  We  applied  the  favor  weak  FLS  to  situation  understanding  based  on  heterogeneous  sensor 
network,  and  it  shows  that  our  favor  weak  fuzzy  logic  system  has  clear  advantage  comparing  to 
the  type-1  FLS.  The  favor  weak  FLS  can  increase  the  probability  of  threat  detection,  and  provides 
timely  I&W. 

2  Cognitive  Radar  Sensor  Networks:  A  Human-Inspired  Informa¬ 
tion  Integration  and  Fuzzy  Logic  System-Based  Approach 

Inspired  by  human’s  innate  ability  to  process  and  integrate  information  from  network-based  sources, 
we  applied  [2]  [4]  human-inspired  information  integration  mechanisms  to  target  detection  in  cognitive 
radar  sensor  network.  Humans’  information  integration  mechanisms  have  been  modelled  using 
maximum-likelihood  estimation  (MLE)  or  soft-max  approaches.  In  [2]  [4],  we  applied  these  two 
algorithms  to  cognitive  radar  sensor  networks  target  detection.  Discrete-cosine-transform  (DCT) 
was  used  to  process  the  integrated  data  from  MLE  or  soft-max.  We  applied  fuzzy  logic  system 
(FLS)  to  automatic  target  detection  based  on  the  AC  power  values  from  DCT.  Simulation  results 
showed  that  our  MLE-DCT-FLS  and  soft-max-DCT-FLS  approaches  performed  very  well  in  the 
radar  sensor  network  target  detection,  whereas  the  existing  2-D  construction  algorithm  didn’t  work 
in  this  study. 

3  Wireless  Channel  Modeling  in  Foliage  Environment:  UWB  ver¬ 
sus  Narrowband 

In  [3],  we  studied  the  wireless  channel  modeling  in  foliage  environment,  a  rich  scattering  and  time- 
varying  environment,  based  on  extensive  data  collected  using  UWB  and  narrowband  (200MHz 
and  400MHz)  radars.  We  applied  two  approaches  to  the  wireless  channel  modeling:  Saleh  and 
Valenzuela  (S-V)  method  for  UWB  channel  modeling  and  CLEAN  method  for  narrowband  and 
UWB  channel  modeling.  We  validated  that  UWB  echo  signals  (within  a  burst)  don’t  hold  self¬ 
similarity,  which  means  the  future  signals  can’t  be  forecasted  based  on  the  received  signals  and 
channel  modeling  is  necessary  from  statistical  point  of  view.  Based  on  the  S-V  method  for  UWB 
channel  modeling,  in  foliage  UWB  channel,  the  multi-path  contributions  arrive  at  the  receiver  are 
grouped  into  clusters.  The  time  of  arrival  of  clusters  can  be  modeled  as  a  Poisson  arrival  process, 
while  within  each  cluster,  subsequent  multipath  contributions  or  rays  also  arrive  according  to  a 
Poisson  process.  At  different  field  (near  field,  medium  field,  and  far  field),  we  observe  that  the 
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Poisson  process  parameters  are  quite  different.  We  also  observe  that  the  amplitude  of  channel 
coefficient  at  each  path  follows  Rician  distribution  for  medium  and  far  field,  and  it’s  non-st  at  ionary 
for  paths  from  near  field  (one  of  two  Rician  distributions),  and  these  observations  are  quite  different 
with  the  IEEE  indoor  UWB  channel  model  and  S-V  indoor  channel  model.  Based  on  the  CLEAN 
method,  the  narrowband  (200MHz  and  400MHz)  and  UWB  channel  impulse  responses  have  many 
similarities:  both  can  be  modeled  as  linear  time-variant  filter  channel. 

In  [5]  [6],  we  studied  the  statistical  modeling  for  outdoor  non  line-of-sight  (NLOS)  channel 
in  rich  scattering  and  time-varying  foliage  environment  based  on  extensive  data  collected  using 
both  narrowband  and  ultra-wideband  (UWB)  radars.  The  multi-path  contributions  arrive  at  the 
receiver  are  grouped  into  clusters.  The  time  of  arrival  of  clusters  can  be  modeled  as  a  Poisson  arrival 
process,  while  within  each  cluster,  subsequent  multipath  contributions  or  rays  also  arrive  according 
to  a  Poisson  process.  However,  the  parameters  are  quite  different  along  with  the  frequency.  We  also 
observe  that  the  amplitude  of  channel  coecient  at  each  path  can  be  more  accurately  characterized 
by  log-logistic  distribution  (LLD)  other  than  log-normal,  Weibull  or  Rayleigh  due  to  the  best 
goodness-of-fit  and  smallest  root-mean-square  (RMS). 

4  Tri-Phase  Coded  Waveforms:  Design  and  Applications  to  Radar 
Systems 

A  phase  coded  waveform  has  a  constant  RF  frequency,  but  an  absolute  phase  that  is  switched 
between  one  of  N  fixed  values  at  regular  intervals  within  the  pulse  length.  In  this  project,  we 
proposed  two  tri-phase  coded  waveform  ( N  =  3)  design  methods. 

In  [7],  we  presented  [7]  the  definition  and  properties  of  ZCZ  sequence-pair  set  based  on  zero  cor¬ 
relation  zone  (ZCZ)  concept,  and  proposed  to  use  the  perfect  punctured  sequence-pair  together  with 
Hadamard  matrix  in  the  zero  correlation  zone  to  construct  the  perfect  punctured  ZCZ  sequence- 
pair  set  for  tri-phase  coded  waveform  design.  According  to  performance  analysis,  perfect  punctured 
ZCZ  sequence-pair  set  has  good  autocorrelation  and  cross  correlation  properties  when  doppler  shift 
is  not  large.  The  radar  target  detection  system  simulation  results  also  show  that  perfect  punc¬ 
tured  ZCZ  sequence-pair  outperforms  other  conventional  pulse  compression  codes,  such  as  the  well 
known  polyphase  codeP4  code.  As  a  result,  perfect  punctured  ZCZ  sequence-pair  set  can  be  good 
candidate  for  pulse  compression  code. 

In  [8],  we  proposed  new  code  punctured  binary  sequence-pair  and  applied  it  to  tri-phase  coded 
waveforms.  The  definitions  and  the  autocorrelation  properties  of  the  proposed  code  are  given. 
Doppler  shift  performance  is  also  investigated.  The  significant  advantages  of  punctured  binary 
sequence-pair  over  conventional  pulse  compression  codes,  such  as  the  widely  used  Barker  codes,  are 
zero  autocorrelation  sidelobes  and  the  longer  length  of  the  code  which  can  be  as  long  as  31  so  fax. 
In  the  radar  target  detection  system  simulation,  punctured  binary  sequencepair  also  outperforms 
other  conventional  pulse  compression  codes.  Therefore,  our  proposed  code  can  be  used  as  one  of 
the  candidates  for  pulse  compression  code. 
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Abstract 

Humans  use  multiple  sources  of  sensory  information  to  estimate  environmental  proper¬ 
ties  and  has  innate  ability  to  integrate  information  from  heterogeneous  data  sources.  How 
the  multi-sensory  and  multimodal  information  are  integrated  in  human  brain?  There  is 
consensus  that  it  depends  on  the  prefrontal  cortex  (PFC).  The  PFC  has  top-down  control 
(favor  weak)  and  rule-based  mechanisms,  and  we  propose  to  incorporate  the  favor  weak 
mechanism  into  rule-based  fuzzy  logic  systems  (FLS)  via  using  upper  and  lower  membership 
functions.  The  inference  engine  of  favor  weak  fuzzzy  logic  system  is  proposed  under  three 
different  categories  based  on  fuzzifiers.  We  observe  that  the  favor  weak  FLS  is  a  special 
type-1  FLS  which  is  embeded  in  an  interval  type-2  FLS,  so  it’s  much  simpler  in  computing 
than  an  interval  type-2  FLS.  We  apply  the  favor  weak  FLS  to  situation  understanding 
based  on  heterogeneous  sensor  network,  and  it  shows  that  our  favor  weak  fuzzy  logic  sys¬ 
tem  has  clear  advantage  comparing  to  the  type-1  FLS.  The  favor  weak  FLS  can  increase 
the  probability  of  threat  detection,  and  provides  timely  indication  &;  warning  (I&W). 
Index  Terms  :  Fuzzy  logic  systems,  prefrontal  cortex  (PFC),  favor  weak,  upper  and  lower 
membership  functions,  situation  understanding,  heterogeneous  sensor  networks. 
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1  Introduction  and  Motivation 


Humans  use  multiple  sources  of  sensory  information  to  estimate  environmental  properties.  For 
example,  the  eyes  and  hands  both  provide  relevant  information  about  an  objects  shape.  The 
eyes  estimate  shape  using  binocular  disparity,  perspective  projection,  etc.  The  hands  supply 
haptic  shape  information  by  means  of  tactile  and  proprioceptive  cues.  Combining  informa¬ 
tion  across  cues  can  improve  estimation  of  object  properties  but  may  come  at  a  cost: 
loss  of  single-cue  information.  Recent  studies  [12]  showed  that  single-cue  information  is  indeed 
lost  when  cues  from  within  the  same  sensory  modality  (e.g.,  disparity  and  texture  gradients 
in  vision)  are  combined,  but  not  when  different  modalities  (vision  and  haptics)  are 
combined.  In  another  study  on  human  [7],  gaze  shifts  are  coordinated  movements  of  the  eyes 
(eyes-r e-head)  and  head  (head-re-space)  that  rapidly  reorient  the  visual  axis  (eyes-re-space)  to 
a  target  of  interest.  Reaction  latencies  for  gaze  shifts  to  combined  auditory  and  visual  stimuli 
presented  in  close  spatial  and  temporal  register  are  less  than  those  to  either  stimulus 
presented  alone,  suggesting  that  the  integration  of  multisensory  information  may  play  an 
important  role  in  forming  appropriate  motor  behaviors.  These  studies  demonstrate  that  human 
has  innate  ability  to  integrate  information  from  heterogeneous  data  sources  and  multi-sensory 
and  multimodal  information  integration  has  clear  advantage. 

In  this  paper,  as  a  product  of  multidisciplinary  collaborative  research,  we  incorporate  hu¬ 
man  brain  mechanisms  to  a  new  fuzzy  logic  system  design  and  apply  it  to  situation  under¬ 
standing  based  on  heterogeneous  sensor  network.  A  heterogeneous  sensor  network  consists  of 
multiple  networked  sensors  with  different  modality  (video,  audio,  acoustic,  radar,  etc),  and 
such  networks  are  necessary  in  different  applications.  For  example,  in  an  emergency  natural 
disaster  scenario,  information  integration  for  first  responders  is  critical  for  search  and  rescue. 
Besides,  the  first  responders  need  to  be  situation-aware.  Danger  may  appear  anywhere  at 
any  time,  therefore,  first  responders  must  monitor  a  large  area  continuously  in  order  to  iden- 
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tify  potential  danger  and  take  actions.  Due  to  the  dynamic  and  complex  nature  of  natural 
disaster,  some  victims  may  not  be  found  with  a  single  type  of  sensor  modality,  for  example, 
image/video  sensors  can’t  be  used  to  find  a  buried/foleage  victim,  UWB  radar  sensors  need  to 
be  used  for  penetrating  the  ground  or  sense-through-wall,  and  acoustic  sensors  are  needed  to 
collect  the  voice  from  victims.  Similarly,  some  potential  dangers  may  not  be  identified  using 
a  single  modality  sensor.  More  modalities  are  required  to  search  victims  and  identify  poten¬ 
tial  dangers  and  that  means  large-scale  Heterogeneous  Sensor  Networks  (HSN)  are  needed  for 
search,  rescue,  and  situation  awareness.  However,  information  integration  algorithms  (espe¬ 
cially  for  situation  awareness)  for  heterogeneous  sensor  networks  don’t  exist.  Motivated  by  the 
above  challenges,  we  study  human-inspired  information  integration  for  heterogeneous  sensor 
networks. 

The  remaining  of  the  paper  is  organized  as  follows.  In  Section  2,  we  present  the  heteroge¬ 
neous  information  integration  in  Human  brain  and  challenges  to  fuzzy  logic  system  design.  In 
Section  3,  we  give  an  overview  on  upper  and  lower  membership  design  and  describe  how  it  can 
be  used  to  the  new  fuzzy  logic  system  design.  In  Section  4,  we  present  the  inference  engine 
for  favor  weak  fuzzy  logic  system,  and  its  relations  with  the  interval  type-2  fuzzy  logic  systems 
are  described  in  Section  5.  The  application  to  situation  understanding  based  on  heterogeneous 
sensor  network  is  presented  in  Section  6.  Section  7  concludes  this  paper. 

2  Heterogeneous  Information  Integration  in  Human  Brain  and 
Challenges  to  Fuzzy  Logic  System  Design 

One  of  the  great  mysteries  of  the  brain  is  cognitive  control.  How  can  the  interactions  between 
millions  of  neurons  result  in  behavior  that  is  coordinated  and  appears  willful  and  voluntary? 
There  is  consensus  that  it  depends  on  the  prefrontal  cortex  (PFC)  [31] [35].  A  schematic  dia¬ 
gram  of  some  of  the  extrinsic  and  intrinsic  connections  of  the  PFC  is  depicted  in  Fig.  1  [31]. 
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Many  PFC  areas  receive  converging  inputs  from  at  least  two  sensory  modalities  [4]  [15].  For 
example,  the  dorsolateral  (DL)  (areas  8,  9,  and  46)  and  ventrolateral  (12  and  45)  PFC  both 
receive  projections  from  visual,  auditory,  and  somatosensory  cortex.  Furthermore,  the  PFC 
is  connected  with  other  cortical  regions  that  are  themselves  sites  of  multimodal  convergence. 
Many  PFC  areas  (9,  12,  46,  and  45)  receive  inputs  from  the  rostral  superior  temporal  sul¬ 
cus,  which  has  neurons  with  bimodal  or  trimodal  (visual,  auditory,  and  somatosensory)  re¬ 
sponses  [1]  [36] .  The  arcuate  sulcus  region  (areas  8  and  45)  and  area  12  seem  to  be  particularly 
multimodal.  They  contain  zones  that  receive  overlapping  inputs  from  three  sensory  modalities 
[36].  Observe,  for  example,  that  mid-dorsal  area  9  directly  processes  and  integrates  visual, 
auditory,  and  multimodal  information. 

2.1  PFC  Top-Down  Control  (Favor  Weak)  and  Rule-Based  Mechanisms 

According  to  [31]  [35],  the  PFC  is  modulatory  rather  than  transmissive.  That  is,  the  pathway 
from  input  to  output  does  not  “run  through”  the  PFC.  Instead,  the  PFC  guides  activity 
flow  along  task-relevant  pathways  in  more  posterior  and/or  subcortical  areas.  The  PFC  is 
important  when  “top-down”  processing  is  needed;  that  is,  when  behavior  must  be  guided  by 
internal  states  or  intentions.  The  PFC  is  critical  in  situations  when  the  mappings  between 
sensory  inputs,  thoughts,  and  actions  either  are  weakly  established  relative  to  other  existing 
ones  or  are  rapidly  changing.  This  is  when  we  need  to  use  the  “ rules  of  the  game,”  internal 
representations  of  goals  and  the  means  to  achieve  them  [31].  Several  investigators  have  argued 
that  this  is  a  cardinal  function  of  the  PFC  [6]  [37]  [11]  [41]  [30].  The  top-down  control  and  favor 
weak  mechanism  can  be  illustrated  using  the  Stroop  task  Wisconsin  card  sort  task  (WCST). 

In  the  Stroop  task  [40]  [27],  subjects  either  read  words  or  name  the  color  in  which  they 
are  written.  To  perform  this  task,  subjects  must  selectively  attend  to  one  attribute.  This  is 
especially  so  when  naming  the  color  of  a  conflict  stimulus  (e.g.  the  word  GREEN  displayed  in 
red),  because  there  is  a  strong  prepotent  tendency  to  read  the  word  (“green”),  which  competes 

4 


Page  8  of  134 


with  the  response  to  the  color  (“red”).  This  illustrates  one  of  the  most  fundamental  aspects 
of  cognitive  control  and  goal-directed  behavior:  the  ability  to  select  a  weaker,  task-relevant 
response  (or  source  of  information)  in  the  face  of  competition  from  an  otherwise  stronger,  but 
task-irrelevant  one.  Patients  with  frontal  impairment  have  difficulty  with  this  task,  especially 
when  the  instructions  vary  frequently,  which  suggests  that  they  have  difficulty  adhering  to  the 
goal  of  the  task  or  its  rules  in  the  face  of  a  competing  stronger  (i.e.  more  salient  or  habitual) 
response  [31].  Similar  findings  are  evident  in  the  WCST  [32].  Subjects  are  instructed  to  sort 
cards  according  to  the  shape,  color,  or  number  of  symbols  appearing  on  them  and  the  sorting 
rule  varies  periodically.  Thus,  any  given  card  can  be  associated  with  several  possible  actions, 
no  single  stimulus-response  mapping  will  work,  and  the  correct  one  changes  and  is  dictated  by 
whichever  rule  is  currently  in  effect.  Humans  with  PFC  damage  show  stereotyped  deficits  in 
the  WCST.  They  are  able  to  acquire  the  initial  mapping  without  much  difficulty  but  are  unable 
to  adapt  their  behavior  when  the  rule  varies  [32].  Monkeys  with  PFC  lesions  are  impaired  in 
an  analog  of  this  task  and  in  others  when  they  must  switch  between  different  rules  [31]. 

The  Stroop  task,  naming  the  color  of  a  conflict  stimulus,  and  WCST  [32]  are  variously 
described  as  tapping  the  cognitive  functions  of  either  selective  attention,  behavioral  inhi¬ 
bition,  top-down  control,  working  memory,  or  rule-based  or  goal-directed  behavior 
[31].  As  suggested  by  Desimone  and  Duncan  [9],  selective  attention  and  behavioral  inhibition 
are  two  sides  of  the  same  coin:  attention  is  the  effect  of  biasing  competition  in  favor  of  task¬ 
relevant  information,  and  inhibition  is  the  consequence  that  this  has  for  the  irrelevant  informa¬ 
tion.  In  this  project,  we  will  study  human  brain  top-down  control  and  rule-based  mechanisms 
inspired  information  integration.  In  current  HSN  design,  the  clusterhead  only  serves  a  “trans¬ 
missive”  (data  collection  and  relay  to  gateway)  function.  In  this  project,  the  HSN  clusterhead 
will  also  provide  a  “modulatory”  function,  i.e.,  multimodal  information  integration.  In  natural 
disaster  or  terrorist  attack  recovery,  the  most  dangerous  potential  threat  factors  (stimuli)  are 
weak  or  hidden  but  are  highly  correlated  with  the  situation  understanding  task,  and  of  course, 
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they  are  time  sensitive  (rapidly  changing),  which  indicates  that  the  PFC  “top-down”  control 
mechanism  can  be  applied  to  HSN-based  situation  understanding. 

The  PFC  top-down  control  signals  favor  weak  ( but  task-relevant)  stimulus-response  map¬ 
pings  when  they  are  in  competition  with  more  habitual,  stronger  ones  (as  in  the  Stroop  task, 
where  the  word  GREEN  is  stronger  and  the  color  red  is  weak),  especially  when  flexibility  is 
needed  (such  as  in  the  WCST)  [31].  Moreover,  all  of  the  PFC  neural  mechanisms  depend  on  the 
representation  of  goals  and  rules  in  the  form  of  patterns  of  activity  in  the  PFC,  which  configure 
processing  in  other  parts  of  the  brain  in  accordance  with  current  task  demands  [31]  [35].  Such 
mechanisms  motivate  us  to  heavily  revisit  a  rule-based  approach:  fuzzy  logic  systems  (FLS), 
mimicking  the  rule-based  PFC  neural  mechanism,  and  subsequently  applying  it  to  HSN-based 
situation  awareness. 

2.2  Overview  of  Fuzzy  Logic  Systems  and  Its  Shortfall 

The  current  type-1  FLS  designs  doesn’t  have  “favor  strong  or  favor  weak  control”.  In  a  type-1 
FLS  with  a  rule  base  of  M  rules,  in  which  each  rule  has  p  antecedents,  let  the  Zth  rule  be 
denoted  by  Rl,  where  Rl:  IF  xi  is  F,,  and  x2  is  Fj,  and,  and  xp  is  Flp  THEN  y  is  GJ. 
The  membership  function,  /iB/(y),  of  a  fired  rule  can  be  expressed  by  the  following  sup-star 
composition  [29]: 

HB,(y)  =  supx€A<  [/iA>  (x)  *  /rAi_Bi  (x,  j/)]  (1) 

where  A*  is  a  p-dimensional  Cartesian  product  space,  A*  =  AJ  x  •  •  •  x  A* ,  A*k  is  the  measurement 
domain  of  input  xk,  (k  =  1,  •  •  •  ,p);  and,  A*  is  given  by 

Ma*(x)  =  Ma;x-xa;(x)  =  HXiixi)  *  ■  ■  •  *  nxp(xp)  (2) 

In  the  current  type-1  FLS  design, 

^Al—>Bl  ^  (®2)  *  ’  ‘ *  A*f|,  (®p)  *  Mg*  (j/)  (3) 
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For  the  Stroop  task,  to  name  the  color  of  a  conflict  stimulus  (e.g.  the  word  GREEN 
displayed  in  red),  needs  the  favor  weak  mechanism  because  there  is  a  strong  prepotent  tendency 
to  read  the  word  (“green”),  which  competes  with  the  response  to  the  color  (“red”).  In  Fig.  2, 
we  illustrate  the  schematic  of  the  Stroop  model  using  the  example  to  name  the  color  when  the 
word  GREEN  displayed  in  red  [31]. 

So  in  the  FLS  design,  the  firing  degree  for  “red”  should  be  boosted  to  reflect  the  “favor 
weak”  mechanism  if  a  conflict  stimulus  presents,  but  no  control  should  be  taken  if  no  con¬ 
flict  stimulus  presents  (e.g.,  name  the  word  when  the  word  GREEN  displayed  in  red).  This 
motivates  us  to  use  different  membership  degrees  under  different  scenarios.  We  propose  to 
use  interval  type-2  fuzzy  membership  function  for  this  favor  weak  (using  upper  membership 
function)  or  no  control  (using  lower  membership  function)  mechanism. 


3  Upper  and  Lower  Membership  Functions 


An  upper  MF  and  a  lower  MF  are  two  type-1  MFs  which  are  bounds  for  the  footprint  of 
uncertainty  of  an  interval  type-2  MF.  The  upper  MF  is  a  subset  which  has  the  maximum 
membership  grade  of  the  footprint  of  uncertainty;  and,  the  lower  MF  is  a  subset  which  has  the 
minimum  membership  grade  of  the  footprint  of  uncertainty  [17]. 

Same  as  that  in  [17],  we  use  an  overbar  (underbar)  to  denote  the  upper  (lower)  MF.  For 
example,  the  upper  and  lower  MFs  of  [i-i  ( X*. )  are  Jlxi  ( x* )  and  /ixi  (x^),  respectively,  so  that 

Ait  Afc  — Ait 


=  f 


(xk),n-i  (xfc)] 
Afc  Afc 


1/9' 


(4) 


Example  1:  Gaussian  Primary  MF  with  Uncertain  Standard  Deviation 

Consider  the  case  of  a  Gaussian  primary  MF  having  a  fixed  mean,  mj^,  and  an  uncertain 
standard  deviation  that  takes  on  values  in  [<7^, c^],  i.e., 


Mfc(zfc)  =  exp  --( 


1  ,xk-mlk2 


)2  . 


uk 
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where:  k  =  1, . . .  ,p;  p  is  the  number  of  antecedents;  l  =  1, . . . ,  M;  and,  M  is  the  number  of 
rules.  The  upper  MF,  JLlk(xk),  is  (see  Fig.  3b) 

A(xk)  =  M(mlk,  <t£2;  xk),  (6) 

and  the  lower  MF,  ^(x*),  is  (see  Fig.  3) 

/ffc(zfc)  =  <rlkl;xk)  (7) 

□ 

This  example  illustrates  how  to  define  Jl  and  /i,  so  that  it  is  clear  how  to  define  these 
membership  functions  for  other  situations  (e.g.,  triangular,  trapezoidal,  bell  MFs). 

4  The  Inference  Engine  for  Favor  Weak  FLSs 

In  a  favor  weak  FLS  with  a  rule  base  of  M  rules,  in  which  each  rule  has  p  antecedents,  let 
the  Zth  rule  be  denoted  by  Rl,  where  Rl :  IF  x\  is  Fj,  and  X2  is  F2,  and,  and  xv  is  Fp 
THEN  y  is  Although  the  rule  may  look  the  same  as  that  of  the  interval  type-2  FLS,  but 
the  inference  engine  is  different. 

In  a  favor  weak  FLS  with  p  antecedents,  without  loss  of  generality,  assume  the  first  w  (w  < 
p)  antecedents  are  weak  and  should  be  favored,  and  all  other  antecedents  (w  +  1,  iw+ 2,  •  •  •  ,  ip) 
should  be  in  no  control.  If  the  antecedents  are  not  in  this  order,  they  can  be  re-orderred.  Then 
we  can  obtain  the  following  Theorem.  Our  major  result  for  favor  weak  FLSs  is  given  in: 

Theorem  1  In  a  favor  weak  nonsingleton  FLS  (the  first  w  antecedents  are  weak  and  should  be 
favored)  with  type-2  fuzzification  and  meet  under  minimum  or  product  t-norm:  (a)  the  firing 
degree  for  rule  l,  i.e .,  the  result  of  the  input  and  antecedent  operations,  is 

f  —  sup  /  /  [tt'Xi (^i)  *  0^1 )]  *  *  *  *  *  [flx *  1^]^  (®ti/)]  (8) 

J X 1  JXP  1 

*^xu,+i(x"’+1)  *^+1(x«'+1)]  *  •  •  •  *  \Hxp(.xp)  (9) 
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the  supremum  is  attained  when  each  term  in  brackets  attains  its  supremum; 

(b)  the  rule  Rl  fired  output  consequent  set ,  is 

i*Bf(v)  =  (10) 

where  [£^{y)  and  are  the  lower  and  upper  membership  grades  of  p^(y);  and, 

(c)  the  output  fuzzy  set,  PB{y),  is 

M 

vijiy)  =  {Jv&iy)  (ii) 

/=1 

When  the  input  is  fuzzified  to  a  type-1  fuzzy  set,  so  that  — ►  pxk  (A:  =  1, .. .  ,p),  the 

upper  and  lower  MFs  of  merge  into  one  MF,  pxk(xk),  in  which  case  Theorem  1  simplifies 
to: 

Corollary  1  In  a  favor  weak  FLS  (the  first  w  antecedents  are  weak  and  should  be  favored) 
with  nonsingleton  type-1  fuzzification  and  meet  under  minimum  or  product  t-norm,  fl  in  (9) 
simplifies  to:  the  firing  degree  for  rule  l  is 

fl  =  sup  /  •••  /  (x„,)]  (12) 

x€XJXi  Jxp 

*[vxw+1{xw+i)  *  Hpi  (Xw+i)]*---*[nxp(xp)*nu(xp)]/x-,  (13) 

1  i/j-f- 1  — r  p 

the  supremum  is  attained  when  each  term  in  brackets  attains  its  supremum; 

When  a  singleton  fuzzifier  is  used,  the  upper  and  lower  MFs  of  A*xfc(xfc)  merge  into  one 
crisp  value,  namely  1,  in  which  case  Theorem  1  simplifies  further  to: 

Corollary  2  In  a  favor  weak  FLS  (the  first  w  antecedents  are  weak  and  should  be  favored) 
with  singleton  fuzzification  and  meet  under  minimum  or  product  t-norm,  fl  in  (9)  simplifies 
to:  the  firing  degree  for  rule  l,  i.e the  result  of  the  input  and  antecedent  operations, 

f  A^j-A  (*^i)  *  ’  *  *  *  A^  (xw)  *  A^  (*Etu-fi)  *  *  *  *  *  Hu  (*^p)  (14) 

1  rw  "u>+l  — rp 

where  Xi  (i  =  1, . . .  ,p)  denotes  the  location  of  the  singleton. 
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5  Relations  between  Favor  Weak  FLS  and  Interval  Type-2  FLS 


In  decision  theory,  ambiguity  about  probabilities  should  not  affect  choices.  However,  recent 
experiments  [13]  showed  that  many  people  are  more  willing  to  bet  on  risky  outcomes  (e.g., 
gambling  on  a  roulette  wheel)  than  on  an  ambiguous  one  (e.g.,  chance  of  a  terrorist  attack 
based  on  meager  or  conflicting  evidence),  holding  the  judged  probability  of  outcomes  constant. 
So,  the  confidence  in  judged  probability  can  vary  widely  for  “risky”  and  “ambiguous”.  Using 
functional  brain  imaging,  Hsu  et  al  [13]  showed  that  the  level  of  ambiguity  in  choices  correlates 
positively  with  activation  in  the  amygdala  and  orbitofrontal  cortex,  and  negatively  with  a  stri¬ 
atal  system.  This  suggests  that  degree  of  uncertainty  should  be  considered  in  decision  making, 
contrary  to  traditional  decision  theory.  Type-2  fuzzy  sets  and  FLSs  are  successful  in  handling 
the  uncertainties  [17].  Type-2  fuzzy  sets  have  grades  of  membership  that  are  themselves  fuzzy. 
A  type-2  membership  grade  can  be  any  subset  in  [0, 1]  -  the  primary  membership ;  and,  corre¬ 
sponding  to  each  primary  membership,  there  is  a  secondary  membership  (which  can  also  be  in 
[0, 1])  that  defines  the  possibilities  for  the  primary  membership.  Figure  4  shows  an  example  of 
a  type-2  set.  The  domain  of  the  membership  grade  corresponding  to  x  =  4  is  also  shown. 

In  [17]  [20]  [22]  [23]  [24],  Liang  and  Mendel  proposed  the  theory  and  design  of  interval  type-2 
FLS.  With  their  pioneering  works,  people  are  able  to  efficiently  handle  uncertainties.  They 
applied  interval  type-2  FLS  to  a  number  of  very  important  applications  where  uncertainties 
abound,  such  as  fading  channel  equalization  [18]  and  co-channel  interference  elimination  [19], 
network  video  traffic  modeling  and  classification  [21],  connection  admission  control  for  ATM 
network  [16].  Recently,  Liang  and  his  students  applied  interval  type-2  FLSs  to  sensor  network 
lifetime  estimation[39][38],  event  forecasting  in  wireless  sensor  networks [25]  [26],  and  cross-layer 
optimization  in  ad  hoc  networks  [43]  [42]. 

In  this  paper,  PFC’s  “top-down  control  (favor  weak)”,  “rule-based”,  and  “brain  handling 
ambiguity”  mechanisms  are  incorporated  into  the  favor  weak  FLS  design  to  mimic  human 
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brain  learning  and  decision  making  for  associating,  integrating  and  understanding/inferencing 
discovered  knowledge  from  disparate  sources.  The  favor  weak  FLS  we  proposed  in  this  paper 
is  a  type-1  FLS  which  is  embeded  in  the  interval  type-2  FLS. 

6  Application  to  Situation  Understanding  Based  on  Heteroge¬ 
neous  Sensor  Network 

Some  valuable  work  has  been  reported  on  situation  understanding  (situation  awareness  and 
threat  assessment).  In  [10],  an  intelligent  threat  assessment  processor  using  genetic  algorithms 
and  fuzzy  logic  was  proposed.  In  [33],  threat  assessment  was  studied  in  tactical  airborne 
environments.  In  [14],  a  neural  network  was  applied  to  threat  assessment  for  automated  visual 
surveillance.  In  [8],  an  intelligent  assistant  was  proposed  to  provide  automatic  situation  and 
threat  advice  in  the  Air  Defence  Ground  Environment.  In  [2],  a  situation  and  threat  assessment 
model  based  on  group  analysis  was  proposed.  A  situation/ threat  assessment  fusion  system 
was  proposed  in  [5].  Other  approaches  that  have  appeared  include  multiple  attribute  decision 
making  [3],  bayesian  networks  [34],  etc.  Unfortunately,  none  of  these  approaches  is  appropriate 
for  the  multimodal  sensor  network  scenario  and  none  of  them  have  used  models  derived  from 
human  or  biological  system  mechanisms,  although  some  computational  intelligence  models 
(such  as  neural  networks  and  fuzzy  logic)  were  used. 

In  situation  awareness,  the  “weak”  can  be  interprated  as,  for  example,  1)  a  target  (or 
entity)  rarely  appears  in  the  sensor  field;  or  2)  the  behavior  pattern  of  this  target  (or  entity) 
has  low  match  with  the  existing  ones  in  database;  or  3)  space/time  correlation  of  data/entities 
to  events  is  low,  etc.  All  the  above  information  can  be  obtained  based  on  the  assumption  that 
high  quality  information  about  objects  and  events  is  available  as  a  contributor  to  situation 
awareness.  In  general,  however,  such  information  is  insufficient  to  provide  adequate  situation 
awareness.  Actually,  some  lessons  gained  from  terrorist  attacks  have  already  demonstrated 
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that  the  above  “weaks”  should  be  favored  in  situation  awareness.  Situation  awareness  needs 


the  “favor  weak”  and  rule-based  mechanisms  from  PFC. 

The  new  favor  weak  FLS  can  be  used  for  level  2/3  fusion  (situation  awareness  and  threat 
assessment)  for  inferring  activities,  relationships,  and  intentions  of  objects  and  people  in  the 
battlespace  based  on  retrieved  knowledge  consisting  of  behavioral  patterns,  new  activities, 
and  anticipated  behavior,  and  also  taking  into  account  contextual  information  (terrain,  roads, 
weather,  etc).  The  level  1  data  fusion  results  (traditional  automatic  target  recognition  and 
pattern  recognition)  from  multimodal  sensors  will  be  used  as  antecedents.  For  example,  con¬ 
sidering  a  heterogeneous  sensor  network  with  radar,  image/video  sensors,  and  GPS  sensor,  we 
can  choose  the  following  three  antecedents: 

1.  The  number  of  switches  from  non-maneuvering  set  (constant  behavior  in  speed,  accelera¬ 
tion,  and  direction,  etc)  to  the  maneuvering  set  (varying  behavior  in  speed,  acceleration, 
and  direction,  etc).  When  a  target  is  beginning  a  maneuver  from  a  non-maneuvering 
class,  the  tracking  system  can  switch  the  algorithms  applied  to  the  problem  from  a  non¬ 
maneuvering  set  to  the  maneuvering  set.  The  errors  in  distance  from  where  the  tracker 
estimates  the  position  of  a  target  between  the  actual  position  can  be  very  large  when 
the  incorrect  motion  models  are  applied  to  the  problem.  Additionally,  when  the  tracker 
does  finally  catch  up  to  the  target  after  the  maneuver,  the  track  will  “jump”  across  the 
operator’s  scope  giving  a  very  unrealistic  and  unreliable  picture  of  what  that  target  is 
actually  doing.  So  a  threat  target  will  quite  often  switch  from  a  non-maneuvering  set  to 
the  maneuvering  set,  and  vice  versa,  to  avoid  being  tracked  all  the  time.  This  knowledge 
can  be  used  as  an  antecedent  for  situation  awareness. 

2.  The  frequency  of  appearance  of  such  type  of  target  based  on  some  a  priori  knowledge 
such  as  archival  radar  data.  Generally  threat  targets  axe  new  comparing  to  archival  radar 
data. 
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3.  The  importance  of  geolocation  of  this  target  based  on  the  geographical  information  sys¬ 
tems  (GISs).  Examples  of  important  geolocations  include  large  metroplex,  landmarks, 
military  bases,  airport,  etc.  Threats  happen  quite  often  in  such  areas. 

Of  the  above  three  antecedents,  the  frequency  of  appearance  of  such  type  of  target  is  weak 
stimulus  and  should  be  favored.  A  typical  rule  using  the  above  three  antecedents  can  be 


IF  the  number  of  switches  from  non-maneuvering  set  to  the  maneuvering  set  is  High ,  and  the 
frequency  of  appearance  of  such  target  is  Low ,  and  the  importance  of  geolocation  of  such 
type  of  target  is  High,  THEN  the  possibility  that  an  I&W  needs  to  be  issued  is  Very  Strong. 


The  linguistic  variables  used  to  represent  each  antecedent  are  divided  into  three  levels: 
Low,  Moderate,  and  High.  The  consequent  -  the  possibility  that  an  indication  and  warning 
(I&W)  needs  to  be  issued  -  is  divided  into  5  levels,  Very  Strong ,  Strong,  Medium,  Weak,  Very 
Weak.  So  we  need  to  set  up  33  =  27  (because  every  antecedent  has  3  fuzzy  sub-sets,  and  there 
are  3  antecedents)  rules  for  this  FLS.  Table  1  summarizes  the  fuzzy  rules  we  use  in  this  paper. 
We  show  these  MFs  in  Fig.  5. 

For  input  (xi,X2,xz),  the  output  is  computed  using 


S/=l  Mp1  {xl)fJ'p^{x2)f£p^{x3)cavg 
y  Xl,X2,X3  Y^lLl  MpJ  (*1  )Pp*  (X2)Mp3  (Xs) 


(15) 


By  repeating  these  calculations  for  Vxj  e  [0, 10],  we  obtain  a  hypersurface  y(x\,  £2,  £3).  This 
equation  represents  the  nonlinear  mapping  between  three  inputs  and  one  output  of  the  FLS. 
Since  it’s  a  4-D  surface  (x\,X2,x$ ,y),  it’s  impossible  to  be  plotted  visually. 

If  we  have  £3  =  8,  and  two  other  antecedents,  x\  and  £2  are  variables,  we  obtain  a  hyper¬ 
surface  y(x i,£2,8)  based  on  the  favor  weak  FLS,  as  plotted  in  Fig.  6(a).  In  contrast,  we  use 
a  type-1  FLS  where  the  antecedent  membership  functions  are  the  lower  membership  functions 
in  Fig.  5a  since  favor  weak  mechanism  is  not  used,  and  its  output  hypersurface  y(x i,£2,8)  is 
plotted  in  Fig.  6(b).  Observe  that  from  Fig.  6,  the  favor  weak  FLS  provides  a  higher  possibility 
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that  this  target  is  a  threat,  which  makes  sense  because  the  weak  factor,  frequency  of  appear¬ 
ance  of  such  type  of  target,  has  been  favored.  So  our  proposed  favor  weak  FLS  can  increase 
the  probability  of  threat  detection,  and  provides  timely  I&W. 

7  Conclusions 

Humans  use  multiple  sources  of  sensory  information  to  estimate  environmental  properties  and 
has  innate  ability  to  integrate  information  from  heterogeneous  data  sources.  There  is  consensus 
that  it  depends  on  the  brain  PFC.  The  PFC  has  top-down  control  (favor  weak)  and  rule-based 
mechanisms,  which  can  be  illustrated  using  the  Stroop  model.  In  this  paper,  we  proposed 
to  incorporate  the  favor  weak  mechanism  into  rule-based  fuzzy  logic  systems  (FLS)  via  using 
upper  and  lower  membership  functions.  The  inference  engine  of  favor  weak  fuzzzy  logic  system 
was  proposed  under  three  different  categories  based  on  fuzzifiers.  We  analyzed  that  the  favor 
weak  FLS  is  a  special  type-1  FLS  which  is  embeded  in  an  interval  type-2  FLS,  so  it’s  much 
simpler  in  computing  than  an  interval  type-2  FLS.  We  apply  the  favor  weak  FLS  to  situation 
understanding  based  on  heterogeneous  sensor  network,  and  it  shows  that  our  favor  weak  fuzzy 
logic  system  has  clear  advantage  comparing  to  the  type-1  FLS.  The  favor  weak  FLS  can  increase 
the  probability  of  threat  detection,  and  provides  timely  indication  &  warning  (I&W). 
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Table  1:  Fuzzy  rules  used  in  the  application.  Ante  1  is  the  number  of  switches  from  non¬ 
maneuvering  set  to  the  maneuvering  set  or  vice  versa;  Ante  2  the  frequency  of  appearance  of 
such  type  of  target;  Ante  3  is  the  importance  of  geolocation  of  this  target;  and  Consequent  is 
the  possibility  that  this  target  is  a  threat. 


Rule  # 

Ante  1 

Ante  2 

Ante  3 

Consequent 

i 

low 

low 

low 

Weak 

2 

low 

low 

moderate 

Medium 

3 

low 

low 

high 

Strong 

4 

low 

moderate 

low 

Very  Weak 

5 

low 

moderate 

moderate 

Weak 
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low 

moderate 

high 

Medium 
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low 

high 

low 

Very  Weak 
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low 

high 

moderate 

Weak 

9 

low 

high 

high 

Medium 

10 

moderate 

low 

low 

Medium 

11 

moderate 

low 

moderate 

Strong 

12 

moderate 

low 

high 

Very  Strong 

13 

moderate 

moderate 

low 

Weak 

14 

moderate 

moderate 

moderate 

Medium 

15 

moderate 

moderate 

high 

Strong 

16 

moderate 

high 

low 

Very  Weak 

17 

moderate 

high 

moderate 

Weak 

18 

moderate 

high 

high 

Medium 

19 

high 

low 

low 

Medium 

20 

high 

low 

moderate 

Strong 

21 

high 

low 

high 

Very  Strong 

22 

high 

moderate 

low 

Weak 

23 

high 

moderate 

moderate 

Medium 

24 

high 

moderate 

high 

Strong 

25 
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high 
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high 

high 

moderate 

Weak 

27 

high 

high 

high 

Moderate 

20 


Page  25  of  134 


List  of  Figures 


1  Schematic  diagram  of  some  of  the  extrinsic  and  intrinsic  connections  of  the  PFC. 

Most  connections  are  reciprocal;  the  exceptions  are  indicated  by  arrows.  The 
frontal  eye  field  (FEF)  has  variously  been  considered  either  adjacent  to,  or  part 

of,  the  PFC .  21 

2  Schematic  of  the  Stroop  model,  (a)  No  control.  Activation  of  conflicting  inputs 

in  the  two  pathways  produces  a  response  associated  with  the  word,  due  to  the 
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Figure  1:  Schematic  diagram  of  some  of  the  extrinsic  and  intrinsic  connections  of  the  PFC. 
Most  connections  are  reciprocal;  the  exceptions  are  indicated  by  arrows.  The  frontal  eye  field 
(FEF)  has  variously  been  considered  either  adjacent  to,  or  part  of,  the  PFC. 
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No  Control 


(b) 


Figure  2:  Schematic  of  the  Stroop  model,  (a)  No  control.  Activation  of  conflicting  inputs  in  the 
two  pathways  produces  a  response  associated  with  the  word,  due  to  the  stronger  connections  in 
the  word  reading  pathway,  (b)  Presentation  of  a  conflict  stimulus.  The  color  unit  is  activated 
(indicated  by  the  orange  fill),  representing  the  current  intent  to  name  the  color.  This  passes 
activation  to  the  intermediate  units  in  the  color  naming  pathway  (indicated  by  arrows),  which 
primes  those  units  (indicated  by  larger  size),  and  biases  processing  in  favor  of  activity  flowing 
along  this  pathway.  This  biasing  effect  favors  activation  of  the  response  unit  corresponding  to 
the  color  input,  even  though  the  connection  weights  in  this  pathway  are  weaker  than  in  the 
word  pathway. 
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Figure  3:  The  type-2  MFs  for  Example  1.  The  thick  solid  lines  denote  upper  MFs,  and  the 
thick  dashed  lines  denote  lower  MFs.  The  shaded  regions  are  the  footprints  of  uncertainty  for 
interval  secondaries.  The  center  of  the  Gaussian  MFs  is  5,  and  the  variance  varies  from  1.0  to 
2.0. 
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Figure  4:  (a)  Pictorial  representation  of  a  type-2  fuzzy  set.  The  secondary  memberships  are 
shown  in  (b)  to  represent  the  confidence  of  the  primary  memberships. 
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Figure  5:  The  MFs  used  to  represent  the  linguistic  labels,  (a)  MFs  for  antecedents,  and  (b) 
MFs  for  consequent. 
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Figure  6:  The  threat  assessment  surface  for  different  FLSs  when  X3  =  8.  (a)  Favor  weak  FLS, 
and  (b)  Traditional  type-1  FLS. 
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Abstract 

Inspired  by  human’s  innate  ability  to  process  and  integrate  information  from  disparate, 
network-based  sources,  we  apply  human-inspired  information  integration  mechanisms  to  tar¬ 
get  detection  in  cognitive  radar  sensor  network.  Humans’  information  integration  mechanisms 
have  been  modelled  using  maximum-likelihood  estimation  (MLE)  or  soft-max  approaches.  In 
this  paper,  we  apply  these  two  algorithms  to  cognitive  radar  sensor  networks  target  detection. 
Discrete-cosine-transform  (DCT)  is  used  to  process  the  integrated  data  from  MLE  or  soft- 
max.  We  apply  fuzzy  logic  system  (FLS)  to  automatic  target  detection  based  on  the  AC  power 
values  from  DCT.  Simulation  results  show  that  our  MLE-DCT-FLS  and  soft-max-DCT-FLS 
approaches  perform  very  well  in  the  radar  sensor  network  target  detection,  whereas  the  existing 
2-D  construction  algorithm  doesn’t  work  in  this  study. 

Index  Terms  :  Human-inspired,  cognitive  radar  sensor  networks,  fuzzy  logic  systems,  auto¬ 
matic  target  recognition,  maximum-likelihood  estimation. 


1 


Page  33  of  134 


1  Introduction  and  Motivation 


A  cognitive  network  is  one  that  is  aware  of  changes  in  user  needs  and  its  environment,  adapts  its 
behavior  to  those  changes,  learns  from  its  adaptations,  and  exploits  knowledge  to  improve  its  future 
behavior.  A  cognitive  radar  sensor  network  consists  of  multiple  networked  radar  sensors  and  radar 
sensors  sense  and  communicate  with  each  other  collaboratively  to  complete  a  mission.  In  real  world, 
cognitive  radar  sensor  network  information  integration  is  necessary  in  different  applications.  For 
example,  in  an  emergency  natural  disaster  scenario,  such  as  Utah  Mine  Collapse  in  August  2007  or 
West  Virginia  Sago  mine  disaster  in  January  2006,  cognitive  radar  sensor  network-based  information 
integration  for  first  responders  is  critical  for  search  and  rescue.  Danger  may  appear  anywhere  at 
any  time,  therefore,  first  responders  must  monitor  a  large  area  continuously  in  order  to  identify 
potential  danger  and  take  actions.  Due  to  the  dynamic  and  complex  nature  of  natural  disaster,  some 
buried/foleage  victims  may  not  be  found  with  image/video  sensors,  and  UWB  radar  sensors  are 
needed  for  penetrating  the  ground  or  sense-through-wall.  Unfortunately,  the  radar  data  acquired 
are  often  limited  and  noisy.  Unlike  medical  imaging  or  synthetic  aperture  radar  imaging  where 
abundance  of  data  is  generally  available  through  multiple  looks  and  where  processing  time  may  not 
be  crucial,  practical  cognitive  radar  sensor  networks  are  typically  the  opposite:  availability  of  data 
is  limited  and  required  processing  time  is  short.  This  need  is  also  motivated  by  the  fact  that  humans 
display  a  remarkable  capability  to  quickly  perform  target  recognition  despite  noisy  sensory  signals 
and  conflicting  inputs.  Humans  are  adept  at  network  visualization,  and  at  understanding  subtle 
implications  among  the  network  connections.  To  date,  however,  human’s  innate  ability  to  process 
and  integrate  information  from  disparate,  network-based  sources  for  situational  understanding  has 
not  translated  well  to  automated  systems.  In  this  paper,  we  apply  human  information  integration 
mechanisms  to  information  fusion  in  cognitive  radar  sensor  network. 
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2  Sense-through-Foliage  Radar  Sensor  Networks  Data  Measure¬ 
ment  and  Collection 

Our  work  is  based  on  the  sense-through-foliage  UWB  radar  sensor  networks.  The  foliage  experiment 
was  constructed  on  a  seven-ton  man  lift,  which  had  a  total  lifting  capacity  of  450  kg.  The  limit  of 
the  lifting  capacity  was  reached  during  the  experiment  as  essentially  the  entire  measuring  apparatus 
was  placed  on  the  lift,  (as  shown  in  Fig.  1).  The  principle  pieces  of  equipment  secured  on  the  lift  are: 
Barth  pulser,  Tektronix  model  7704  B  oscilloscope,  dual  antenna  mounting  stand,  two  antennas, 
rack  system,  IBM  laptop,  HP  signal  Generator,  Custom  RF  switch  and  power  supply  and  Weather 
shield  (small  hut).  The  target  is  a  trihedral  reflector  (as  shown  in  Fig.  2).  Throughout  this  work, 
a  Barth  pulse  source  (Barth  Electronics,  Inc.  model  732  GL)  was  used.  The  pulse  generator  uses 
a  coaxial  reed  switch  to  discharge  a  charge  line  for  a  very  fast  rise  time  pulse  outputs.  The  model 
732  pulse  generator  provides  pulses  of  less  than  50  picoseconds  (ps)  rise  time,  with  amplitude  from 
150  V  to  greater  than  2  KV  into  any  load  impedance  through  a  50  ohm  coaxial  line.  The  generator 
is  capable  of  producing  pulses  with  a  minimum  width  of  750  ps  and  a  maximum  of  1  microsecond. 
This  output  pulse  width  is  determined  by  charge  line  length  for  rectangular  pulses,  or  by  capacitors 
for  1/e  decay  pulses. 

For  the  data  we  used  in  this  paper,  each  sample  is  spaced  at  50  picosecond  interval,  and  16,000 
samples  were  collected  for  each  collection  for  a  total  time  duration  of  0.8  microseconds  at  a  rate  of 
approximately  20  Hz.  We  plot  the  transmitted  pulse  (one  realization)  in  Fig.  3a)  and  the  received 
echos  in  one  collection  in  Fig.  3b  (averaged  over  35  pulses).  The  data  collections  were  extensive. 
20  different  positions  were  used,  and  35  collections  were  performed  at  each  position  using  UWB 
radar  sensor  networks. 
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3  Human  Information  Integration  Mechanisms 

Recently,  a  maximum-likelihood  estimation  (MLE)  approach  was  proposed  for  multi-sensory  data 
fusion  in  human  [4].  In  the  MLE  approach  [4],  sensory  estimates  of  an  environmental  property  can 
be  represented  by  Sj  =  /;(S)  where  S  is  the  physical  property  being  estimated,  /  is  the  operation 
the  nervous  system  performs  to  derive  the  estimate,  and  S  is  the  perceptual  estimate.  Sensory 
estimates  are  subject  to  two  types  of  error:  random  measurement  error  and  bias.  Thus,  estimates 
of  the  same  object  property  from  different  cues  usually  differ.  To  reconcile  the  discrepancy,  the 
nervous  system  must  either  combine  estimates  or  choose  one,  thereby  ignoring  the  other  cues. 
Assuming  that  each  single-cue  estimate  is  unbiased  but  corrupted  by  independent  Gaussian  noise, 
the  statistically  optimal  strategy  for  cue  combination  is  a  weighted  average  [4] 

M 

Sc  =  ^  ^  WjSj  (1) 

i=i 

where  Wi  =  and  is  the  weight  given  to  the  zth  single-cue  estimate,  of  is  that  estimates 

variance,  and  M  is  the  total  number  of  cues.  Combining  estimates  by  this  MLE  rule  yields  the 
least  variable  estimate  of  S  and  thus  more  precise  estimates  of  object  properties. 

Besides,  some  other  summation  rules  have  been  proposed  in  perception  and  cognition  such  as 
soft-max  rule:  y  =  (YliL l  x?)”  [3]  where  X{  denotes  the  input  from  an  input  source  i,  and  M  is  the 
total  number  of  sources.  In  this  paper,  we  will  apply  MLE  and  soft-max  human  brain  information 
integration  mechanisms  to  cognitive  radar  sensor  network  information  integration. 


4  Human-Inspired  Sense-through-Foliage  Target  Detection 

In  Figs.  4a  and  4b,  we  plot  two  collections  of  UWB  radars.  Fig.  4a  has  no  target  on  range,  and  Fig. 
4b  has  target  at  samples  around  13,900.  We  plot  the  echo  differences  between  Figs.  4a  and  4b  in 
Fig.  4c.  However,  it  is  impossible  to  identify  whether  there  is  any  target  and  where  there  is  target 
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based  on  Fig.  4c.  Since  significant  pulse-to-pulse  variability  exists  in  the  echos,  this  motivate  us  to 
explore  the  spatial  and  time  diversity  using  Radar  Sensor  Networks  (RSN). 

In  RSN,  the  radar  sensors  are  networked  together  in  an  ad  hoc  fashion.  They  do  not  rely  on  a 
preexisting  fixed  infrastructure,  such  as  a  wireline  backbone  network  or  a  base  station.  They  are  self¬ 
organizing  entities  that  are  deployed  on  demand  in  support  of  various  events  surveillance,  battlefield, 
disaster  relief,  search  and  rescue,  etc.  Scalability  concern  suggests  a  hierarchical  organization 
of  radar  sensor  networks  with  the  lowest  level  in  the  hierarchy  being  a  cluster.  As  argued  in 
[7]  [6]  [5]  [9],  in  addition  to  helping  with  scalability  and  robustness,  aggregating  sensor  nodes  into 
clusters  has  additional  benefits: 

1.  conserving  radio  resources  such  as  bandwidth; 

2.  promoting  spatial  code  reuse  and  frequency  reuse; 

3.  simplifying  the  topology,  e.g.,  when  a  mobile  radar  changes  its  location,  it  is  sufficient  for 
only  the  nodes  in  attended  clusters  to  update  their  topology  information; 

4.  reducing  the  generation  and  propagation  of  routing  information;  and, 

5.  concealing  the  details  of  global  network  topology  from  individual  nodes. 

In  RSN,  each  radar  can  provide  their  pulse  parameters  such  as  timing  to  their  clusterhead  radar, 
and  the  clusterhead  radar  can  combine  the  echos  (RF  returns)  from  the  target  and  clutter.  In  this 
paper,  we  propose  a  RAKE  structure  for  combining  echos,  as  illustrated  by  Fig.  5.  The  integration 
means  time-average  for  a  sample  duration  T  and  it’s  for  general  case  when  the  echos  are  not  in 
discrete  values.  It  is  quite  often  assumed  that  the  radar  sensor  platform  will  have  access  to  Global 
Positioning  Service  (GPS)  and  Inertial  Navigation  Unit  (INU)  timing  and  navigation  data  [1].  In 
this  paper,  we  assume  the  radar  sensors  are  synchronized  in  RSN.  In  Fig.  5,  the  echo,  i.e.,  RF 
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response  by  the  pulse  of  each  cluster-member  sensor,  will  be  combined  by  the  clusterhead  using  a 
weighted  average,  and  the  weight  is  determined  by  the  two  human-inspired  mechanisms. 

We  applied  the  human-inspired  MLE  algorithm  to  combine  the  sensed  echo  collection  from 
M  =  30  UWB  radars,  and  then  the  combined  data  are  processed  using  discrete-cosine  transform 
(DCT)  to  obtain  the  AC  values.  Based  on  our  experiences,  echo  with  a  target  generally  has 
high  and  nonfluctuating  AC  values  and  the  AC  values  can  be  obtained  using  DCT.  We  plot  the 
power  of  AC  values  in  Figs.  6a  and  6b  using  MLE  and  DCT  algorithms  for  the  two  cases  (with 
target  and  without  target)  respectively.  Observe  that  in  Fig.  6b,  the  power  of  AC  values  (around 
sample  13,900)  where  the  target  is  located  is  non-fluctuating  (somehow  monotonically  increase  then 
decrease).  Although  some  other  samples  also  have  very  high  AC  power  values,  it  is  very  clear  that 
they  are  quite  fluctuating  and  the  power  of  AC  values  behaves  like  random  noise  because  generally 
the  clutter  has  Gaussian  distribution  in  the  frequency  domain. 

Similarly,  we  applied  the  soft-max  algorithm  (n  =  2)  to  combine  the  sensed  echo  collection  from 
M  =  30  UWB  radars,  and  then  used  DCT  to  obtain  the  AC  values.  We  plot  the  power  of  AC  values 
in  Figs.  6a  and  6b  using  soft-max  and  DCT  algorithms  for  the  two  cases  (with  target  and  without 
target)  respectively.  Observe  that  in  Fig.  7b,  the  power  of  AC  values  (around  sample  13,900)  where 
the  target  is  located  is  non-fluctuating  (somehow  monotonically  increase  then  decrease). 

We  made  the  above  observations.  However,  in  real  world  application,  automatic  target  detection 
is  necessary  to  ensure  that  our  algorithms  could  be  performed  in  real  time.  In  Section  5,  we  apply 
fuzzy  logic  systems  to  automatic  target  detection  based  on  the  power  of  AC  values  (obtained  via 
MLE-DCT  or  soft-max-DCT). 

We  compared  our  approaches  to  the  scheme  proposed  in  [10].  In  [10],  2-D  image  was  created 
via  adding  voltages  with  the  appropriate  time  offset.  In  Figs.  8a  and  8b,  we  plot  the  2-D  image 
created  based  on  the  above  two  data  sets  (from  samples  13,800  to  14,200).  The  sensed  data  from 
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30  radars  are  averaged  first,  then  plotted  in  2-D  [10].  However,  it’s  not  clear  which  image  shows 
there  is  target  on  range. 


5  Fuzzy  Logic  System  for  Automatic  Target  Detection 

5.1  Overview  of  Fuzzy  Logic  Systems 

Figure  9  shows  the  structure  of  a  fuzzy  logic  system  (FLS)  [8].  When  an  input  is  applied  to  a 
FLS,  the  inference  engine  computes  the  output  set  corresponding  to  each  rule.  The  defuzzifer 
then  computes  a  crisp  output  from  these  rule  output  sets.  Consider  a  p-input  1-output  FLS,  using 
singleton  fuzzification,  center- of- sets  defuzzification  [8]  and  “IF-THEN”  rules  of  the  form 

Rl  :  IF  x\  is  F,  and  X2  is  Fj  and  •  •  •  and  xp  is  Flp,  THEN  y  is  G*. 

Assuming  singleton  fuzzification,  when  an  input  x'  =  {x'j, . . .  ,x'p }  is  applied,  the  degree  of  firing 
corresponding  to  the  Zth  rule  is  computed  as 


Mf<  (*'i)  *  Mf' (4)  *  ’  •  •  *  Mf| ,(xp)  =  VLiWiA) 


(2) 


where  *  and  T  both  indicate  the  chosen  f-norm.  There  are  many  kinds  of  defuzzifiers.  In  this  paper, 
we  focus,  for  illustrative  purposes,  on  the  center-of-sets  defuzzifier  [8].  It  computes  a  crisp  output 
for  the  FLS  by  first  computing  the  centroid,  cGi,  of  every  consequent  set  G/,  and,  then  computing  a 
weighted  average  of  these  centroids.  The  weight  corresponding  to  the  Zth  rule  consequent  centroid 
is  the  degree  of  firing  associated  with  the  Zth  rule,  l^=1nFi(x'i),  so  that 


,  #x  _ 


(3) 


where  M  is  the  number  of  rules  in  the  FLS.  In  this  paper,  we  design  a  FLS  for  automatic  target 
recognition  based  on  the  AC  values  obtained  using  MLE-DCT  or  soft-max-DCT. 
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5.2  FLS  for  Automatic  Target  Detection 

Observe  that  in  Figs.  6  and  7,  the  power  of  AC  values  are  quite  fluctuating  and  have  lots  of 
uncertainties.  FLS  is  well  known  to  handle  the  uncertainties.  For  convenience  in  describing  the 
FLS  design  for  Automatic  Target  Detection  (ATD),  we  first  give  the  definition  of  footprint  of 
uncertainty  of  AC  power  values  and  region  of  interest  in  the  footprint  of  uncertainty. 

Definition  1  (Footprint  of  Uncertainty)  Uncertainty  in  the  AC  power  values  and  time  index 
consists  of  a  bounded  region,  that  we  call  the  footprint  of  uncertainty  of  AC  power  values.  It  is  the 
union  of  all  AC  power  values. 

Definition  2  (Region  of  Interest  (Rol))  An  Rol  in  the  footprint  of  uncertainty  is  a  contour 
consisting  of  a  large  number  (greater  than  50)  of  AC  power  values  where  AC  power  values  increase 
then  decrease. 

Definition  3  (Fluctuating  Point  in  Rol)  P{i)  is  called  a  fluctuating  point  in  the  Rol  if  P(i  - 
1),  P(i),  P[i  +  1)  are  non-monotonically  increasing  or  decreasing. 

Our  FLS  for  automatic  target  detection  will  classify  each  ROI  (with  target  or  no  target)  based 
on  two  antecedents:  the  centroid  of  the  ROI  and  the  number  of  fluctuating  points  in  the  ROI. 
The  linguistic  variables  used  to  represent  these  two  antecedents  were  divided  into  three  levels:  low, 
moderate ,  and  high.  The  consequent  —  the  possibility  that  there  is  a  target  at  this  Rol  —  was  divided 
into  5  levels,  Very  Strong,  Strong,  Medium,  Weak,  Very  Weak.  We  used  trapezoidal  membership 
functions  (MFs)  to  represent  low,  high,  very  strong,  and  very  weak;  and  triangle  MFs  to  represent 
moderate,  strong,  medium,  and  weak.  All  inputs  to  the  antecedents  are  normalized  to  0-10. 

Based  on  the  fact  the  AC  power  value  of  target  is  non-fluctuating  (somehow  monotonically 
increase  then  decrease),  and  the  AC  power  value  of  clutter  behaves  like  random  noise  because 
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generally  the  clutter  has  Gaussian  distribution  in  the  frequency  domain,  we  design  a  fuzzy  logic 
system  using  rules  such  as: 

Rl  :  IF  centroid  of  a  Rol  (#i)  is  F*,  and  the  number  of  fluctuating  points  in  the  ROI  (£2)  is  F^, 
THEN  the  possibility  that  there  is  a  target  at  this  Rol  ( y )  is  Gl. 


where  l  =  1, . . . ,  9.  We  summarize  all  the  rules  in  Table  1.  For  every  input  (xi,  £2)?  the  output  is 


computed  using 


y(x  1,2:2) 


Ef=l  % }  (Xi)^(X2)(jvg 

EiLiMFifoW’te) 


(4) 


We  ran  simulations  to  1000  collections  in  the  real  world  sense-through-foliage  experiment,  and 


found  that  our  FLS  performs  very  well  in  the  automatic  target  detection  based  on  the  AC  power 


values  obtained  from  MLE-DCT  or  soft-max-DCT,  and  achieve  probability  of  detection  pd  =  100% 
and  false  alarm  rate  pja  =  0. 


6  Conclusions 

Inspired  by  human’s  innate  ability  to  process  and  integrate  information  from  disparate,  network- 
based  sources,  we  applied  human-inspired  information  integration  mechanisms  to  target  detection  in 
cognitive  radar  sensor  network.  Humans’  information  integration  mechanisms  have  been  modelled 
using  maximum-likelihood  estimation  (MLE)  or  soft-max  approaches.  In  this  paper,  we  applied 
these  two  algorithms  to  cognitive  radar  sensor  networks  target  detection.  Discrete-cosine-transform 
(DCT)  was  used  to  process  the  integrated  data  from  MLE  or  soft-max.  We  applied  fuzzy  logic 
system  (FLS)  to  automatic  target  detection  based  on  the  AC  power  values  from  DCT.  Simulation 
results  showed  that  our  MLE-DCT-FLS  and  soft-max-DCT-FLS  approaches  performed  very  well  in 
the  radar  sensor  network  target  detection,  whereas  the  existing  2-D  construction  algorithm  couldn’t 
work  in  this  study. 
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Table  1:  The  rules  for  target  detection.  Antecedent  1  is  centroid  of  a  Rol ,  Antecedent  2  is  the 
number  of  fluctuating  points  in  the  ROI,  and  Consequent  is  the  possibility  that  there  is  a  target  at 
this  Rol. 


Rule  # 

Antecedent  1 

Antecedent  2 

Consequent 

1 

low 

low 

medium 

2 

low 

moderate 

weak 

3 

low 

high 

very  weak 

4 

moderate 

low 

strong 

5 

moderate 

moderate 

medium 

6 

moderate 

high 

weak 

7 

high 

low 

very  strong 

8 

high 

moderate 

strong 

9 

high 

high 

medium 
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1  This  figure  shows  the  lift  with  the  experiment.  The  antennas  are  at  the  far  end  of  the 
lift  from  the  viewer  under  the  roof  that  was  built  to  shield  the  equipment  from  the 
elements.  This  picture  was  taken  in  September  with  the  foliage  largely  still  present. 

The  cables  coming  from  the  lift  are  a  ground  cable  to  an  earth  ground  and  one  of  4 
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Figure  1:  This  figure  shows  the  lift  with  the  experiment.  The  antennas  are  at  the  far  end  of  the 
lift  from  the  viewer  under  the  roof  that  was  built  to  shield  the  equipment  from  the  elements.  This 
picture  was  taken  in  September  with  the  foliage  largely  still  present.  The  cables  coming  from  the 
lift  are  a  ground  cable  to  an  earth  ground  and  one  of  4  tethers  used  in  windy  conditions. 
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Figure  2:  The  target  (a  trihedral  reflector)  is  shown  on  the  stand  at  300  feet  from  the  lift. 
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Figure  3:  Transmitted  pulse  and  received  echos  in  one  experiment,  (a)  Transmitted  pulse,  (b) 
Received  echos. 
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Figure  4:  Measurement  with  35  pulses  average,  (a)  Expanded  view  of  traces  (no  target)  from 
sample  13,001  to  15,000.  (b)  Expanded  view  of  tl&ces  (with  target)  from  samples  13,001  to  15,000. 
(c)  The  differences  between  (a)  and  (b).  Page  50  of  134 


Figure  5:  Echo  combining  by  clusterhead  in  RSN. 
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Figure  6:  Power  of  AC  values  using  MLE-based  information  integration  and  DCT.  (a)  No  target 
(b)  With  target  in  the  field. 
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(b) 


Figure  7:  Power  of  AC  values  using  soft-max  based  information  integration  and  DCT.  (a)  No  target 
(b)  With  target  in  the  field. 
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Figure  8:  2-D  image  created  via  adding  voltages  with  the  appropriate  time  offset,  (a)  No  target 
(b)  With  target  in  the  field. 
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FUZZY  LOGIC  SYSTEM 


Figure  9:  The  structure  of  a  fuzzy  logic  system. 
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Wireless  Channel  Modeling  in  Foliage  Environment: 


UWB  versus  Narrowband 
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Abstract 

In  this  paper,  we  study  the  wireless  channel  modeling  in  foliage  environment,  a  rich  scattering 
and  time- varying  environment,  based  on  extensive  data  collected  using  UWB  and  narrowband 
(200MHz  and  400MHz)  radars.  We  apply  two  approaches  to  the  wireless  channel  modeling: 
Saleh  and  Valenzuela  (S-V)  method  for  UWB  channel  modeling  and  CLEAN  method  for  nar¬ 
rowband  and  UWB  channel  modeling.  We  validated  that  UWB  echo  signals  (within  a  burst) 
don’t  hold  self-similarity,  which  means  the  future  signals  can’t  be  forecasted  based  on  the  re¬ 
ceived  signals  and  channel  modeling  is  necessary  from  statistical  point  of  view.  Based  on  the 
S-V  method  for  UWB  channel  modeling,  in  foliage  UWB  channel,  the  multi-path  contributions 
arrive  at  the  receiver  are  grouped  into  clusters.  The  time  of  arrival  of  clusters  can  be  modeled 
as  a  Poisson  arrival  process,  while  within  each  cluster,  subsequent  multipath  contributions  or 
rays  also  arrive  according  to  a  Poisson  process.  At  different  field  (near  field,  medium  field,  and 
far  field),  we  observe  that  the  Poisson  process  parameters  are  quite  different.  We  also  observe 
that  the  amplitude  of  channel  coefficient  at  each  path  follows  Rician  distribution  for  medium 
and  far  field,  and  it’s  non-stationary  for  paths  from  near  field  (one  of  two  Rician  distributions), 
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and  these  observations  are  quite  different  with  the  IEEE  indoor  UWB  channel  model  and  S-V 
indoor  channel  model.  Based  on  the  CLEAN  method,  the  narrowband  (200MHz  and  400MHz) 
and  UWB  channel  impulse  responses  have  many  similarities:  both  can  be  modeled  as  linear 
time- variant  filter  channel. 


1  Introduction  and  Motivation 

The  true  challenge  for  new  communication  technologies  is  to  “make  the  thing  work”  in  real-world 
wireless  channels.  System  designers  classically  focus  on  the  impact  of  the  radio  channel  on  the 
received  signals  and  use  propagation  models  for  testing  and  evaluation  of  receiver  designs  and 
transmission  schemes.  Yet,  the  needs  for  such  models  evolve  as  new  applications  emerge  with 
different  bandwidths,  terminal  mobility,  higher  carrier  frequencies,  new  antennas,  and  so  forth. 
Furthermore,  channel  characterization  also  yields  the  fundamental  ties  to  classical  electromagnetics 
and  physics,  as  well  as  the  answers  to  some  crucial  questions  in  communication  and  information 
theory  [8].  While  many  efforts  have  been  spent  on  indoor  channel  modeling  as  well  as  outdoor  open 
space  (or  less  scattering)  environment,  wireless  channel  modeling  in  foliage  (forest)  environment 
has  not  been  studied.  In  Vietnam,  Bosnia,  Kosovo,  and  in  the  jungles  of  Columbia,  air  to  ground 
communnications  have  been  thwarded  by  the  presence  of  foliage  that  protects  and  hides  the  enemy. 
In  this  paper,  we  study  UWB  and  narrowband  channel  modeling  in  foliage  environment. 

In  July  2003,  the  Channel  Modeling  sub-committee  of  study  group  IEEE  802.15.SG3a  published 
the  final  report  regarding  the  UWB  indoor  multipath  channel  model  [6].  It  is  a  modified  version 
of  the  indoor  Saleh  and  Valenzuela  (S-V)  channel  model  [10].  The  S-V  model  was  developped 
for  NLOS  channel,  and  it  has  also  been  applied  to  LOS  channels  where  it  is  perhaps  less  valid, 
unless  LOS  components  are  specifically  added  [9]  The  IEEE  suggested  an  initial  set  of  values  for 
the  indoor  UWB  channel  model  which  has  range  less  than  10  meters.  However,  lots  of  applications 
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of  UWB  axe  for  outdoor  activities  such  as  sense-through- foliage  target  detection.  Forests  favor 
asymmetric  threats  because  the  warfighter  has  a  limited  sensing  capability.  Forests  provide  excellent 
concealment  from  observation,  ambush,  and  escape,  as  well  as  provide  secure  bases  for  enemy 
Command  &  Control  (C2),  weapons  caches,  and  Improvised  Explosive  Device  (IED)/  Weapon  of 
Mass  Destruction  (WMD)  assembly.  These  have  become  “the  high  ground”  in  fourth-generation 
warfare,  providing  a  significant  strategic  advantage.  Unfortunately,  no  work  has  been  done  on  the 
outdoor  UWB  channel  modeling. 

In  this  paper,  we  will  model  the  UWB  and  narrowband  channels  using  UWB  and  narrowband 
radars  in  foliage  environment  which  is  a  rich  scattering  and  time- varying  environment.  UWB  radar 
emissions  are  at  a  relatively  low  frequency-typically  between  100  MHz  and  3  GHz.  Additionally, 
the  fractional  bandwidth  of  the  signal  is  very  large  (greater  than  0.2).  Such  radar  sensor  has 
exceptional  range  resolution  that  also  has  an  ability  to  penetrate  many  common  materials  (e.g., 
walls).  Law  enforcement  personnel  have  used  UWB  ground  penetrating  radars  (GPRs)  for  at  least 
a  decade.  Like  the  GPR,  sense-through- foliage  radar  takes  advantage  of  UWB’s  very  fine  resolution 
(time  gating)  and  low  frequency  of  operation. 

The  rest  of  this  paper  is  organized  as  follows.  In  Section  2,  we  summarize  the  measurement  and 
collection  of  data  we  used  in  this  paper.  In  Section  3,  we  demonstrate  that  the  UWB  reflected  signal 
in  foliage  environment  does  not  hold  self-similarity,  and  validate  that  outdoor  channel  modeling 
is  necessary.  In  Section  4,  we  present  our  outdoor  UWB  channel  model  in  rich  scattering  and 
time- varying  environment.  In  Section  5,  we  study  the  channel  impulse  response  based  on  CLEAN 
method  for  narrowband  and  UWB  channels.  We  conclude  this  paper  in  Section  6. 
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2  Experiment  Setup  and  Data  Collection 


Our  work  is  based  on  the  data  collected  in  UWB  radar-based  sense-through-foliage  experiment  in 
late  summer  and  fall.  Late  summer  foliage,  because  of  the  limited  rainfall,  involved  foliage  with 
decreased  water  content.  Late  fall  and  winter  measurements  involved  largely  defoliated  but  dense 
forest,  so  it’s  a  rich  scattering  environment.  Because  of  wind  or  different  temperature  in  dense 
forest,  it’s  also  a  time- varying  environment.  The  UWB  radar-based  experiment  was  constructed  on 
a  seven- ton  man  lift,  which  had  a  total  lifting  capacity  of  450  kg.  The  limit  of  the  lifting  capacity 
was  reached  during  the  experiment  as  essentially  the  entire  measuring  apparatus  was  placed  on  the 
lift  (as  shown  in  Fig.  1).  The  principle  pieces  of  equipment  secured  on  the  lift  are:  Barth  pulser, 
Tektronix  model  7704  B  oscilloscope,  dual  antenna  mounting  stand,  two  antennas,  rack  system, 
IBM  laptop,  HP  signal  Generator,  Custom  RF  switch  and  power  supply  and  Weather  shield  (small 
hut).  Throughout  this  work,  a  Barth  pulse  source  (Barth  Electronics,  Inc.  model  732  GL)  was 
used.  The  pulse  generator  uses  a  coaxial  reed  switch  to  discharge  a  charge  line  for  a  very  fast  rise 
time  pulse  outputs.  The  model  732  pulse  generator  provides  pulses  of  less  than  50  picoseconds  (ps) 
rise  time,  with  amplitude  from  150  V  to  greater  than  2  KV  into  any  load  impedance  through  a  50 
ohm  coaxial  line.  The  generator  is  capable  of  producing  pulses  with  a  minimum  width  of  750  ps 
and  a  maximum  of  1  microsecond.  This  output  pulse  width  is  determined  by  charge  line  length  for 
rectangular  pulses,  or  by  capacitors  for  1/e  decay  pulses. 

For  the  UWB  data  we  used  in  this  paper,  each  sample  is  spaced  at  50  picosecond  interval,  and 
16,000  samples  were  collected  for  each  collection  for  a  total  time  duration  of  0.8  microseconds  at  a 
rate  of  approximately  20  Hz.  The  Barth  pulse  source  was  operated  at  low  amplitude  and  35  pulses 
reflected  signal  were  averaged  for  each  collection.  Significant  pulse-to-pulse  variability  was  noted 
for  these  collections.  We  plot  the  transmitted  pulse  (one  realization)  in  Fig.  2a)  and  the  received 
echos  in  one  collection  in  Fig.  2b  (averaged  over  35  pulses). 
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For  comparison,  we  also  studied  narrowband  (200MHz  and  400MHz)  radar  signal  propagation. 
Henry  Radio  preamplifiers  were  model  50B-200  and  50B-400  for  200MHz  and  400MHz  respectively. 
The  Henry  Radio  power  amplifiers  were  TEMPO-2002A  and  TEMPO-2400A  for  1  KW  pulsed  at 
200MHz  and  400MHz  respectively.  The  amplifier  was  speced  at  a  minimum  bandwidth  of  2  MHz 
around  it  center  frequency.  The  source  for  200MHz  and  400MHz  narrow  band  wave  signal  was  an 
Agilent  8648A  signal  generator.  For  the  data  we  used  in  this  paper,  each  sample  is  spaced  at  50 
picosecond  interval,  and  16,000  samples  were  collected  for  each  collection  for  a  total  time  duration 
of  0.8  microseconds  at  a  rate  of  approximately  20  Hz.  Fig.  3a  shows  the  transmitted  signal  and 
Fig.  3b  shows  the  received  echos  (averaged  over  35  pulses)  for  200MHz  narrowband  radar.  Fig.  4a 
shows  the  transmitted  signal  and  Fig.  4b  shows  the  received  echos  (averaged  over  35  pulses)  for 
400MHz  narrowband  radar.  The  data  collections  were  extensive.  20  different  positions  were  used, 
and  35  collections  were  performed  at  each  position  for  UWB,  200MHZ,  and  400MHz  radars. 

3  Self-Similarity  Properties  of  UWB  Reflected  Signals 

It  has  been  observed  that  ethernet  video/ voice/data  traffic  have  self-similarity  [7]  [4]  [13].  According 
to  Stallings  [12],  “Self-similarity  is  such  an  important  concept  that,  in  a  way,  it  is  surprising  that 
only  recently  has  it  been  applied  to  data  communications  traffic  analysis.”,  and  “Since  1993,  a 
number  of  studies  reported  in  the  literature  have  documented  that  the  pattern  of  data  traffic 
is  well  modeled  by  self-similar  processes  in  a  wide  variety  of  real-world  networking  situations.” 
Such  self-similarity  is  quite  common  in  both  natural  and  human-made  phenomena  [12]  such  as  the 
distribution  of  earthquakes,  ocean  waves,  fluctuation  of  the  stock  market.  But  the  self-similarity 
of  UWB  signals  has  not  been  studied. 

For  a  detailed  discussion  on  self-similarity  in  time-series,  see  [13]  [12].  Here  we  briefly  present 
its  definition  [3].  Given  a  zero-mean,  stationary  time-series  X  =  ( Xt\t  =  1,2,3,  •••),  we  define 
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the  ra-aggregated  series  X =  (X^]k  =  1,2, 3,  •••)  by  summing  the  original  series  X  over 
nonoverlapping  blocks  of  size  m.  Then  it’s  said  that  X  is  H- self-similar,  if,  for  all  positive  m,  X^m>j 
has  the  same  distribution  as  X  rescaled  by  mH .  That  is, 

tm 

Xt  =  m~H  Xi  VmzN  (1) 

i=(t—  l)m+l 

If  X  is  i/-self-similar,  it  has  the  same  autocorrelation  function  r(k)  =  E[(Xt  —  p)(Xt+k  —  /i)]/cr2  as 
the  series  X ^  for  all  m,  which  means  that  the  series  is  distributionally  self-similar:  the  distribution 
of  the  aggregated  series  is  the  same  as  that  of  the  original. 

Self-similar  processes  can  show  long-range  dependence.  A  process  with  long-range  dependence 
has  an  autocorrelation  function  r(k)  ~  k~&  as  k  — ►  oo,  where  0  <  (3  <  1.  The  degree  of  self¬ 
similarity  can  be  expressed  using  Hurst  parameter  H  =  1  -  (3/2.  For  self-similar  series  with 
long-range  dependence,  1/2  <  H  <  1.  As  H  — ►  1,  the  degree  of  both  self-similarity  and  long-range 
depence  increases. 

One  method  that  has  been  widely  used  to  verify  self-similarity  is  the  variance-time  plot ,  which 
relies  on  the  slowly  decaying  variance  of  a  self-similar  series.  The  variance  of  is  plotted 

against  m  on  a  log-log  plot,  and  a  straight  line  with  slope  (— (3)  greater  than  -1  is  indicative  of 
self-similarity,  and  the  parameter  H  is  given  by  H  =  1  -(3/2.  We  use  this  method  in  this  paper.  In 
Fig.  5,  we  plot  the  variance  of  X against  m  on  a  log-log  plot  for  10  different  UWB  data  collections. 
From  this  figure,  it’s  very  clear  that  the  UWB  signal  does  not  has  self-similarity  because  its 
trace  has  slope  lower  than  —1.  This  conclusion  means  that  we  can’t  use  current  received  signals 
to  forecast  future  reflected  signals  within  one  collection,  so  channel  modeling  is  very  important  to 
UWB  outdoor  channel  because  the  charateristics  of  the  future  reflected  signal  could  be  known  in 
advance  if  its  channel  can  be  modelled. 
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4  UWB  Channel  Modeling  Based  on  S-V  Approaches 


4.1  Introduction  to  Channel  Modeling  for  Indoor  UWB  Channel 

In  the  S-V  model  [10],  the  arrival  of  clusters  is  modelled  as  a  Poisson  arrival  process  with  a  rate 
A,  while  within  each  cluster,  subsequent  multipath  contributions  or  rays  also  arrive  according  to 
a  Poisson  process  with  a  rate  A  (see  Fig.  7).  In  the  S-V  model,  the  magnitude  of  the  fc-th  path 
within  the  Z-th  cluster  follows  a  Rayleigh  distribution,  and  the  phase  of  each  path  is  assumed  to 
be  a  statistically  independent  random  variable  over  [0,27 r).  Besides,  the  average  Power  Decay 
Profile  (PDP)  is  characterized  by  an  exponential  decay  of  the  amplitude  of  the  clusters,  and  a 
different  exponential  decay  for  the  amplitude  of  the  received  pulses  within  each  cluster,  as  shown 
in  Fig.  8.  In  the  IEEE  UWB  indoor  channel  model  [6],  the  cluster  approach  was  adopted  (same 
as  S-V  model),  but  a  log-normal  distribution  was  suggested  for  characterizing  the  multi-path  gain 
amplitude,  and  an  additional  log-normal  variable  was  introbuced  for  representing  the  fluctuations 
of  the  total  multipath  gain.  Besides,  the  phase  of  each  path  is  assumed  to  be  either  0  or  n  with 
equal  probability. 

4.2  Outdoor  UWB  Channel  Modeling 
4.2.1  Cluster  Arrival  and  Power  Decay  Profile 

We  study  the  outdoor  UWB  signal  propagation  in  three  cases:  near  field  (less  than  55m),  medium 
field  (55m-85m),  and  far  field  (above  85m  and  up  to  120m  in  this  study).  In  the  data  collection, 
each  sample  is  spaced  at  50  picosecond  interval,  so  these  cases  are  corresponding  to  samples  1-7333 
for  near  field,  samples  7333-11333  for  medium  field,  and  samples  11334-16000  for  far  field.  In  Fig. 
6,  we  plot  the  power  profile  of  the  received  echos  (averaged  over  30  collections  to  eliminate  the 
effect  of  random  noise  and  each  collection  was  averaged  based  on  35  pulses)  for  the  three  different 
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cases.  Since  the  transmitted  pulse  (as  plotted  in  Fig.  2a)  is  a  very  narrow  impulse  pulse  (like 
a  delta  function  in  time  domain),  we  analyzed  the  channel  property  based  on  the  received  echos 
power  profile  plotted  in  Fig.  6,  and  similar  methodology  was  also  used  in  S-V  model  studies  [10]. 

Observe  Fig.  6,  multi-path  contributions  arrive  at  the  receiver  grouped  into  clusters.  The  time 
of  arrival  of  clusters  can  be  modeled  as  a  Poisson  arrival  process  with  a  rate  A,  while  within  each 
cluster,  subsequent  multipath  contributions  or  rays  also  arrive  according  to  a  Poisson  process  with 
a  rate  A  (see  Fig.  7).  We  define: 

•  Ti  =  the  arrival  time  of  the  first  path  of  the  Z-th  cluster; 

•  rk,l  =  the  delay  of  the  k- th  path  within  the  Z-th  cluster  relative  to  the  first  path  arrival  time 

Ti ; 

•  A  =  the  cluster  arrival  rate; 

•  A  =  the  ray  arrival  rate,  i.e.,  the  arrival  rate  of  the  paths  within  each  cluster. 

By  definition,  we  have  r0/  =  7}.  The  distributions  of  the  cluster  arrival  time  and  the  ray  arrival 
time  are  given  by 

p(7j|Ti_i)  =  Aexp(-A(7}  -  Ti_x),  Z  >  0 
P(Tk,l\T(k-l),l)  =  *exp  (— A(rte,i  -  r(fc_1)tJ)),  k  >  0  (2) 

The  above  observations  are  very  similar  as  that  for  the  indoor  UWB  channel.  Specifically,  we  also 
observed  the  A  and  A  are  quite  different  for  three  different  cases. 

•  Observe  Fig.  6a  for  near  field,  A  (1/ns)  is  around  0.02  (one  cluster  in  every  50ns  or  1000 
samples),  and  A  (1/ns)  is  around  0.4  (one  path  in  every  2.5ns  or  50  samples).  Perhaps  it’s 
because  some  major  scatters  in  near  field  (such  as  tree  stems)  reflected  signals,  so  some  paths 
are  quite  dominant. 
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•  Observe  Fig.  6b  for  medium  field,  clusters  arrive  quite  often.  A  (1/ns)  is  around  0.05  (one 
cluster  in  every  20ns  or  400  samples),  and  A  (1/ns)  is  around  1  (one  path  in  every  Ins  or  20 
samples). 

•  Observe  Fig.  6c  for  far  field,  clusters  almost  always  arrive  (because  of  rich  scattering),  so  A 
(1/ns)  is  around  0.5  (one  cluster  in  every  2 ns  or  20  samples),  and  A  (1/ns)  is  around  4  (one 
path  in  every  250ps  or  5  samples).  Perhaps  it’s  because  of  rich  scattering,  every  path  has 
very  similar  power  level. 

Besides,  the  average  PDP  can  be  represented  by  an  exponential  decay  of  the  amplitude  of  the 
clusters,  and  a  different  exponential  decay  for  the  amplitude  of  the  received  pulses  within  each 
cluster,  as  shown  in  Fig.  8. 

4.2.2  Statistical  Distribution  of  Channel  Coefficients 

We  also  study  the  statistical  distributions  of  each  given  path.  We  plot  the  histogram  for  some 
sample  values  of  the  above  three  cases  based  on  30  collections  and  each  collection  is  averaged  over 
35  pulses.  Near  field  samples  are  based  on  samples  5001-6000;  medium  field  samples  are  based  on 
samples  8001-9000;  and  far  field  samples  are  based  on  samples  12001-13000.  Since  the  samples  are 
very  close  (within  7.5m  distance),  so  their  path-loss  effect  can  be  ignored.  For  each  case,  we  have 
30000  samples,  and  we  plot  their  histogram  in  Fig.  9. 

First,  observe  Fig.  9c  for  fax  field,  the  histogram  can  be  almost  perfectly  modelled  by  a  non¬ 
zero-mean  Gaussian  distribution,  which  means  the  amplitude  of  the  channel  coefficient  follows  a 
Rician  distribution, 

Pa(x)  =  }/0(^)  X  >  0  (3) 

where  s  is  the  mean  value  of  Gaussian  and  /o(-)  is  the  zero  order  modified  Bessel  function.  This 
kind  of  channel  is  known  as  Rician  fading  channel.  A  Rician  channel  is  characterized  by  two 
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parameters,  Rician  factor  K  which  is  the  ratio  of  the  direct  path  power  to  that  of  the  multipath, 

i.e.,  K  =  s2/2ct2,  and  the  Doppler  spread  (or  single-sided  fading  bandwidth)  /<*.  Similarly,  Fig. 
9b  for  medium  field,  the  histogram  can  be  approximately  modelled  by  a  non-zero-mean  Gaussian 
distribution,  which  means  the  amplitude  of  the  channel  coefficient  follows  a  Rician  distribution. 
Observe  Fig.  9a  for  near  field,  the  histogram  can  be  approximately  modelled  by  two  non-zero-mean 
Gaussian  distributions,  which  means  it’s  noil-stationary,  and  the  amplitude  of  the  channel  coefficient 
follows  one  of  two  Rician  distributions.  The  above  observations  are  quite  different  with 
the  indoor  UWB  channel  model  (log-normal  distribution)  and  S-V  model  (Rayleigh 
distribution).  The  sign  of  channel  coefficient  is  either  -|-1  or  -1,  i.e.,  its  phase  is  either  0  or  7r, 
which  matches  the  IEEE  indoor  UWB  channel  model. 

5  Wireless  Channel  Modeling  Based  on  CLEAN  Method 

We  apply  the  CLEAN  algorithm  to  obtain  the  UWB  channel  model  based  on  the  transmitted  pulses 
and  received  echos.  The  CLEAN  algorithm  was  first  introduced  in  [5]  and  has  been  applied  to  UWB 
measurements  [2]  [11]  and  it  assumes  that  the  channel  is  a  series  of  impulses  which  is  consistent 
with  the  tapped-delay  line  channel  model.  This  algorithm  searches  the  received  echos  iteratively 
with  the  template  to  find  the  maximum  correlation  [1].  The  steps  are  [9]: 

1.  Calculate  the  autocorrelation  of  the  template  rss(t)  and  the  cross-correlation  of  the  template 
with  the  received  waveform  rsy{t)- 

2.  Find  the  largest  correlation  peak  in  rsy(£),  record  the  normalized  amplitudes  a/t  and  relative 
time  delay  r ^  of  the  correlation  peak. 

3.  Subtract  rss(t)  scaled  by  from  rsy(t)  at  the  time  delay  r^. 

4.  If  a  stopping  criterion  (e.g.,  a  minimum  threshold  on  the  peak  correlation)  is  not  met,  go  to 
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step  2.  Otherwise  stop. 


Based  on  the  CLEAN  method,  we  successfully  obtained  the  channel  impulse  responses  for  all 
transmit  waveforms  and  receive  echoes.  For  illustration  purposes,  in  Figs.  10,  11,  and  12,  we  plot 
the  channel  impulse  responses  for  200MHz,  400MHz,  and  UWB  channels  using  CLEAN  method  in 
two  experiments.  Observe  that  for  all  channels,  channel  impulse  responses  have  many  similarities: 
all  can  be  modeled  as  linear  time-variant  filter  channel,  which  is  a  more  general  case  of  the  S-V 
model. 


6  Conclusions 

In  this  paper,  we  studied  the  statistical  modeling  for  outdoor  wireless  channels  (200MHz,  400MHz, 
and  UWB)  in  rich  scattering  and  time- varying  environment  based  on  extensive  data  collected  using 
narrowband  and  UWB  radars.  We  validated  that  UWB  echo  signals  (within  a  burst)  don’t  hold 
self-similarity,  which  means  the  future  signals  can’t  be  forecasted  based  on  the  received  signals  and 
channel  modeling  is  necessary  from  statistical  point  of  view.  In  outdoor  UWB  channel,  the  multi- 
path  contributions  arrive  at  the  receiver  are  grouped  into  clusters.  The  time  of  arrival  of  clusters 
can  be  modeled  as  a  Poisson  arrival  process,  while  within  each  cluster,  subsequent  multipath 
contributions  or  rays  also  arrive  according  to  a  Poisson  process.  At  different  field  (near  field, 
medium  field,  and  far  field),  we  observed  that  the  Poisson  process  parameters  are  quite  different. 
We  also  observed  that  the  amplitude  of  channel  coefficient  at  each  path  follows  Rician  distribution 
for  medium  and  far  field,  and  it’s  non-stationary  for  paths  from  near  field  (one  of  two  Rician 
distributions),  and  these  observations  are  quite  different  with  the  IEEE  indoor  UWB  channel 
model  and  S-V  model.  Using  CLEAN  method,  we  observed  that  for  all  channels,  channel  impulse 
responses  have  many  similarities:  all  can  be  modeled  as  linear  time-variant  filter  channel,  which  is 
a  more  general  case  of  the  S-V  model. 
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Figure  1:  This  figure  shows  the  lift  with  the  experiment.  The  antennas  are  at  the  far  end  of  the 
lift  from  the  viewer  under  the  roof  that  was  built  to  shield  the  equipment  from  the  elements.  This 
picture  was  taken  in  September  with  the  foliage  largely  still  present.  The  cables  coming  from  the 
lift  are  a  ground  cable  to  an  earth  ground  and  one  of  4  tethers  used  in  windy  conditions. 
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Figure  2:  UWB  radar  Transmitted  pulse  and  received  echos  in  one  experiment,  (a)  Transmitted 
pulse,  (b)  Received  echos. 
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Figure  3:  Narrowband  radar  (200MHz)  transmitted  pulse  and  received  echos  in  one  experiment, 
(a)  Transmitted  pulse,  (b)  Received  echos. 
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Figure  4:  Narrowband  radar  (400MHz)  transmitted  pulse  and  received  echos  in  one  experiment, 
(a)  Transmitted  pulse,  (b)  Received  echos. 
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Figure  5:  The  variance-time  plot  of  10  UWB  data  collections,  which  demonstrates  that  UWB 
reflected  signals  are  not  self-similar  within  each  collection. 
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Figure  6:  The  power  profile  for  three  different  cases:  (a)  near  field,  (b)  medium  field,  and  (c)  far 
field. 
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Path  Magnitude 


Figure  7:  An  illustration  of  the  channel  impulse  response  in  S-V  model. 


Power  Delay 
Profile 


Figure  8:  An  illustration  of  the  double  exponential  decay  of  the  mean  cluster  power  and  the  ray 
power  within  clusters  in  S-V  model. 
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(c) 


Figure  9:  The  histograms  and  their  approximation  using  Gaussian  distributions  (dashed  lines). 
The  histograms  are  based  on  30  collections  and  each  collection  is  averaged  over  35  pulses,  (a)  near 
field  samples,  (b)  medium  field  samples,  and  (c)  far  field  samples. 
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Figure  10:  The  channel  impulse  responses  for  200MHz  channel  using  CLEAN  method  in  two 
experiments. 


22 


Page  79  of  134 


,  x  10' 


(b) 


Figure  11:  The  channel  impulse  responses  for  400MHz  channels  using  CLEAN  method  in  two 
experiments. 
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Figure  12:  The  channel  impulse  responses  for  UWB  channels  using  CLEAN  method  in  two  exper¬ 
iments. 
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Abstract — Inspired  by  human's  innate  ability  to  process 
and  integrate  information  from  disparate,  network-based 
sources,  we  apply  human-inspired  information  integration 
mechanisms  to  target  detection  in  UWB  radar  sensor  network. 
Humans'  information  integration  mechanisms  have  been  mod¬ 
elled  using  maximum-likelihood  estimation  (MLE)  or  soft- 
max  approaches.  In  this  paper,  we  apply  these  two  algorithms 
to  cognitive  radar  sensor  networks  target  detection.  Discrete- 
cosine-transform  (DCT)  is  used  to  process  the  integrated  data 
from  MLE  or  soft-max.  We  apply  fuzzy  logic  system  (FLS) 
to  automatic  target  detection  based  on  the  AC  power  values 
from  DCT.  Simulation  results  show  that  our  MLE-DCT-FLS 
and  soft-max-DCT-FLS  approaches  perform  very  well  in  the 
radar  sensor  network  target  detection,  whereas  the  existing 
2-D  construction  algorithm  doesn't  work  in  this  study. 

Index  Terms  :  Human-inspired,  UWB  radar  sensor  networks, 
fuzzy  logic  systems,  automatic  target  recognition,  maximum- 
likelihood  estimation. 


I.  Introduction  and  Motivation 

In  real  world,  UWB  radar  sensor  network  information 
integration  is  necessary  in  different  applications.  For  ex¬ 
ample,  in  an  emergency  natural  disaster  scenario,  such  as 
Utah  Mine  Collapse  in  August  2007  or  West  Virginia  Sago 
mine  disaster  in  January  2006,  UWB  radar  sensor  network- 
based  information  integration  for  first  responders  is  critical 
for  search  and  rescue.  Danger  may  appear  anywhere  at  any 
time,  therefore,  first  responders  must  monitor  a  large  area 
continuously  in  order  to  identify  potential  danger  and  take 
actions.  Due  to  the  dynamic  and  complex  nature  of  natural 
disaster,  some  buried/foleage  victims  may  not  be  found 
with  image/video  sensors,  and  UWB  radar  sensors  are 
needed  for  penetrating  the  ground  or  sense-through-wall. 
Unfortunately,  the  radar  data  acquired  are  often  limited  and 
noisy.  Unlike  medical  imaging  or  synthetic  aperture  radar 
imaging  where  abundance  of  data  is  generally  available 
through  multiple  looks  and  where  processing  time  may 
not  be  crucial,  practical  cognitive  radar  sensor  networks 
are  typically  the  opposite:  availability  of  data  is  limited 
and  required  processing  time  is  short.  This  need  is  also 
motivated  by  the  fact  that  humans  display  a  remarkable 
capability  to  quickly  perform  target  recognition  despite 
noisy  sensory  signals  and  conflicting  inputs.  Humans  are 


adept  at  network  visualization,  and  at  understanding  subtle 
implications  among  the  network  connections.  To  date, 
however,  human’s  innate  ability  to  process  and  integrate 
information  from  disparate,  network-based  sources  for  sit¬ 
uational  understanding  has  not  translated  well  to  automated 
systems.  In  this  paper,  we  apply  human  information  inte¬ 
gration  mechanisms  to  information  fusion  in  UWB  radar 
sensor  network. 

II.  Sense-through-Foliage  Radar  Sensor 
Networks  Data  Measurement  and  Collection 

Our  work  is  based  on  the  sense-through-foliage  UWB 
radar  sensor  networks.  The  foliage  experiment  was  con¬ 
structed  on  a  seven-ton  man  lift,  which  had  a  total  lifting 
capacity  of  450  kg.  The  limit  of  the  lifting  capacity  was 
reached  during  the  experiment  as  essentially  the  entire 
measuring  apparatus  was  placed  on  the  lift.  The  principle 
pieces  of  equipment  secured  on  the  lift  are:  Barth  pulser, 
Tektronix  model  7704  B  oscilloscope,  dual  antenna  mount¬ 
ing  stand,  two  antennas,  rack  system,  IBM  laptop,  HP 
signal  Generator,  Custom  RF  switch  and  power  supply  and 
Weather  shield  (small  hut).  The  target  is  a  trihedral  reflector 
(as  shown  in  Fig.  1).  Throughout  this  work,  a  Barth  pulse 
source  (Barth  Electronics,  Inc.  model  732  GL)  was  used. 
The  pulse  generator  uses  a  coaxial  reed  switch  to  discharge 
a  charge  line  for  a  very  fast  rise  time  pulse  outputs.  The 
model  732  pulse  generator  provides  pulses  of  less  than  50 
picoseconds  (ps)  rise  time,  with  amplitude  from  1 50  V  to 
greater  than  2  KV  into  any  load  impedance  through  a  50 
ohm  coaxial  line.  The  generator  is  capable  of  producing 
pulses  with  a  minimum  width  of  750  ps  and  a  maximum 
of  1  microsecond.  This  output  pulse  width  is  determined  by 
charge  line  length  for  rectangular  pulses,  or  by  capacitors 
for  1/e  decay  pulses. 

For  the  data  we  used  in  this  paper,  each  sample  is 
spaced  at  50  picosecond  interval,  and  16,000  samples  were 
collected  for  each  collection  for  a  total  time  duration  of 
0.8  microseconds  at  a  rate  of  approximately  20  Hz.  We 
plot  the  transmitted  pulse  (one  realization)  in  Fig.  2a)  and 
the  received  echos  in  one  collection  in  Fig.  2b  (averaged 
over  35  pulses). 
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Fig.  1.  The  target  (a  trihedral  reflector)  is  shown  on  the  stand  at  300 
feet  from  the  lift. 


III.  Human  Information  Integration 
Mechanisms 

Recently,  a  maximum-likelihood  estimation  (MLE)  ap¬ 
proach  was  proposed  for  multi-sensory  data  fusion  in  hu¬ 
man  [4].  In  the  MLE  approach  [4],  sensory  estimates  of  an 
environmental  property  can  be  represented  by  Sj  =  fi(S) 
where  S  is  the  physical  property  being  estimated,  /  is  the 
operation  the  nervous  system  performs  to  derive  the  esti¬ 
mate,  and  S  is  the  perceptual  estimate.  Sensory  estimates 
are  subject  to  two  types  of  error:  random  measurement 
error  and  bias.  Thus,  estimates  of  the  same  object  property 
from  different  cues  usually  differ.  To  reconcile  the  discrep¬ 
ancy,  the  nervous  system  must  either  combine  estimates 
or  choose  one,  thereby  ignoring  the  other  cues.  Assuming 
that  each  single-cue  estimate  is  unbiased  but  corrupted 
by  independent  Gaussian  noise,  the  statistically  optimal 
strategy  for  cue  combination  is  a  weighted  average  [4] 

M 

sc=y:  wjSj  (i) 

t=i 

1  /a2 

where  Wi  =  ^  and  is  the  weight  given  to  the  zth 
single-cue  estimate,  erf  is  that  estimates  variance,  and  M 
is  the  total  number  of  cues.  Combining  estimates  by  this 
MLE  rule  yields  the  least  variable  estimate  of  S  and  thus 
more  precise  estimates  of  object  properties. 

Besides,  some  other  summation  rules  have  been  pro¬ 
posed  interception  and  cognition  such  as  soft-max  rule: 
V  =  (Ei=i  x? )n  [3]  where  Xi  denotes  the  input  from  an 
input  source  i ,  and  M  is  the  total  number  of  sources.  In  this 
paper,  we  will  apply  MLE  and  soft-max  human  brain  in¬ 
formation  integration  mechanisms  to  cognitive  radar  sensor 
network  information  integration. 

IV.  Human-Inspired  Sense-through-Foliage 
Target  Detection 

In  Figs.  3a  and  3b,  we  plot  two  collections  of  UWB 
radars.  Fig.  3a  has  no  target  on  range,  and  Fig.  3b  has 
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Fig.  2.  Transmitted  pulse  and  received  echos  in  one  experiment,  (a) 
Transmitted  pulse,  (b)  Received  echos. 


target  at  samples  around  13,900.  We  plot  the  echo  differ¬ 
ences  between  Figs.  3a  and  3b  in  Fig.  3c.  However,  it  is 
impossible  to  identify  whether  there  is  any  target  and  where 
there  is  target  based  on  Fig.  3c.  Since  significant  pulse-to- 
pulse  variability  exists  in  the  echos,  this  motivate  us  to 
explore  the  spatial  and  time  diversity  using  Radar  Sensor 
Networks  (RSN). 

In  RSN,  the  radar  sensors  are  networked  together  in  an 
ad  hoc  fashion.  They  do  not  rely  on  a  preexisting  fixed 
infrastructure,  such  as  a  wireline  backbone  network  or  a 
base  station.  They  are  self-organizing  entities  that  are  de¬ 
ployed  on  demand  in  support  of  various  events  surveillance, 
battlefield,  disaster  relief,  search  and  rescue,  etc.  Scalability 
concern  suggests  a  hierarchical  organization  of  radar  sensor 
networks  with  the  lowest  level  in  the  hierarchy  being  a 
cluster.  As  argued  in  [7]  [6]  [5]  [9],  in  addition  to  helping 
with  scalability  and  robustness,  aggregating  sensor  nodes 
into  clusters  has  additional  benefits: 
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Fig.  3.  Measurement  with  poor  signal  quality  and  35  pulses  average,  (a) 
Expanded  view  of  traces  (no  target)  from  sample  13,001  to  15,000.  (b) 
Expanded  view  of  traces  (with  target)  from  samples  13,001  to  15,000.  (c) 
The  differences  between  (a)  and  (b). 

1 )  conserving  radio  resources  such  as  bandwidth; 

2)  promoting  spatial  code  reuse  and  frequency  reuse; 

3)  simplifying  the  topology,  e.g.,  when  a  mobile  radar 
changes  its  location,  it  is  sufficient  for  only  the 
nodes  in  attended  clusters  to  update  their  topology 
information; 

4)  reducing  the  generation  and  propagation  of  routing 
information;  and, 

5)  concealing  the  details  of  global  network  topology 
from  individual  nodes. 

In  RSN,  each  radar  can  provide  their  pulse  parameters  such 
as  timing  to  their  clusterhead  radar,  and  the  clusterhead 
radar  can  combine  the  echos  (RF  returns)  from  the  target 
and  clutter.  In  this  paper,  we  propose  a  RAKE  structure  for 


combining  echos,  as  illustrated  by  Fig.  4.  The  integration 
means  time-average  for  a  sample  duration  T  and  it’s  for 
general  case  when  the  echos  are  not  in  discrete  values.  It  is 
quite  often  assumed  that  the  radar  sensor  platform  will  have 
access  to  Global  Positioning  Service  (GPS)  and  Inertial 
Navigation  Unit  (INU)  timing  and  navigation  data  [1].  In 
this  paper,  we  assume  the  radar  sensors  are  synchronized 
in  RSN.  In  Fig.  4,  the  echo,  i.e.,  RF  response  by  the  pulse 
of  each  cluster-member  sensor,  will  be  combined  by  the 
clusterhead  using  a  weighted  average,  and  the  weight  Wi 
is  determined  by  the  two  human-inspired  mechanisms. 


Fig.  4.  Echo  combining  by  clusterhead  in  RSN. 

We  applied  the  human-inspired  MLE  algorithm  to  com¬ 
bine  the  sensed  echo  collection  from  M  =  30  UWB  radars, 
and  then  the  combined  data  are  processed  using  discrete- 
cosine  transform  (DCT)  to  obtain  the  AC  values.  Based 
on  our  experiences,  echo  with  a  target  generally  has  high 
and  nonfluctuating  AC  values  and  the  AC  values  can  be 
obtained  using  DCT.  We  plot  the  power  of  AC  values  in 
Figs.  5a  and  5b  using  MLE  and  DCT  algorithms  for  the  two 
cases  (with  target  and  without  target)  respectively.  Observe 
that  in  Fig.  5b,  the  power  of  AC  values  (around  sam¬ 
ple  13,900)  where  the  target  is  located  is  non-fluctuating 
(somehow  monotonically  increase  then  decrease).  Although 
some  other  samples  also  have  very  high  AC  power  values,  it 
is  very  clear  that  they  are  quite  fluctuating  and  the  power  of 
AC  values  behaves  like  random  noise  because  generally  the 
clutter  has  Gaussian  distribution  in  the  frequency  domain. 

Similarly,  we  applied  the  soft-max  algorithm  (n  =  2)  to 
combine  the  sensed  echo  collection  from  M  =  30  UWB 
radars,  and  then  used  DCT  to  obtain  the  AC  values.  We 
plot  the  power  of  AC  values  in  Figs.  5a  and  5b  using  soft- 
max  and  DCT  algorithms  for  the  two  cases  (with  target 
and  without  target)  respectively.  Observe  that  in  Fig.  6b, 
the  power  of  AC  values  (around  sample  13,900)  where  the 
target  is  located  is  non-fluctuating  (somehow  monotonically 
increase  then  decrease). 

We  made  the  above  observations.  However,  in  real 
world  application,  automatic  target  detection  is  necessary  to 
ensure  that  our  algorithms  could  be  performed  in  real  time. 
In  Section  V,  we  apply  fuzzy  logic  systems  to  automatic 
target  detection  based  on  the  power  of  AC  values  (obtained 
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via  MLE-DCT  or  soft-max-DCT). 


Sample  Index  x 


(b) 


Fig.  6.  Power  of  AC  values  using  soft-max  based  information  integration 
and  DCT.  (a)  No  target  (b)  With  target  in  the  field. 


Fig.  5.  Power  of  AC  values  using  MLE-based  information  integration 
and  DCT.  (a)  No  target  (b)  With  target  in  the  field. 

We  compared  our  approaches  to  the  scheme  proposed  in 
[10].  In  [10],  2-D  image  was  created  via  adding  voltages 
with  the  appropriate  time  offset.  In  Figs.  7a  and  7b,  we  plot 
the  2-D  image  created  based  on  the  above  two  data  sets 
(from  samples  13,800  to  14,200).  The  sensed  data  from  30 
radars  are  averaged  first,  then  plotted  in  2-D  [10].  However, 
it’s  not  clear  which  image  shows  there  is  target  on  range. 

V.  Fuzzy  Logic  System  for  Automatic  Target 
Detection 

A.  Overview  of  Fuzzy  Logic  Systems 

Figure  8  shows  the  structure  of  a  fuzzy  logic  system 
(FLS)  [8].  When  an  input  is  applied  to  a  FLS,  the  inference 
engine  computes  the  output  set  corresponding  to  each  rule. 
The  defuzzifer  then  computes  a  crisp  output  from  these 
rule  output  sets.  Consider  a  p-input  1 -output  FLS,  using 
singleton  fuzzification,  center-of-sets  defuzzification  [8] 
and  “IF-THEN”  rules  of  the  form 


Rl  :  IF  x\  is  F[  and  X2  is  F2  and  •  •  •  and  xp  is  Fj,, 
THEN  y  is  G*. 


Assuming  singleton  fuzzification,  when  an  input  x'  = 
{xi, . . . ,  x^}  is  applied,  the  degree  of  firing  corresponding 
to  the  /th  rule  is  computed  as 


Mf‘  04)  *  •  •  •  */*F|>p)  =  7;= iMf;(*'<)  (2) 

where  ★  and  T  both  indicate  the  chosen  f-norm.  There 
are  many  kinds  of  defuzzifiers.  In  this  paper,  we  focus,  for 
illustrative  purposes,  on  the  center-of-sets  defuzzifier  [8].  It 
computes  a  crisp  output  for  the  FLS  by  first  computing  the 
centroid,  cGi,  of  every  consequent  set  G/,  and,  then  com¬ 
puting  a  weighted  average  of  these  centroids.  The  weight 
corresponding  to  the  /th  rule  consequent  centroid  is  the 
degree  of  firing  associated  with  the  /th  rule,  T'i^1pF/ (x'J, 
so  that 


,  „  £,=iCg<7?=1Mf;(z;) 

2/“»(x  )  =  y>M  - 77T 


(3) 


where  M  is  the  number  of  rules  in  the  FLS.  In  this  paper, 
we  design  a  FLS  for  automatic  target  recognition  based  on 
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Fig.  7.  2-D  image  created  via  adding  voltages  with  the  appropriate  time 
offset,  (a)  No  target  (b)  With  target  in  the  field. 


the  AC  values  obtained  using  MLE-DCT  or  soft-max-DCT. 


FUZZY  LOGIC  SYSTEM 


Fig.  8.  The  structure  of  a  fuzzy  logic  system. 


B.  FLS  for  Automatic  Target  Detection 

Observe  that  in  Figs.  5  and  6,  the  power  of  AC  values 
are  quite  fluctuating  and  have  lots  of  uncertainties.  FLS  is 
well  known  to  handle  the  uncertainties.  For  convenience  in 
describing  the  FLS  design  for  Automatic  Target  Detection 
(ATD),  we  first  give  the  definition  of  footprint  of  uncer¬ 
tainty  of  AC  power  values  and  region  of  interest  in  the 
footprint  of  uncertainty. 


Definition  1  (Footprint  of  Uncertainty):  Uncertainty  in 
the  AC  power  values  and  time  index  consists  of  a  bounded 
region,  that  we  call  the  footprint  of  uncertainty  of  AC 
power  values.  It  is  the  union  of  all  AC  power  values. 

Definition  2  (Region  of  Interest  (Rol)):  An  Rol  in  the 
footprint  of  uncertainty  is  a  contour  consisting  of  a  large 
number  (greater  than  50)  of  AC  power  values  where  AC 
power  values  increase  then  decrease. 

Definition  3  (Fluctuating  Point  in  Rol):  P(i)  is  called  a 
fluctuating  point  in  the  Rol  if  P(i  -  1),  P(i),  P(i  4-1)  are 
non-monotonically  increasing  or  decreasing. 

Our  FLS  for  automatic  target  detection  will  classify  each 
ROI  (with  target  or  no  target)  based  on  two  antecedents: 
the  centroid  of  the  ROI  and  the  number  of  fluctuating 
points  in  the  ROI.  The  linguistic  variables  used  to  represent 
these  two  antecedents  were  divided  into  three  levels:  low, 
moderate ,  and  high.  The  consequent  -  the  possibility  that 
there  is  a  target  at  this  Rol  -  was  divided  into  5  levels, 
Very  Strong,  Strong,  Medium,  Weak,  Very  Weak.  We  used 
trapezoidal  membership  functions  (MFs)  to  represent  low, 
high ,  very  strong ,  and  very  weak ;  and  triangle  MFs  to 
represent  moderate ,  strong ,  medium ,  and  weak.  All  inputs 
to  the  antecedents  are  normalized  to  0-10. 

Based  on  the  fact  the  AC  power  value  of  target  is 
non-fluctuating  (somehow  monotonically  increase  then  de¬ 
crease),  and  the  AC  power  value  of  clutter  behaves  like 
random  noise  because  generally  the  clutter  has  Gaussian 
distribution  in  the  frequency  domain,  we  design  a  fuzzy 
logic  system  using  rules  such  as: 

Rl  :  IF  centroid  of  a  Rol  (xi)  is  Fj,  and  the  number  of 
fluctuating  points  in  the  ROI  ( xv )  is  F^,  THEN  the 
possibility  that  there  is  a  target  at  this  Rol  (y)  is  Gz. 

where  /  =  1, . . . ,  9.  We  summarize  all  the  rules  in  Table  I. 
For  every  input  (xi,X2),  the  output  is  computed  using 


y{xux2) 


£LlMF‘(si)MFjKs2)Cqv9 

£?=1MF‘(Z1)MF?(Z2) 


(4) 


We  ran  simulations  to  1000  collections  in  the  real  world 
sense-through-foliage  experiment,  and  found  that  our  FLS 
performs  very  well  in  the  automatic  target  detection  based 
on  the  AC  power  values  obtained  from  MLE-DCT  or  soft- 
max-DCT,  and  achieve  probability  of  detection  pd  =  100% 
and  false  alarm  rate  pja  =  0. 


VI.  Conclusions 

Inspired  by  human’s  innate  ability  to  process  and  inte¬ 
grate  information  from  disparate,  network-based  sources, 
we  applied  human-inspired  information  integration  mecha¬ 
nisms  to  target  detection  in  UWB  radar  sensor  network. 
Humans’  information  integration  mechanisms  have  been 
modelled  using  maximum-likelihood  estimation  (MLE)  or 
soft-max  approaches.  In  this  paper,  we  applied  these  two 
algorithms  to  UWB  radar  sensor  networks  target  detection. 
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TABLE  I 

The  rules  for  target  detection.  Antecedent  1  is  centroid  of  a 
Rol,  ANTECEDENT  2  IS  the  number  of  fluctuating  points  in  the  ROI, 
AND  CONSEQUENT  IS  the  possibility  that  there  is  a  target  at  this  Rol. 


Rule  # 

Antecedent  1 

Antecedent  2 

Consequent 

1 

low 

low 

medium 

2 

low 

moderate 

weak 

3 

low 

high 

very  weak 

4 

moderate 

low 

strong 

5 

moderate 

moderate 

medium 

6 

moderate 

high 

weak 

7 

high 

low 

very  strong 

8 

high 

moderate 

strong 

9 

high 

high 

medium 

Discrete-cosine-transform  (DCT)  was  used  to  process  the 
integrated  data  from  MLE  or  soft-max.  We  applied  fuzzy 
logic  system  (FLS)  to  automatic  target  detection  based 
on  the  AC  power  values  from  DCT.  Simulation  results 
showed  that  our  MLE-DCT-FLS  and  soft-max-DCT-FLS 
approaches  performed  very  well  in  the  radar  sensor  network 
target  detection,  whereas  the  existing  2-D  construction 
algorithm  couldn’t  work  in  this  study. 
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Abstract 

In  this  paper,  we  study  the  statistical  modeling  for  outdoor  non  line-of-sight  (NLOS)  channel 
in  rich  scattering  and  time-varying  foliage  environment  based  on  extensive  data  collected  using 
both  narrowband  and  ultra-wideband  (UWB)  radars.  The  multi-path  contributions  arrive  at 
the  receiver  are  grouped  into  clusters.  The  time  of  arrival  of  clusters  can  be  modeled  as  a 
Poisson  arrival  process,  while  within  each  cluster,  subsequent  multipath  contributions  or  rays 
also  arrive  according  to  a  Poisson  process.  However,  the  parameters  are  quite  different  along 
with  the  frequency.  We  also  observe  that  the  amplitude  of  channel  coefficient  at  each  path 
can  be  more  accurately  characterized  by  log-logistic  distribution  (LLD)  other  than  log-normal, 
Weibull  or  Rayleigh  due  to  the  best  goodness-of-fit  and  smallest  root-mean-square  (RMS). 

Index  Terms  : Propagation  channel,  outdoor,  NLOS,  log-logistic,  goodness-of-fit,  narrow- 
band,  UWB 

1  Introduction  and  Motivation 

The  area  of  indoor  radio  propagation  channel  measurement  and  modeling  has  been  intensively 
investigated  [1]  [2]  [3]  because  if  the  characterization  of  propagation  channel  is  known,  efforts  can  be 
employed  at  both  transmitters  and  receivers  to  overcome  the  degrading  effects  and  therefore  improve 
the  communication  performance.  As  far  as  military  warefare  is  concerned,  lots  of  applications  are  for 
outdoor  activities  such  as  sense-through- foliage  target  detection.  It  is  believed  that  understanding 
the  characterization  of  outdoor  non  line-of-sight  (NLOS)  channels  will  assist  both  target  detection 
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as  well  as  commercial  wireless  communications. 


Like  indoor  environment,  the  foliage  contains  a  wealth  of  multiple  scattering  other  than  LOS 
free  space.  In  addition,  the  movement  of  leaves,  branches  and  even  the  tree  trunks  contribute  to 
the  time-variance  fading  phenomenon.  Since  the  foliage  medium  can  be  completely  described  by 
its  time  and  space  varying  feature,  one  can  investigate  the  channel  model  based  on  characterization 
of  the  channel  impulse  response  (CIR). 

In  our  investigation,  we  will  apply  both  narrowband  and  ultra-wide  band  (UWB)  radar  to  model 
the  propagation  channel,  as  we  believe  that  foliage  is  composed  of  intervening  materials  that  are 
electromagnet ically  dispersive,  which  contributes  to  the  strong  frequency  dependence  of  foliage, 
and  thus  a  narrowband-wideband  study  would  assist  with  the  better  understanding  of  statistic 
property  of  the  channel.  Narrowband  signals  have  been  tried  at  200,  400  and  600  megahertz 
respectively,  while  UWB  radar  emissions  are  at  a  relatively  low  frequency- typically  between  100 
MHz  and  3  GHz.  Each  frequency  component  in  a  radar  signal  will  sense  the  foliage  clutter  in  a 
slightly  different  manner,  therefore  provides  difference  in  multipath. 

The  rest  of  this  paper  is  organized  as  follows.  In  Section  2,  we  summarize  the  measurement 
and  collection  of  the  data.  In  Section  3,  we  apply  CLEAN  algorithm  to  extract  CIR  for  200MHz, 
400MHZ,  600MHz  and  UWB  signals.  Section  4  presents  the  channel  model  in  view  of  temporal 
characterization  as  well  as  statistic  model  comparison.  We  conclude  our  work  in  Section  5. 

2  Measurement  Setup 

The  foliage  penetration  measurement  effort  began  in  August  2005  and  continued  through  December 
2005.  Working  in  August  through  the  fall  of  2005,  the  foliage  measured  included  late  summer 
foliage  and  fall  and  early  winter  foliage.  Late  summer  foliage,  because  of  the  limited  rainfall, 
involved  foliage  with  decreased  water  content.  Late  fall  and  winter  measurements  involved  largely 
defoliated  but  dense  forest,  so  it’s  a  rich  scattering  environment.  Because  of  wind  or  different 
temperature  in  dense  forest,  it’s  also  a  time-varying  environment. 

The  radar  experiment  was  constructed  on  a  seven- ton  man  lift,  which  had  a  total  lifting  capacity 
of  450  kg.  The  limit  of  the  lifting  capacity  was  reached  during  the  experiment  as  essentially  the 
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entire  measuring  apparatus  was  placed  on  the  lift  (as  shown  in  Fig.  1).  The  principle  pieces  of 
equipment  secured  on  the  lift  are:  dual  antenna  mounting  stand;  two  antennas;  200MHz,  lKw 
Amplifier,  power  supply,  pre-amp;  400MHz,  lKw  Amplifier,  power  supply,  pre-amp;  600MHz, 
lKw  Amplifier,  power  supply,  pre-amp;  Tektronix  model  7704  B  oscilloscope;  rack  system;  IBM 
laptop;  HP  signal  Generator;  Custom  RF  switch  and  power  supply  and  weather  shield  (small  hut). 
Particularly  for  UWB  signal,  a  Barth  pulse  source  (Barth  Electronics,  Inc.  model  732  GL)  was 
used.  The  pulse  generator  uses  a  coaxial  reed  switch  to  discharge  a  charge  line  for  a  very  fast  rise 
time  pulse  outputs.  The  model  732  pulse  generator  provides  pulses  of  less  than  50  picoseconds  (ps) 
rise  time,  with  amplitude  from  150  V  to  greater  than  2  KV  into  any  load  impedance  through  a  50 
ohm  coaxial  line.  The  generator  is  capable  of  producing  pulses  with  a  minimum  width  of  750  ps 
and  a  maximum  of  1  microsecond.  This  output  pulse  width  is  determined  by  charge  line  length  for 
rectangular  pulses,  or  by  capacitors  for  1/e  decay  pulses. 

The  system  was  pointing  at  the  specified  250  feet  one  way  distance.  For  the  data  we  used  in  this 
paper,  each  sample  is  spaced  at  50  picoseconds  interval,  and  16,000  samples  were  collected  for  each 
collection  for  a  total  time  duration  of  0.8  microseconds.  The  accomplished  data  structure  is  shown 
in  Fig.  2.  For  each  frequency  band,  12  positions  named  ’0,,,2’,’4,,’5,,,6’,,8,,T0,,T2’,T4,,T6,,T8’ 
and  ’20’  have  been  tested.  35  pulses  have  been  obtained  for  the  transmitted  and  received  signal 
respectively  at  each  position.  ’EMCO’  is  the  brand  name  of  the  antenna  horn.  ’Background’  refers 
to  the  background  data  taken  along  the  track  with  the  antenna  bore  sighted  to  the  targets. 

In  the  foliage-target  data,  we  fired  from  each  of  twelve  positions  along  the  track.  The  firing 
was  into  the  foliage  at  right  angles  to  the  track  and  these  right  angle  shots  are  called  track  data, 
i.e.,  ’Track’.  At  all  twelve  positions  along  the  track,  the  azimuth  angle  of  the  transmit  antenna  was 
changed  from  plus  10  to  -  10  degrees  to  provide  radial  measurements,  i.e.,  ’Radial’. 
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3  Channel  Impulse  Response  Based  on  the  Measured  Data  and 
CLEAN  Algorithm 

The  average  transmitted  and  received  pulses  at  position  4  for  different  frequencies  have  been  il¬ 
lustrated  from  Figs.  3-6.  The  purpose  of  average  is  to  remove  white  Gaussian  noise.  Note  that 
at  a  different  position  the  result  will  be  slightly  different.  However,  illustration  at  one  position  is 
sufficient  enough  to  describe  the  characterization. 

The  complicated  multipath  and  time- varying  CIR  can  be  modeled  as  follows  [4] 

r(t)  »  X]  "  T")  (!) 

n 

where  an  and  rn  is  referred  to  as  the  amplitude  and  delay  of  the  nth  propagation  path.  In  order 
to  extract  the  CIR  from  our  measurement,  the  CLEAN  algorithm  have  been  used.  It  was  initially 
introduced  in  [5]  to  enhance  radio  astronomical  maps  of  the  sky,  and  has  also  been  employed  in 
narrowband  channel  modeling  [6]  [7];  also  applied  in  UWB  channel  characterization  problems  [4]  [8]. 
Our  steps  involved  [8]  are: 

1.  Calculate  the  autocorrelation  of  the  transmitted  signal  Rss(t )  and  the  cross-correlation  of  the 
transmitted  with  the  received  waveform  Rsy(t). 

2.  Find  the  largest  correlation  peak  in  Rsy(t ),  record  the  normalized  amplitudes  a k  and  the 
relative  time  delay  r k  of  the  correlation  peak. 

3.  Subtract  Rss(t)  scaled  by  a k  from  Rsy(t)  at  the  time  delay  r^. 

4.  If  a  stopping  criterion  (a  minimum  threshold  on  the  peak  correlation  is  not  met,  go  step  2. 
otherwise  stop.) 

Given  the  transmission,  reception  and  the  CLEAN  processing  described  above,  the  obtained  CIR 
are  illustrated  from  Fig.  7  to  Fig.  10.  Note  that  we  plot  the  absolute  value  of  the  UWB  channel 
for  the  comparison  between  the  outdoor  UWB  channel  with  the  indoor  S-V  model  [9]  (see  Fig.  10). 
It  is  shown  that 
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1.  Both  narrowband  and  UWB  channels  are  made  up  of  multipath  components.  The  time- 
varying  path  magnitude  implies  that  fading  generally  exists,  therefore  the  received  pulses  are 
fairly  random  from  one  time  to  another. 

2.  The  UWB  channel  we  obtained  looks  similar  as  the  CIR  in  S-V  model  shown  in  Fig.  10. 
However,  they  are  different  in  the  arrival  of  cluster  and  subsequent  rays  as  well  as  envelope 
decay. 

3.  In  outdoor  environment,  the  largest  scattering,  i.e.,  the  highest  magnitude  does  not  always 
appear  at  the  first  path.  This  phenomenon  is  clearly  illustrated  in  Fig.  7  and  10. 

4.  Channels  are  frequency  dependent.  It  has  been  observed  that  the  intervening  materials,  such 
as  foliage  and  soil,  have  dielectric  properties  that  are  strongly  frequency  dependent.  This  in 
part  explains  the  difference  among  those  channels.  We  will  further  analyze  the  detail  in  the 
following  section. 

4  Outdoor  Channel  Modeling 

4.1  Temporal  Characterization 

In  the  S-V  model,  the  arrival  of  clusters  is  modelled  as  a  Poisson  arrival  process  with  a  rate  A, 
while  within  each  cluster,  subsequent  multipath  contributions  or  rays  also  arrive  according  to  a 
Poisson  process  with  a  rate  A  (see  Fig.  12).  Observe  Fig.  7~10,  like  in  S-V  model,  multipath 
contributions  arrive  at  the  receiver  grouped  into  clusters  and  therefore  similar  methodology  used 
in  S-V  model  studies  may  be  also  applied  to  200MHz,  400MHz,  600MHz  and  UWB  CIR.  The  time 
of  arrival  of  clusters  can  be  modeled  as  a  Poisson  arrival  process  with  a  rate  A,  while  within  each 
cluster,  subsequent  multipath  contributions  or  rays  also  arrive  according  to  a  Poisson  process  with 
a  rate  A.  We  define: 

•  T[  =  the  arrival  time  of  the  first  path  of  the  l-th  cluster; 

•  rk,l  =  the  delay  of  the  k-th  path  within  the  l-th  cluster  relative  to  the  first  path  arrival  time 

Tv, 
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•  A  =  the  cluster  arrival  rate; 


•  A  =  the  ray  arrival  rate,  i.e.,  the  arrival  rate  of  the  paths  within  each  cluster. 

By  definition,  we  have  tq/  =  7).  The  distributions  of  the  cluster  arrival  time  and  the  ray  arrival 
time  are  given  by 


p(2/|Tj_i)  =  Aexp  (— A(7J  -  7]_ i),Z  >  0 
p(Tk,i\T{k-i),l)  =  *exp  (-A(rfc,i  -  T(fc_1)t,)),  k>  0  (2) 

Specifically,  we  also  observed  that  the  A  and  A  are  quite  different  for  200MHz,  400MHz,  600MHz 
and  UWB  CIR.  We  listed  observed  parameters  in  Table  1.  As  for  indoor  UWB  data,  we  refer  [11]. 
The  higher  A  and  A  of  UWB  implies  its  exceptional  range  resolution.  Lower  A  and  A  of  outdoor 
UWB  than  those  of  indoor  means  indoor  environment  typically  have  a  richer  multiple  scattering. 

4.2  Statistical  Distribution  of  Channel  Amplitude 

In  the  S-V  model,  the  average  Power  Decay  Profile  (PDP)  is  characterized  by  an  exponential  decay 
of  the  clusters  and  a  different  exponential  decay  for  the  pulses  within  each  cluster.  In  other  words, 
the  amplitude  follows  rayleigh  distribution.  In  the  IEEE  UWB  indoor  channel  model  [10],  the 
clutter  approach  was  adopted  (same  as  S-V  model),  but  a  lognormal  distribution  was  suggested  for 
characterizing  the  multipath  gain  amplitude,  and  an  additional  log-normal  variable  was  introduced 
for  representing  the  fluctuations  of  the  total  multipath  gain.  In  this  Section  we  propose  that  log- 
logistic  model  that  would  better  characterize  the  multipath  gain  amplitude  for  both  outdoor  NLOS 
narrowband  and  UWB. 

4.2.1  Statistic  Models 

Log-logistic  distribution  (LLD)  [12]  is  a  special  case  of  Burr’s  type-XII  distribution  [13]  as  well  as  a 
special  case  of  the  kappa  distribution  proposed  by  Mielke  and  Jonson  [14].  Lee  et  al  employed  the 
LLD  for  frequency  analysis  of  multiyear  drought  durations  [15],  whereas  Shoukri  et  al  employed 
LLD  to  analyse  extensive  Canadian  precipitation  data  [16],  and  Narda  &  Malik  used  LLD  to 
develop  a  model  of  root  growth  and  water  uptake  in  wheat  [17].  In  spite  of  intensive  application  in 
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precipitation  and  stream-flow  data,  so  far  the  log-logistic  distribution  (LLD)  statistical  model  has 
never  been  applied  to  foliage  channel  model  to  the  best  our  knowledge. 

This  model  is  intended  to  be  employed  on  a  basis  of  higher  kurtosis  and  longer  tails,  as  well  as 
its  shape  similarity  to  log-normal  and  Weibull  distributions.  PDF  for  LLD  on  a  basis  of  different 
of  /i  and  a  are  illustrated  in  Fig.  13.  As  we  shall  see  that  this  model  provides  the  best  curve  fit 
compared  to  lognormal,  Weibull  and  Rayleigh. 

Here  we  use  the  two-parameter  distribution  with  parameters  fi  and  a.  The  PDF  for  this 
distribution  is  given  by 

lnx—n 

f(x)  =  — -C  —1^—,  x  >  0,  cr  >  0  (3) 

ax(l  +  e *  )2 

where  \i  is  scale  parameter  and  cr  is  shape  parameter.  The  mean  of  the  the  LLD  is 


E{x)  =  e»T(l  +  <t)T(1  -  a) 


(4) 


The  variance  is  given  by 


Var{x}  =  e2"{r(l  +  2*)r(l  -  2<r)  -  [T(l  +  *)r(l  -  a)]2} 


(5) 


while  the  moment  of  order  k  is 


where 


E{xk}  =  cre^B(ka)  1  —  kc r),  k  <  - 

(7 


B(m,  n)  =  f  xm  ldx 

Jo 


(6) 


(7) 


Similarly,  the  log-normal  distribution  [18]  is  a  two-parameter  distribution  with  parameters  [i 
and  cr.  The  PDF  for  this  distribution  is  given  by 


1  (lnx-ji)2 

f{x)  = - 7=e  2*2  ,  a:  >  0,  cr  >  0 

XCTV27T 


(8) 


where  /i  is  the  scale  parameter  and  <r  is  the  shape  parameter. 

The  Weibull  distribution,  which  is  named  after  Waloddi  Weibull,  can  be  made  to  fit  measure¬ 
ments  that  lie  between  the  Rayleigh  and  log-normal  distribution  [19]. 
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The  Weibull  distribution  is  also  a  two-parameter  distribution  with  parameters  a  and  b.  The 
PDF  for  this  distribution  is  given  by 

f(x)  =  6a_6x6_1e“^a^,  x  >  0,a  >  0,  b  >  0  (9) 

where  b  is  the  shape  parameter  and  a  is  the  scale  parameter. 

The  Rayleigh  distribution,  whose  real  and  imaginary  components  are  Gaussian,  has  the  PDF 
as  follows: 

/(*)  =  pe"£,  b>  0  (10) 

If  a  and  b  are  the  parameters  of  the  Weibull  distribution,  then  the  Rayleigh  distribution  with 
parameter  b  is  equivalent  to  the  Weibull  distribution  with  parameters  a  =  \/2b  and  b  =  2. 

4.2.2  Maximum  Likelihood  Estimation 

On  a  basis  of  CIR  cluster  amplitude  from  12  different  positions,  we  apply  Maximum  Likelihood 
Estimation  (MLE)  approach  to  estimate  the  parameters  for  log-logistic,  log-normal,  Weibull  and 
Rayleigh  models  respectively.  MLE  is  often  used  when  the  sample  data  are  known  and  parameters 
of  the  underlying  probability  distribution  are  to  be  estimated  [20]  [21].  It  is  generalized  as  follows: 

Let  2/1,  y2,  •  •  • ,  Vn  be  N  independent  samples  drawn  from  a  random  variable  Y  with  m  para¬ 
meters  0i,  02,  •  •  • ,  0m,  where  0*  £  0,  then  the  joint  PDF  of  2/1,  2/2,  •  *  * ,  UN  is 

Ln(  Y|0)  =  /y|0(2/i|01,02,***  ,0m)/y|0(2/2|01,02,  •  •  •  ,^m)-*-/Y|0(yn|01,^2,-**  ,0m)  (11) 

When  expressed  as  the  conditional  function  of  Y  depends  on  the  parameter  0,  the  likelihood 
function  is 


N 

Ln{ Y|0)  =  JJ  fY\o(yk\0\,02,  •  •  •  ,0m)  (12) 

k= 1 

The  maximum  likelihood  estimate  of  0i,  02,  •  ■  • ,  0m  is  the  set  of  values  0i,  02,  *  •  • ,  0m  that  maximize 
the  likelihood  function  L^(Y\9). 

As  the  logarithmic  function  is  monotonically  increasing,  maximizing  L^(Y|0)  is  equivalent  to 
maximizing  ln(Z,;v(Y|0)).  Hence,  it  can  be  shown  that  a  necessary  but  not  sufficient  condition  to 
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obtain  the  ML  estimate  6  is  to  solve  the  likelihood  equation 

^ln(LN(Y\0))  =  O  (13) 

We  obtain  /}  and  a  for  log-logistic,  ji  and  a  for  log-normal,  a  and  b  for  weibull  and  b  for  Rayleigh 
respectively,  which  are  shown  in  table  2.  Note  that  due  to  the  very  small  amount  of  channel  sample 
of  600MHz,  its  analysis  have  to  be  ignored.  We  also  explore  the  standard  deviation  (STD)  error  of 
each  parameter.  These  descriptions  are  also  shown  in  table  ??  in  the  form  of  ex ,  where  x  denotes 
different  parameter  for  each  model.  It  is  obvious  that  log-logistic  model  provides  smaller  STD 
errors  than  those  of  log-normal. 


4.2.3  Goodness-of-fit  in  curve  and  RMSE 

We  may  also  observe  that  to  what  extend  does  the  PDF  curve  of  the  statistic  model  match  that 
of  CIR  cluster  amplitude  by  root  mean  square  error  (RMSE).  Let  i  (i=l,2,  •  •  •  ,n)  be  the  sample 
index  of  CIR  amplitude,  c*  is  the  corresponding  PDF  value  whereas  c*  is  the  PDF  value  of  the 
statistical  model  with  estimated  parameters  by  means  of  MSE.  RMSE  is  obtained  through 


RMSE  = 


1=1 


(14) 


where  n  is  the  amount  of  sample  index.  The  RMSE  for  200MHz,  400Hz  and  UWB  have  been  listed 
in  Table  3.  It  demonstrates  that  LLD  turns  out  to  be  the  model  that  fit  the  channel  data  best. 

One  may  also  draw  the  above  conclusion  from  the  Fig.  14,  which  describes  the  the  goodness- 
of-fit  in  curve.  The  absolute  amplitude  of  clusters  have  been  plotted  in  terms  of  histogram.  It  can 
be  easily  seen  that  Rayleigh  model  provide  the  worst  goodness-of-fit  compared  to  LLD,  log-normal 
and  Weibull,  so  that  exponential  PDP  of  the  clusters  adopted  in  S-V  model  can  not  be  applied 
in  outdoor  NLOS  environment.  Also,  Weibull  is  not  a  good  choice  due  to  the  inaccurate  kurtosis 
and  high  tails.  Compare  LLD  with  log-normal,  it  is  obvious  that  LLD  is  able  to  provide  shaper 
kurtosis,  shaper  slope, and  lower  tail.  In  other  word,  LLD  provides  better  goodness-of-fit  than  that 
of  log-normal. 

Since  the  above  investigations  have  shown  that  LLD  can  better  characterize  the  multipath  gain 
amplitude  for  outdoor  NLOS  narrowband  and  UWB  channels  than  log-normal,  we  may  suggest  an 
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additional  LLD  variable  to  represent  the  fluctuations  of  the  total  multi-path  gain.  Moreover,  similar 
to  the  IEEE  UWB  model,  the  phase  of  each  path  in  outdoor  NLOS  environment  may  assumed  to 
be  either  0  or  7r  with  equal  probability. 

5  Conclusion 

In  our  investigation,  we  accomplished  following  conclusions:  1)  Outdoor  NLOS  channels  are  fre¬ 
quency  dependent  as  intervening  materials  have  dielectric  properties  that  are  strongly  frequency 
dependent.  2)  Both  narrowband  and  UWB  channels  are  made  up  of  multipath  and  time-varying 
components.  3)  In  outdoor  NLOS  environment,  the  largest  scattering,  i.e.,  the  highest  magnitude 
does  not  always  appear  at  the  first  path.  4)  The  outdoor  UWB  channels  we  obtained  are  similar 
in  their  basic  features  as  indoor  models.  However,  they  are  different  in  the  arrival  of  cluster  and 
subsequent  rays  as  well  as  envelope  decay.  5)  The  amplitude  of  channel  coefficient  at  each  path  can 
be  more  accurately  characterized  by  log-logistic  distribution  (LLD)  other  than  log-normal,  Weibull 
or  Rayleigh  due  to  the  best  goodness-of-fit  and  smallest  root- mean-square  (RMS). 
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Table  1:  Temporal  Parameters  for  Channel  Models 


Scenario 

A(l/ns) 

X(l/ns) 

200MHz 

0.012 

0.4 

400MHz 

0.004 

0.128 

600MHz 

0.002 

0.06 

Outdoor  UWB 

0.04 

0.8 

Indoor  UWB  Extreme  NLOS 

0.0667 

2.1 

Table  2:  Estimated  Parameters  for  Statistic  Model 


PDF 

Log-Logistic 

Log-normal 

Weibull 

Rayleigh 

200MHz 

A  =  -3.79907 
a  =  0.43948 
eM  =  0.0517626 
=  0.0250518 

A  =  -3.69473 
<t  =  0.811659 
e„  =  0.0550099 

E<j  =  0.0390963 

a  =  0.0388139 
b  =  1.00543 

Ea  =  0.0027934 

Eb  =  0.00456447 

6=0.0474046 

400MHz 

A  =  -3.75666 
a  =  0.482505 
e, „  =  0.071783 

Ea  =  0.035901 

A  =  -3.61265 
a  =  0.917049 
=  0.0795182 

Ea  =  0.0565477 

a  =  0.0447926 
b  =  0.903163 

Ea  =  0.00458706 

Eb  =  0.0536079 

6=0.0609159 

Outdoor  UWB 

A  =  -3.30616 
a  =  0.590192 
=  0.202988 

Ea  =  0.101636 

A  =  -3.13344 
<7  =  1.12623 
=  0.225245 

Ea  =  0.164277 

a  =  0.080002 
b  =  0.765597 

Ea  =  0.0222858 

Eb  =  0.106023 

6=0.141188 

Table  3:  root  mean  square  error  (RMSE)  comparison  between  Statistic  Models 


PDF 

Log-Logistic 

Log-normal 

Weibull 

Rayleigh 

200MHz 

5.7016 

6.2850 

8.8810 

9.7562 

400MHz 

5.9023 

6.5635 

9.7056 

10.3359 

UWB 

2.1867 

2.4756 

3.0136 

4.8975 
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Figure  1:  This  figure  shows  the  lift  with  the  experiment.  The  antennas  are  at  the  far  end  of  the 
lift  from  the  viewer  under  the  roof  that  was  built  to  shield  the  equipment  from  the  elements.  This 
picture  was  taken  in  September  with  the  foliage  largely  still  present.  The  cables  coming  from  the 
lift  are  a  ground  cable  to  an  earth  ground  and  one  of  4  tethers  used  in  windy  conditions. 
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Figure  2:  Data  file  structure 
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Figure  3:  Measurement  of  200MHz  and  35  pulses  average  :  (a)  transmitted  pulse  (b)  received 
echoes 
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Figure  4:  Measurement  of  400MHz  and  35  pulses  average:  (a)  transmitted  pulse  (b)  received  echoes 
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Figure  5:  Measurement  of  600MHz  and  35  pulses  average:  (a)  transmitted  pulse  (b)  received  echoes 
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Tranmistted  UWB  signal 


(b) 


Figure  6:  Measurement  of  UWB  and  35  pulses  average:  (a) transmitted  pulse  (b)  received  echoes 
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Figure  7:  200MHz  Channel 


Figure  8:  400MHz  Channel 
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Figure  9:  600MHz  Channel 


Figure  10:  UWB  Channel 
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Path  Magnitude 


Figure  11:  An  illustration  of  the  channel  impulse  in  S-V  model. 


Power  Delay 
Profile 


Figure  12:  An  illustration  of  the  double  exponential  decay  of  the  mean  cluster  power  and  the  ray 
power  within  clusters  in  S-V  model. 
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Figure  13:  An  illustration  of  log-logistic  PDF. 
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Figure  14:  Goodness-of-fit  (a)200MHz  (b)400MHz  (c)UWB 
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Outdoor  Propagation  Channel  Modeling  in  Foliage 

Environment 

Jing  Liang  and  Qilian  Liang 
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Abstract — In  this  paper,  we  study  the  statistical  modeling  for 
outdoor  non  line-of-sight  (NLOS)  channels  in  rich  scattering 
and  time-varying  foliage  environment  based  on  extensive  data 
collected  by  both  narrowband  and  ultra-wideband  (UWB)  radars. 
The  multi-path  contributions  arrive  at  the  receiver  are  grouped 
into  clusters.  The  time  of  arrival  of  clusters  can  be  modeled  as 
a  Poisson  arrival  process,  while  within  each  cluster,  subsequent 
multipath  contributions  or  rays  also  arrive  according  to  a  Poisson 
process.  However,  the  parameters  are  quite  different  along  with 
the  frequency.  We  also  observe  that  the  amplitude  of  channel 
coefficient  at  each  path  can  be  more  accurately  characterized  as 
log-logistic  distribution  (LLD)  other  than  log-normal,  Weibull  or 
Rayleigh  due  to  the  best  goodness-of-fit  and  smallest  root-mean- 
square  (RMS). 

I.  Introduction 

The  indoor  radio  propagation  channel  measurement  and 
modeling  has  been  intensively  investigated  [l]-[3].  If  the 
characterization  of  propagation  channel  is  known,  efforts  can 
be  employed  at  both  transmitters  and  receivers  to  overcome 
the  channel  effects  and  therefore  improve  the  communication 
performance.  As  far  as  military  warefare  is  concerned,  lots  of 
applications  are  for  outdoor  activities  such  as  sense-through- 
foliage  target  detection.  It  is  believed  that  understanding  the 
characterization  of  outdoor  non  line-of-sight  (NLOS)  channels 
will  assist  both  target  detection  as  well  as  commercial  wireless 
communications. 

Like  indoor  environment,  the  foliage  contains  a  wealth  of 
multiple  scattering  other  than  LOS  free  space.  In  addition, 
the  movement  of  leaves,  branches  and  even  the  tree  trunks 
contribute  to  the  time-variance  fading  phenomenon.  Since  the 
foliage  medium  can  be  completely  described  by  its  time  and 
space  varying  feature,  one  can  investigate  the  channel  model 
based  on  characterization  of  the  channel  impulse  response 
(CIR). 

In  our  investigation,  we  will  apply  both  narrowband  and 
ultra-wide  band  (UWB)  radars  to  model  the  propagation 
channels,  as  we  believe  that  foliage  is  composed  of  interven¬ 
ing  materials  that  are  electromagnetically  dispersive,  which 
contributes  to  the  strong  frequency  dependence  of  foliage, 
and  thus  a  narrowband-wideband  study  would  assist  with 
the  better  understanding  of  statistic  property  of  the  channel. 
Narrowband  signals  have  been  tried  at  200,  400  and  600 
megahertz  respectively,  while  UWB  radar  emissions  are  at  a 
relatively  low  frequency-typically  between  100  MHz  and  3 
GHz.  Each  frequency  component  in  a  radar  signal  will  sense 


the  foliage  clutter  in  a  slightly  different  manner,  therefore 
provides  differences  in  multipath. 

The  rest  of  this  paper  is  organized  as  follows.  In  Section 
II,  we  summarize  the  measurement  and  collection  of  the 
data.  In  Section  HI,  we  apply  CLEAN  algorithm  to  extract 
CIR  for  200MHz,  400MHZ,  600MHz  and  UWB  signals. 
Section  IV  presents  the  channel  model  in  view  of  temporal 
characterization  as  well  as  statistic  model  comparison.  We 
conclude  our  work  in  Section  V. 

II.  Measurement  Setup 

The  foliage  penetration  measurement  effort  began  in  Au¬ 
gust  2005  and  continued  through  December  2005.  Working 
in  August  through  the  fall  of  2005,  the  foliage  measured 
included  late  summer  foliage  and  fall  and  early  winter  foliage. 
Late  summer  foliage,  due  to  the  limited  rainfall,  involved 
foliage  with  decreased  water  content.  Late  fall  and  winter 
measurements  involved  largely  defoliated  but  dense  forest, 
so  it’s  a  rich  scattering  environment.  Because  of  wind  or 
different  temperature  in  dense  forest,  it’s  also  a  time- varying 
environment. 

The  radar  experiment  was  constructed  on  a  seven-ton  man 
lift,  which  had  a  total  lifting  capacity  of  450  kg.  The  limit 
of  the  lifting  capacity  was  reached  during  the  experiment  as 
essentially  the  entire  measuring  apparatus  was  placed  on  the 
lift  (as  shown  in  Fig.  1).  The  principle  pieces  of  equipment 
secured  on  the  lift  are:  dual  antenna  mounting  stand;  two 
antennas;  200MHz,  lKw  Amplifier,  power  supply,  pre-amp; 
400MHz,  lKw  Amplifier,  power  supply,  pre-amp;  600MHz, 
lKw  Amplifier,  power  supply,  pre-amp;  Tektronix  model  7704 
B  oscilloscope;  rack  system;  IBM  laptop;  HP  signal  Generator; 
Custom  RF  switch  and  power  supply  and  weather  shield  (small 
hut).  Particularly  for  UWB  signal,  a  Barth  pulse  source  (Barth 
Electronics,  Inc.  model  732  GL)  was  used.  The  pulse  generator 
uses  a  coaxial  reed  switch  to  discharge  a  charge  line  for  a  very 
fast  rise  time  pulse  outputs.  The  model  732  pulse  generator 
provides  pulses  of  less  than  50  picoseconds  (ps)  rise  time, 
with  amplitude  from  150  V  to  greater  than  2  KV  into  any 
load  impedance  through  a  50  ohm  coaxial  line.  The  generator 
is  capable  of  producing  pulses  with  a  minimum  width  of  750 
ps  and  a  maximum  of  1  microsecond.  This  output  pulse  width 
is  determined  by  charge  line  length  for  rectangular  pulses,  or 
by  capacitors  for  1/e  decay  pulses. 

The  system  was  pointing  at  the  specified  250  feet  dis¬ 
tance.  For  the  data  we  used  in  this  paper,  each  sample  is 
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Fig.  1.  This  figure  shows  the  lift  with  the  experiment.  The  antennas  are  at 
the  far  end  of  the  lift  from  the  viewer  under  the  roof  that  was  built  to  shield 
the  equipment  from  the  elements.  This  picture  was  taken  in  September  with 
the  foliage  largely  still  present.  The  cables  coming  from  the  lift  are  a  ground 
cable  to  an  earth  ground  and  one  of  4  tethers  used  in  windy  conditions. 


spaced  at  50  picoseconds  interval,  and  16,000  samples  were 
collected  for  each  collection  for  a  total  time  duration  of 
0.8  microseconds.  The  accomplished  data  structure  is  shown 
in  Fig.  2.  For  each  frequency  band,  12  positions  named 

and  ’20’  have  been 

tested.  35  pulses  have  been  obtained  for  the  transmitted  and 
received  signal  respectively  at  each  position.  *EMCO’  is  the 
brand  name  of  the  antenna  horn.  ’Background’  refers  to  the 
background  data  taken  along  the  track  with  the  antenna  bore 
sighted  to  the  targets. 


EMCO 


[ 


Background 


|  200MHz  |  |  400MHz  |  |  60CMHz  |  |  UWB  | 

CgD  [jp  Hn  i~s~i  r~en  i~e~i  non  rfen  nVi  rfei  nai  rzoi 

Target  |  |  No  Target  | 


|  Transmit  |  |  Receive  |  [Transmit  |  [  Receive 


Fig.  2.  Data  file  structure 

In  the  foliage-target  data,  we  fired  from  each  of  twelve 
positions  along  the  track.  The  firing  was  into  the  foliage  at 
right  angles  to  the  track  and  these  right  angle  shots  are  called 
track  data,  i.e.,  ’Track’.  At  all  twelve  positions  along  the  track, 
the  azimuth  angle  of  the  transmit  antenna  was  changed  from 
plus  10  to  -  10  degrees  to  provide  radial  measurements,  i.e., 
’Radial’. 

III.  Channel  Impulse  Response  Based  on  the 
Measured  Data  and  CLEAN  Algorithm 

The  average  transmitted  and  received  pulses  at  position  4  for 
different  frequencies  have  been  illustrated  from  Figs.  3-6.  The 
purpose  of  average  is  to  remove  white  Gaussian  noise.  Note 


that  at  a  different  position  the  result  will  be  slightly  different. 
However,  illustration  at  one  position  is  sufficient  enough  to 
describe  the  characterization. 


0  2000  4000  6000  8000  10000  12000  14000  16000 

Time  sample  index 

(b) 


Fig.  3.  Measurement  of  200MHz  and  35  pulses  average  :  (a)  transmitted 
pulse  (b)  received  echoes 

The  complicated  multipath  and  time-varying  CIR  can  be 
modeled  as  follows  [4] 

r(t)  ~  ^  ^  anPn(t  ~  Tn)  (1) 

n 

where  an  and  rn  is  referred  to  as  the  amplitude  and  delay  of 
the  nth  propagation  path.  In  order  to  extract  the  CIR  from  our 
measurement,  the  CLEAN  algorithm  has  been  used.  It  was 
initially  introduced  in  [5]  to  enhance  radio  astronomical  maps 
of  the  sky,  and  has  also  been  employed  in  narrowband  channel 
modeling  [6][7]and  UWB  channel  characterization  problems 
[4][8]. 

Our  steps  involved  [8]  are: 

1)  Calculate  the  autocorrelation  of  the  transmitted  signal 
Rss{t)  and  the  cross-correlation  of  the  transmitted  with 
the  received  waveform  Rsy(t). 

2)  Find  the  largest  correlation  peak  in  Rsy(t)y  record  the 
normalized  amplitudes  a*  and  the  relative  time  delay  rk 
of  the  correlation  peak. 

3)  Subtract  Rss{t)  scaled  by  ak  from  Rsy{t)  at  the  time 
delay  rk. 
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Fig.  4.  Measurement  of  400MHz  and  35  pulses  average:  (a)  transmitted  Fig.  5.  Measurement  of  600MHz  and  35  pulses  average:  (a)  transmitted 
pulse  (b)  received  echoes  pulse  (b)  received  echoes 


4)  If  a  stopping  criterion  (a  minimum  threshold)  on  the 
peak  correlation  is  not  met,  go  step  2.  otherwise  stop. 

Given  the  transmission,  reception  and  the  CLEAN  process¬ 
ing  described  above,  the  obtained  CIR  are  illustrated  from  Fig. 
7  to  Fig.  10.  Note  that  we  plot  the  absolute  value  of  the  UWB 
channel  for  the  comparison  between  the  outdoor  UWB  channel 
with  the  indoor  S-V  model  [9]  (see  Fig.  11).  It  is  shown  that 

1)  Both  narrowband  and  UWB  channels  are  made  up  of 
multipath  components.  The  time-varying  path  magnitude 
implies  that  fading  generally  exists,  therefore  the  re¬ 
ceived  pulses  are  fairly  random  from  one  time  to  another. 

2)  The  UWB  channel  we  obtained  looks  similar  as  the 
CIR  in  S-V  model  shown  in  Fig.  1 1 .  However,  they  are 
different  in  the  arrival  of  cluster,  subsequent  rays  and 
envelope  decay. 

3)  In  outdoor  environment,  the  largest  scattering,  i.e.,  the 
highest  magnitude  does  not  always  appear  at  the  first 
path.  This  phenomenon  is  clearly  illustrated  in  Fig.  7 
and  10. 

4)  Channels  are  frequency  dependent.  It  has  been  observed 
that  the  intervening  materials,  such  as  foliage  and  soil, 
have  dielectric  properties  that  are  strongly  frequency 
dependent.  This  in  part  explains  the  difference  among 
those  channels.  We  will  further  analyze  the  detail  in  the 
following  section. 


IV.  Outdoor  Channel  Modeling 

A.  Temporal  Characterization 

In  the  S-V  model,  the  arrival  of  clusters  is  modelled  as 
a  Poisson  arrival  process  with  a  rate  A,  while  within  each 
cluster,  subsequent  multipath  contributions  or  rays  also  arrive 
according  to  a  Poisson  process  with  a  rate  A  (see  Fig.  12). 
Observe  Fig.  7~10,  like  in  S-V  model,  multipath  contributions 
arrive  at  the  receiver  grouped  into  clusters  and  therefore 
similar  methodology  used  in  S-V  model  studies  may  be  also 
applied  to  200MHz,  400MHz,  600MHz  and  UWB  CIR.  The 
time  of  arrival  of  clusters  can  be  modeled  as  a  Poisson  arrival 
process  with  a  rate  A,  while  within  each  cluster,  subsequent 
multipath  contributions  or  rays  also  arrive  according  to  a 
Poisson  process  with  a  rate  A. 

We  define: 

•  Ti  =  the  arrival  time  of  the  first  path  of  the  Z-th  cluster; 

•  Tk,i  =  the  delay  of  the  k- th  path  within  the  Z-th  cluster 
relative  to  the  first  path  arrival  time  7); 

•  A  =  the  cluster  arrival  rate; 

•  A  =  the  ray  arrival  rate,  i.e.,  the  arrival  rate  of  the  paths 
within  each  cluster. 

By  definition,  we  have  tqi  =  T}.  The  distributions  of  the 
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(a) 


(b) 


Fig.  6.  Measurement  of  UWB  and  35  pulses  average:  (a)transmitted  pulse 
(b)  received  echoes 


TABLE  I 

Temporal  Parameters  for  Channel  Models 


Scenario 

A(l/ns) 

A(1  /ns) 

200MHz 

0.012 

0.4 

400MHz 

0.004 

0.128 

600MHz 

0.002 

0.06 

Outdoor  UWB 

0.04 

0.8 

Indoor  UWB  Extreme  NLOS 

0.0667 

2.1 

cluster  arrival  time  and  the  ray  arrival  time  are  given  by 

p(Ti\Ti-i)  =  Aexp(— A(7|  >  0 

P('rfc,il'r(fc— i),i)  =  Aexp  (-A (tm  -  T(fc_1M)),  k  >  0  (2) 

Specifically,  we  also  observed  that  the  A  and  A  are  quite 
different  for  200MHz,  400MHz,  600MHz  and  UWB  CIR.  We 
listed  observed  parameters  in  Table  I.  As  for  indoor  UWB 
data,  we  refer  [11].  The  higher  A  and  A  of  UWB  implies  its 
exceptional  range  resolution.  Lower  A  and  A  of  outdoor  UWB 
than  those  of  indoor  means  outdoor  environment  typically 
more  sparse  multiple  scattering  than  that  of  indoor. 

B.  Statistical  Distribution  of  Channel  Amplitude 

In  the  S-V  model,  the  average  Power  Decay  Profile  (PDP) 
is  characterized  by  an  exponential  decay  of  the  clusters  and  a 
different  exponential  decay  for  the  pulses  within  each  cluster. 


Fig.  7.  200MHz  Channel 


Fig.  8.  400MHz  Channel 


In  other  words,  the  amplitude  follows  rayleigh  distribution. 
In  the  IEEE  UWB  indoor  channel  model  [10],  the  clutter 
approach  was  adopted  (same  as  S-V  model),  but  a  log-normal 
distribution  was  suggested  for  characterizing  the  multipath 
gain  amplitude,  and  an  additional  log-normal  variable  was  in¬ 
troduced  for  representing  the  fluctuations  of  the  total  multipath 
gain.  In  this  Section  we  propose  that  log-logistic  model  may 
better  characterize  amplitude  of  the  multipath  for  both  outdoor 
NLOS  narrowband  and  UWB  signals. 

1)  Statistic  Models:  Log-logistic  distribution  (LLD)  [12]  is 
a  special  case  of  Burr’s  type-XII  distribution  [13]  as  well  as 
a  special  case  of  the  kappa  distribution  proposed  by  Mielke 
and  Jonson  [14].  Lee  et  al.  employed  the  LLD  for  frequency 
analysis  of  multiyear  drought  durations  [15],  whereas  Shoukri 
et  al.  employed  LLD  to  analyse  extensive  Canadian  precipi¬ 
tation  data  [16],  and  Narda  &  Malik  used  LLD  to  develop  a 
model  of  root  growth  and  water  uptake  in  wheat  [17].  In  spite 
of  intensive  application  in  precipitation  and  stream-flow  data, 
so  far  LLD  statistical  model  has  never  been  applied  to  foliage 
channel  model  to  the  best  our  knowledge. 

This  model  is  intended  to  be  employed  on  a  basis  of  higher 
kurtosis  and  longer  tails,  as  well  as  its  shape  similarity  to  log¬ 
normal  and  Weibull  distributions.  PDF  for  LLD  on  a  basis  of 
different  of  p  and  a  are  illustrated  in  Fig.  13.  As  we  shall 
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Fig.  9.  600MHz  Channel 


Fig.  12.  An  illustration  of  the  double  exponential  decay  of  the  mean  cluster 
power  and  the  ray  power  within  clusters  in  S-V  model. 


Fig.  13.  An  illustration  of  log-logistic  PDF. 


Fig.  10.  UWB  Channel 


see  that  this  model  provides  the  best  curve  fit  compared  to 
log-normal,  Weibull  and  Rayleigh. 

Here  we  use  the  two-parameter  distribution  with  parameters 
/z  and  a.  The  PDF  for  this  distribution  is  given  by 


The  variance  is  given  by 

Var{x}  =  e2>1{r(l  +  2a)r(l  -  2<r)  -  [r(l  +  <r)r(l  -  cr)]2} 

(5) 

while  the  moment  of  order  k  is 


E{xk}  =  cre^B^cr,  1  —  fccr),  k  <  — 


1 


f(x)  = 


6  "t a-a.  ,  *  >  0,  <T  >  0  (3) 

crx(l  +  e  *  y 

where  /z  is  scale  parameter  and  cr  is  shape  parameter.  The 
mean  of  the  the  LLD  is 


where 


B(m,n)  =  f  xm  1  (1  —  x)n  ldx 

Jo 


(6) 


(7) 


E{x}  =  e*ir(l+a)r(l-cr) 


(4) 


Similarly,  the  log-normal  distribution  [18]  is  a  two- 
parameter  distribution  with  parameters  /z  and  a.  The  PDF  for 
this  distribution  is  given  by 


Path  Magnitude 


f(x)  = 


1 


Ty/2n 


( lnx-u ) 
^2 


2°*  ,  x  >  0,  a  >  0 


(8) 


Fig.  11.  An  illustration  of  the  channel  impulse  in  S-V  model. 


where  /z  is  the  scale  parameter  and  a  is  the  shape  parameter. 

The  Weibull  distribution  can  be  made  to  fit  measurements 
that  lie  between  the  Rayleigh  and  log-normal  distribution  [19]. 

The  Weibull  distribution  is  also  a  two-parameter  distribution 
with  parameters  a  and  b.  The  PDF  for  this  distribution  is  given 
by 

f(x)  =  ba-bx6-1e”(x/o)\  x  >  0,  a  >  0,  b  >  0  (9) 

where  b  is  the  shape  parameter  and  a  is  the  scale  parameter. 


Page  118  of  134 


The  Rayleigh  distribution,  whose  real  and  imaginary  com¬ 
ponents  are  Gaussian,  has  the  PDF  as  follows: 

f(x)  =  £e~&,  b>  0  (10) 

2)  Maximum  Likelihood  Estimation:  On  a  basis  of  CIR 
cluster  amplitude  from  12  different  positions,  we  apply  Max¬ 
imum  Likelihood  Estimation  (MLE)  approach  to  estimate  the 
parameters  for  log-logistic,  log-normal,  Weibull  and  Rayleigh 
models  respectively.  MLE  is  often  used  when  the  sample 
data  are  known  and  parameters  of  the  underlying  probability 
distribution  are  to  be  estimated  [20]  [21].  It  is  generalized  as 
follows: 

Let  2/i>  2/2*  •  *  * »  2/n  be  N  independent  samples  drawn  from  a 
random  variable  Y  with  m  parameters  0i,  02,  •  •  • ,  0m,  where 
0j  €  0,  then  the  joint  PDF  of  yi,  2/2,  *  *  * ,  Vn  is 

Ln{ Y|0)  =  /y|0(j/l|01,02,  •  •  *  >0m)  •  •  •  /y|fl(yn|01,02,  •  •  •  ,0m) 

(ID 

When  expressed  as  the  conditional  function  of  Y  depends 
on  the  parameter  0,  the  likelihood  function  is 

N 

Ln(y\o)  =  Y\  fY\$(yk\0i,02,---  ,0m)  (12) 

k= 1 

The  maximum  likelihood  estimate  of  0X,  02,  •  •  • ,  0m  is  the  set 
of  values  0i,  02,  •  •  • ,  0m  that  maximize  the  likelihood  function 
Ln{Y\0). 

As  the  logarithmic  function  is  monotonically  increas¬ 
ing,  maximizing  Ljy(Y|0)  is  equivalent  to  maximizing 
ln(Z/yv (Y|0)).  Hence,  it  can  be  shown  that  a  necessary  but  not 
sufficient  condition  to  obtain  the  ML  estimate  0  is  to  solve  the 
likelihood  equation 

^ln(Ljv(Y|fl))=0  (13) 

We  obtain  jl  and  <r  for  log-logistic,  /}  and  a  for  log-normal, 
a  and  b  for  weibull  and  b  for  Rayleigh  respectively,  which  are 
shown  in  table  II.  Note  that  due  to  the  very  small  amount  of 
channel  sample  of  600MHz,  its  analysis  have  to  be  ignored. 
We  also  explore  the  standard  deviation  (STD)  error  of  each 
parameter.  These  descriptions  are  also  shown  in  table  II  in 
the  form  of  ex,  where  x  denotes  different  parameter  for  each 
model.  It  is  obvious  that  log-logistic  model  provides  smaller 
STD  errors  than  those  of  log-normal. 

3)  Goodness-of-fit  in  curve  and  RMSE:  We  may  also  ob¬ 
serve  that  to  what  extend  does  the  PDF  curve  of  the  statistic 
model  match  that  of  CIR  cluster  amplitude  by  root  mean 
square  error  (RMSE).  Let  i  (i=l,  2,  •  •  •  ,  n)  be  the  sample  index 
of  CIR  amplitude,  c*  is  the  corresponding  PDF  value  whereas 
Ci  is  the  PDF  value  of  the  statistical  model  with  estimated 
parameters  by  means  of  MSE.  RMSE  is  obtained  through 


RMSE  = 


\ 


t=l 


(14) 


where  n  is  the  amount  of  sample  index.  The  RMSE  for 
200MHz,  400Hz  and  UWB  have  been  listed  in  Table  III.  It 


TABLE  II 

Estimated  Parameters  for  Statistic  Model 


PDF 

Log-Logistic 

Log-normal 

200MHz 

A  =  —3.79907 
<7  =  0.43948 
=  0.0517626 
e*  =  0.0250518 

A  =  -3.69473 
<7  =  0.811659 
e n  =  0.0550099 
=  0.0390963 

400MHz 

A  =  -3.75666 

<7  =  0.482505 
=  0.071783 
e*  =  0.035901 

A  =  -3.61265 
a  =  0.917049 
=  0.0795182 
£a  =  0.0565477 

Outdoor  UWB 

A  =  -3.30616 
<7  =  0.590192 
=  0.202988 
£<r  =  0.101636 

A  =  -3.13344 
<7=  1.12623 
£»  =  0.225245 
£a  =0.164277 

PDF 

Weibull 

Rayleigh 

200MHz 

a  =  0.0388139 

5  =  1.00543 
£a  =  0.0027934 
£b  =  0.00456447 

5=0.0474046 

400MHz 

a  =  0.0447926 

5  =  0.903163 
£a  =  0.00458706 
£b  =  0.0536079 

5=0.0609159 

Outdoor  UWB 

a  =  0.080002 

6  =  0.765597 
£a  =  0.0222858 
£b  =  0.106023 

5=0.141188 

demonstrates  that  LLD  turns  out  to  be  the  model  that  fits  the 
channel  data  best. 


TABLE  ill 

ROOT  MEAN  SQUARE  ERROR  (RMSE)  COMPARISON  BETWEEN  STATISTIC 

Models 


PDF 

Log-Logistic 

Log-normal 

Weibull 

Rayleigh 

200MHz 

5.7016 

6.2850 

8.8810 

9.7562 

400MHz 

5.9023 

6.5635 

9.7056 

10.3359 

UWB 

2.1867 

2.4756 

3.0136 

4.8975 

One  may  also  draw  the  above  conclusion  from  the  Fig.  14, 
which  describes  the  the  goodness-of-fit  in  curve.  The  absolute 
amplitude  of  clusters  have  been  plotted  in  terms  of  histogram. 
It  can  be  easily  seen  that  Rayleigh  model  provides  the  worst 
goodness-of-fit  compared  to  LLD,  log-normal  and  Weibull,  so 
that  exponential  PDP  of  the  clusters  adopted  in  S-V  model  can 
not  be  applied  in  outdoor  NLOS  environment.  Also,  Weibull 
is  not  a  good  choice  due  to  the  inaccurate  kurtosis  and  high 
tails.  Compare  LLD  with  log-normal,  it  is  obvious  that  LLD 
is  able  to  provide  shaper  kurtosis,  shaper  slope,  and  lower  tail. 
In  other  word,  LLD  provides  better  goodness-of-fit  than  that 
of  log-normal. 

Since  the  above  investigations  have  shown  that  LLD  can 
better  characterize  the  multipath  gain  amplitude  for  outdoor 
NLOS  narrowband  and  UWB  channels  than  log-normal,  we 
may  suggest  an  additional  LLD  variable  to  represent  the 
fluctuations  of  the  total  multi-path  gain.  Moreover,  similar  to 
the  IEEE  UWB  model,  the  phase  of  each  path  in  outdoor 
NLOS  environment  may  assumed  to  be  either  0  or  n  with 
equal  probability. 
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different  in  the  arrival  of  cluster  and  subsequent  rays  as  well 
as  envelope  decay.  5)  The  amplitude  of  channel  coefficient  at 
each  path  can  be  more  accurately  characterized  as  log-logistic 
distribution  (LLD)  other  than  log-normal,  Weibull  or  Rayleigh 
due  to  the  best  goodness-of-fit  and  smallest  root-mean-squarc 
(RMS). 
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Fig.  14.  Goodness-of-fit  (a)200MHz  (b)400MHz  (c)UWB 


V.  Conclusion 

In  our  investigation,  we  accomplished  following  conclu¬ 
sions:  1)  Outdoor  NLOS  channels  are  frequency  dependent 
as  intervening  materials  have  dielectric  properties  that  are 
strongly  frequency  dependent.  2)  Both  narrowband  and  UWB 
channels  are  made  up  of  multipath  and  time-varying  compo¬ 
nents.  3)  In  outdoor  NLOS  environment,  the  largest  scattering, 
i.e.,  the  highest  magnitude  does  not  always  appear  at  the  first 
path.  4)  The  outdoor  UWB  channels  we  obtained  are  similar 
in  their  basic  features  as  indoor  models.  However,  they  are 
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I.  Abstract 

This  paper  presents  new  developed  code  —  punctured  bi¬ 
nary  sequence -pair.  The  definitions  and  the  autocorrelation 
properties  of  the  proposed  code  are  given.  Doppler  shift 
performance  is  also  investigated.  The  significant  advantages 
of  punctured  binary  sequence-pair  over  conventional  pulse 
compression  codes,  such  as  the  widely  used  Barker  codes, 
are  zero  autocorrelation  side  lobes  and  the  longer  length  of 
the  code  which  can  be  as  long  as  31  so  far.  In  the  radar 
target  detection  system  simulation,  punctured  binary  sequence- 
pair  also  outperforms  other  conventional  pulse  compression 
codes.  Therefore,  our  proposed  code  can  be  used  as  one  of 
the  candidates  for  pulse  compression  code. 

II.  Introduction 

Pulse  compression,  which  allows  a  radar  to  simultaneously 
achieve  the  energy  of  a  long  pulse  and  the  resolution  of  a 
short  pulse  without  the  high  peak  power  required  by  a  high 
energy  short  duration  pulse  [1],  is  generally  used  in  modem 
radar  system.  The  main  purpose  of  this  technique  is  to  raise 
the  signal  to  maximum  sidelobe  (signal-to-sidelobe)  ratio  to 
improve  the  target  detection  and  range  resolution  abilities  of 
the  radar  system.  The  range  sidelobes  are  harmful  because  they 
can  mask  main  peaks  caused  by  small  targets  situated  near 
large  targets.  The  lower  the  sidelobes,  relative  to  the  mainlobe 
peak,  the  better  the  main  peak  can  be  distinguished. 

In  the  industrial  world,  pulse  compression  is  one  of  the 
significant  factors  to  determine  the  performance  of  high  de¬ 
tection  and  high  resolution  radar.  For  example,  a  satellite- 
borne  rain  radar  demands  very  stringent  requirements  on 
range  sidelobe  level  of  -60dB  [2],  a  downward  looking  rain 
measuring  radar  requires  a  range  sidelobe  of  55dB  under 
the  mainlobe  level  [3]  [4],  and  the  air  traffic  control  system 
demands  the  sidelobe  level  lower  -55dB  [5]. 

There  are  two  kinds  of  basic  waveform  designs  suitable  for 
pulse  compression:  frequency-codes,  such  as  linear  frequency 
modulation(LFM)  codes  [6]  [7]  and  nonlinear  frequency  mod¬ 
ulation  codes(NLFM)  [7]  [8]  [9];  phase-coded  waveforms, 
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such  as  binary  phase  codes  and  polyphase  codes.  For  a  phase- 
coded  waveform,  a  long  pulse  of  duration  T  is  divided  into  N 
subpulses  each  of  width  Ts.  Each  subpulse  has  a  particular 
phase,  which  is  selected  in  accordance  with  a  given  code 
sequence.  And  the  pulse  compression  ratio  equals  to  the 
number  of  subpulses  N  =  T/Ts  «  BTy  where  the  bandwidth 
is  B  «  1/7V 

A  common  form  of  phase  coding  is  binary  phase  coding ,  in 
which  the  phase  of  each  subpulse  is  selected  to  be  either  0  or  n 
radians.  Since  the  binary  phase  codes  are  easy  to  generate,  they 
are  widely  used  in  modem  radar  system.  However,  when  the 
selection  of  the  phase  is  made  randomly,  the  expected  maxi¬ 
mum  sidelobe  is  only  about  2/N  of  the  peak  of  the  compressed 
pulse.  So  completely  random  selection  of  the  phase,  is  not  a 
good  idea,  and  the  criterion  for  selecting  the  subpulse  phases 
is  that  all  the  time-sidelobe  of  the  compressed  pulse  should 
be  equal  and  as  low  as  possible.  One  family  of  binary  phase 
code  widely  used  nowadays  that  can  produce  compressed 
waveforms  with  constant  sidelobe  levels  equal  to  unity  is  the 
Barker  code.  It  has  special  features  with  which  its  sidelobe 
structure  contains  the  minimum  energy  which  is  theoretically 
possible  for  binary  codes,  and  the  energy  is  uniformly  dis¬ 
tributed  among  the  sidelobes(the  sidelobe  level  of  the  Barker 
codes  is  l/N2  that  of  the  peak  signal)  [10].  Unfortunately,  the 
length  N  of  known  binary  and  complex  Barker  codes  is  limited 
to  13  and  25,  respectively  [11],  which  may  not  be  sufficient 
for  the  desired  radar  applications.  In  [12]  [13]  [14],  polyphase 
codes,  with  better  Doppler  tolerance  and  lower  range  sidelobes 
such  as  the  Frank  and  PI  codes,  the  Butler-matrix  derived 
P2  code  and  the  linear-frequency  derived  P3  and  P4  codes 
were  intensively  analyzed.  However,  the  low  range  sidelobe 
of  the  polyphase  codes  can  not  reach  the  level  zero  either. 
The  structure  of  polyphase  codes  is  more  complicated  and  is 
not  easy  to  generate  comparing  with  binary  codes.  Therefore, 
we  propose  and  analyze  a  new  kind  of  code-punctured  binary 
sequence-pair,  whose  sidelobe  level  is  as  low  as  0  and  the 
longest  length  of  which  is  3 1  in  this  paper,  and  subsequently 
apply  it  to  radar  system.  According  to  the  simulation  results, 
the  new  code  can  be  a  good  alternative  for  the  current  used 
pulse  compression  codes  in  radar  system. 
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The  rest  of  the  paper  is  organized  as  following.  Section  2 
introduces  the  basic  concept  and  properties  of  our  proposed 
code.  In  Section  3,  an  example  of  punctured  binary  sequence- 
pair  is  given  and  its  properties  are  investigated.  In  Section  4, 
the  performance  of  our  proposed  code  is  also  simulated  and 
analyzed  in  radar  targets  detection  system.  In  Section  5,  some 
conclusions  are  drawn  about  the  punctured  binary  sequence- 
pair. 

III.  Design  of  Punctured  Binary  Sequence-pair 

First  of  all,  there  are  some  relating  definitions  listed  here. 

Definition  1  A  sequence-pair  (x,  y)  is  made  up  of  two 
AMength  sequences  x  =  (xo,xi,  •  ■  •  ,xjv_i)  and  y  = 
(l/o,  J/i,  *  *  •  ,2/N-l), 

N- 1 

Rxy(r)  =  ^2  X jy (j+r)modN  1  ^  <  T  <  N  —  \  (1) 

7= o 


Theorem  1  Mapping  property,  if  xi(z)  =  x(-i),y1(i)  = 
y(— i),then  sequence-pair  (xi,yi)  is  perfect  punctured  binary 
sequence-pair. 

Theorem  2  Opposite  to  element  symbol  property,  if  x  i  (i)  = 
-x(i),yi(i)  =  -y(z),then  sequence-pair  (xi,yi)  is  perfect 
punctured  binary  sequence-pair. 

Theorem  3  Cyclic  shift  property,  if  xi(i)  =  -x(i  + 
tO,yi(i)  =  — y(i  +  t*), then  sequence-pair  (xi,yi)  is  perfect 
punctured  binary  sequence-pair. 

Theorem  4  Periodically  sampling  property,  if  xi(i)  = 
-x(ki),yi(i)  =  -y(hi),k  and  N  are  relatively  prime,  then 
sequence-pair  (xi,yi)  is  perfect  punctured  binary  sequence- 
pair. 

In  [15],  the  perfect  punctured  binary  sequence-pairs  of 
length  from  3  to  31  are  presented  in  the  Table  1. 


IV.  Properties 


is  called  the  periodic  autocorrelation  function  of  the  se¬ 
quence  pair,  while  x  =  y,  the  sequence-pair  (x,y)  turns  to 
be  a  one-sequence  code. 

Definition  2  [15]  Sequence  y  =  (yo,yi,  •  •  •  >2W_i)  is  the 
punctured  sequence  for  x  =  (x0, xi,  •  •  • , x^_x), 

{0  j  ep-punctured  bits 
Xj  j  GA on-punctured  bits 

Where  p  is  the  number  of  punctured  bits  in  sequence  jc, 
suppose  x  =  [ — 1, 1],  is  p-punctured  binary  sequence,  (x,y) 
is  called  a  punctured  binary  sequence-pair. 

Definition  3  The  periodic  autocorrelation  of  punctured 
sequence-pair  (x,y)  is  defined 


A.  Autocorrelation  Properties 

The  autocorrelation  function  is  one  of  the  most  important 
properties  that  represents  the  compressed  pulse  in  an  ideal 
pulse  compression  system,  because  it  is  proportional  to  the 
matched  filter  response  in  the  noise-free  condition.  As  pre¬ 
sented  in  the  equation  (4)  in  the  last  section,  the  periodic 
autocorrelation  function  of  the  punctured  binary  sequence-pair 
is 


N- 1 

Rxy(r)  =  ^2  xjy(j+r)modN  = 

7=o 


r  =  0  mod  N 
otherwise 


N- 1 

Rxy(x)  =  ^2  xjy\j+T)modN >  0  <  T  <  N  -  l 
7=o 

When  punctured  sequence-pair  has  the  following  autocor¬ 
relation  property 


R 


xy 


E 

0 


r  =  0  mod  N 
otherwise 


(3) 


it  is  called  perfect  punctured  sequence-pair  [15].  Here,  E  = 
]Cj=o  =  N  —  py  is  the  energy  of  punctured  sequence- 
pair.  Then  binary  sequence-pair  (x,y)  is  called  a  /^-punctured 
sequence-pair.  The  energy  efficiency  of  the  sequence-pair  is 
defined  as 


77  = 


E 

N 


N-p 

N 


(4) 


Definition  4  The  balance  of  the  sequence  x  is  defined  as 
7  =  XJ  —  np  ~  nn,  while  np,nn  are  the  number  of 

'  +  1'  and  '  -  1'  in  x  separately. 


EXAMPLE 

The  autocorrelation  property  of  13-length  punctured  binary 
sequence  pair(x,y),  (x  =  (111  -111-1-1-1-111-1) 
and  y  =  (1010110000110)),  and  that  of  31 -length  punctured 
binary  sequence-pair (x,  y),  (x  =  (1111  —  1  —  1  —  11  —  11  — 

1111  -1-1-1-11-1-11-1-  1111  -  111  -  1) 

and  y  =  (1111000101011100001001001110110))  are  shown 
in  the  Fig.l  and  Fig.2. 

As  it  is  known  that  a  suitable  criterion  for  evaluating  code 
of  length  N  is  the  peak  signal  to  peak  signal  sidelobe  ratio 
(PSR)  of  their  aperiodic  autocorrelation  function,  which  can 
be  bounded  by  [16] 

<  20 logN  =  [PSRmax }dD  (5) 

The  only  uniform  phase  codes  that  can  reach  the  PSRmax 
are  the  Barker  codes  whose  length  is  equal  or  less  than  13. 
However,  the  sidelobe  of  the  new  code  in  both  Fig.l  and  Fig.2 
can  be  as  low  as  0.  In  some  other  words,  the  peak  signal  to 
peak  signal  sidelobe  can  be  as  large  as  infinite.  In  addition,  it 
is  also  obvious  that  the  length  of  the  new  code  can  expend  to 
31  that  is  much  longer  than  the  length  of  the  Barker  code. 
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TABLE  I 

Optimum  Punctured  Binary  Sequence-pair 


Length 

Sequence 

(octet) 

Punctured  positions 

Energy 

Efficiency(%) 

3 

6 

3 

66.67 

5 

32 

3  4  5 

40.00 

5 

34 

2  4  5 

40.00 

7 

162 

4  5  7 

57.14 

7 

164 

4  6  7 

57014 

9 

652 

1234567 

22.22 

9 

760 

1234678 

22.22 

11 

3426 

4  5  6  8  11 

54.54 

11 

3550 

4  7  9  10  11 

54.54 

12 

7426 

1  6  7  12 

66.67 

12 

7550 

4  5  10  11 

66.67 

12 

7624 

3  6  9  12 

66.67 

13 

16606 

2  4  7  8  9  10  13 

46.15 

13 

17124 

5  6  8  9  10  12  13 

46.15 

15 

74232 

5  6  7  9  10  13  15 

53.33 

15 

75TT0 

6  7  10  II  IS  14  15 

5T33 

17 

351134 

4  6  7  8  9  10  12  16  17 

47.06 

17 

372142 

3  6  8  9  10  13  14  15  17 

47.06 

19 

1715412 

5  6  9  12  13  14  15  17  19 

52.63 

19 

1724154 

5  7  9  10  11  12  15  18  19 

52.63 

20 

3433330 

2  5  6  7  8  9  li  15  1<>  17 

18  19 

40.00 

20 

3610556 

1  6  7  8  9  10  11  16  17  18 

19  20 

40.00 

21 

7405316 

2  5  6  7  8  9  11  13  14  16 

17  20  21 

38.10 

21 

7563240 

3  5  6  9  10  12  13  15  17 

18  19  20  21 

38.10 

23 

37024632 

6  7  8  9  11  13  14  17  18 

21  23 

507 

23 

37263120 

6  8  11  12  15  16  18  20  2l 

22  23 

507 

28 

1702164566 

4  5  6  7  10  11  18  19  20  21 

24  25 

57.14 

28 

1734164226 

4  5  8  9  10  11  18  19  22  23 

24  25 

57.14 

28 

1740465534 

4  6  7  9  10  13  18  20  21  23 
24  27 

57.14 

29 

3556415302 

4  7  11  13  14  15  16  19  20 

21  24  25  26  27  29 

48.28 

29 

3642213634 

5  7  8  9  11  12  14  15  16  18 

23  24  26  28  29 

48.28 

31 

17053411166 

5  6  7  9  11  15  16  17  18  20 

21  23  24  28  31 

51.61 

31 

17464412730 

6  7  10  12  13  15  16  17 

18  20  22  26  29  30  31 

51.61 

Note:  The  sequences  in  the  table  arc  presented  by  octets,  *1*  for  *+1*,  ’0’ 
for  *-1’  and  punctured  positions  begins  from  left  to  right. 


Fig.  1.  Periodic  autocorrelation  property  of  punctured  binary  sequcncc-pair: 
(a)  13  length  code  (b)  31  length  code 


The  ambiguity  function  which  is  usually  used  to  analyze 
the  radar  performance  of  Doppler  shift  and  time  delay  can  be 
found  in  [17] 


y(t,FD)=  r 

J  —  oo 


x {s)ej2l’FDSx’{s  -  t)ds  =  A(t,FD) 


(7) 


where  t  is  the  time  delay  and  Fq  is  the  Doppler  shift. 

An  equivalent  definition  can  be  given  in  terms  of  the  signal 
spectrum  by  applying  basic  Fourier  transform  properties 


A(t,FD)  =  j  X'{F)X{F- FDy2”FtdF 

J  — oo 


(8) 


B.  Ambiguity  Function 

When  the  transmitted  impulse  is  reflected  by  a  moving 
target,  the  reflected  echo  signal  includes  a  linear  phase  shift 
which  corresponds  to  a  Doppler  shift  [17].  As  a  result  of  the 
Doppler  shift  /<*,  the  main  peak  of  the  autocorrelation  function 
is  reduced: 


[d]dB  =  10  log 


So  x(s)x*(s)da 
Jq  x(s)ei2l’fdTcx‘  ( s)ds 


(6) 


The  ambiguity  function  is  defined  as  the  magnitude  of 

M*>Fd)  [17] 

A(t,FD)  =  \A(t,FD)\  (9) 

However,  the  code  used  in  our  paper  is  sequence-pair,  so 
the  transmitting  code  and  the  receiving  code  are  not  the  same. 
The  ambiguity  function  can  be  obtained  as  following: 

/oo 

x(s)ej27rFDSy*(s  -  t)ds  (10) 


In  addition,  the  SNR  degraded  and  the  sidelobe  structure  is 
changed  thanks  to  the  Doppler  shift. 


APair{t,FD)  =  j"  X*(F)Y(F  -  FD)ej2*FtdF  (1 1) 

J  —  oo 
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The  ambiguity •  function  is  defined  as  the  magnitude  of 

Apmr  {f  Fd  ) 

^piur{t,Fo)  =  \Apuir(t,  Fp)\  (12) 

In  addition,  assume  the  punctured  binary  sequence-pair  is 

(x,y),  x(t)  =  En=o'  xn{i  -  nT0  and  2/(0  =  2A.(<  “ 

nTs),  since  the  periodic  correlation  is  used  instead  of  aperiodic 
correlation  in  this  paper.  The  Apnu  (£,  Fp)  in  one  period  of 
length  NTs  can  be  expressed  as 

A,,„ir(t,  Fd)  =  [  x(s  -  nTs)y'(s  -  nTs  -  (NTS  -  <)) 
Jo 

/•(N-l)T, 

eii2*FDs)ds  +  /  x(s-nTa)y'(s-nTs-0 

e(  j?*fds)(Is  (13) 

EXAMPLE 

Ambiguity  functions  of  the  punctured  binary  sequence- 
pair  within  length  of  13  and  31  used  in  the  last  section  are 
simulated,  where  maximal  time  delay  is  1  unit(normalized  to 
length  of  the  code,  in  units  of  NTS)  and  maximal  Doppler  shift 
is  5  units(normalized  to  the  inverse  of  the  length  of  the  code, 
in  units  of  1  /NTS).  The  ambiguity  functions  of  13-length 
long  Barker  code  and  31 -length  long  P4  code  are  presented  in 
Fig. 2  and  Fig. 3  in  order  to  compare  with  the  punctured  binary 
sequence-pairs  of  the  same  length 

Fig. 2  and  Fig. 3  show  that  the  sidelobe  improvement  of  the 
punctured  binary  sequence-pair  is  obvious  comparing  with 
those  of  Barker  code  and  P4  code  when  there  is  no  Doppler 
shift.  The  sidelobe  of  punctured  binary  sequence-pair  can 
reach  as  low  as  zero.  Nevertheless,  when  there  are  Doppler 
shift  and  time  delay,  the  ambiguity  functions  of  punctured 
binary  sequence-pair  is  not  as  flat  as  those  of  Barker  code 
or  P4  code.  In  some  other  words,  punctured  binary  sequence- 
pair  is,  to  some  extent,  less  tolerant  of  Doppler  shifts  than 
P4  code.  One  of  the  reasons  why  the  proposed  code  is  not 
tolerant  of  large  Doppler  shift  is  that  periodic  correlation 
property  is  used  for  punctured  binary  sequence-pair  instead 
of  the  aperiodic  correlation  property  which  is  used  for  the 
other  two  codes.  However,  comparing  with  P4  code,  punctured 
binary  sequence-pair  made  up  of  only  three  different  phases, 
is  more  simple  and  easy  to  obtain  in  the  industrial  world. 

C.  Doppler  Shift  Performance  Without  Time  Delay 

The  ambiguity  function  can  be  simplified  when  there  is  no 
time  delay. 

fTs 

Apair(0,FD)=  /  x(s  —  nTs)y‘(s  —  nTa)e^il,FDS^ds  (14) 
Jo 

However,  in  Fig. 2  and  Fig. 3,  it  is  not  obvious  to  see  the 
Doppler  shift  performance  of  punctured  binary  sequence-pair 
and  the  other  two  codes  when  there  is  no  time  delay.  The 
Doppler  shift  performance  without  time  delay  is  presented  in 
Fig. 4  and  Fig. 5. 


Fig.  2.  Ambiguity  function  of  13-lcngth  codes:  (a)  Punctured  binary 
scqucncc-pair  (b)  Barker  code 


Without  time  delay,  while  the  Doppler  shift  is  less  than 
1  unit  (normalized  to  length  of  the  code,  in  units  of  NTS), 
punctured  binary  sequence-pair  has  the  similar  performance  of 
Barker  and  P4  code  that  the  amplitude  has  a  sharp  downward 
trend.  Furthermore,  amplitude  of  punctured  binary  sequence- 
pair  decreases  more  quickly  than  amplitude  of  the  other  two 
codes.  However,  when  the  Doppler  shift  is  larger  than  1 
unit  (normalized  to  length  of  the  code,  in  units  of  NTS ), 
the  performances  of  these  codes  are  distinguished.  On  one 
hand,  the  trend  presented  by  punctured  binary  sequence-pair 
is  not  as  regular  as  the  other  two  kinds  of  code  when  the 
Doppler  shift  is  larger  than  1.  On  the  other  hand,  for  Barker 
and  P4  code,  when  Doppler  frequencies  equal  to  multiples 
of  the  pulse  repetition  frequency  (PRF  =  1/P/?/  =  1/Ts) 
the  ambiguity  value  turns  to  be  zero.  Because  of  these  zeros, 
such  multiples  of  the  pulse  repetition  frequency  will  render 
the  radar  blind  [l]to  their  velocities.  Nevertheless,  referring  to 
punctured  binary  sequence-pair,  ambiguity  values  do  not  go 
to  zero  when  Doppler  frequencies  are  equal  to  multiples  of 
the  PRF.  Therefore,  using  the  punctured  binary  sequence-pair 
as  the  compression  code  could,  to  some  extent,  improve  the 
blind  speed  problem  in  moving  target  detection  system. 
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(b) 


Fig.  3.  Ambiguity  function  of  3 1 -length  codes:  (a)  Punctured  binary 
scqucncc-pair  (b)  P4  code 


Fig.  4.  Doppler  shift  of  13-length  codcs(timc  delayed):  (a)  Punctured  binary 
scqucncc-pair  (b)  Barker  code 


V.  System  Simulation  in  Radar  System 

As  a  rule,  the  following  probabilities  [17]  are  of  most 
interest  in  the  Radar  system. 

1.  Probability  of  Detection,  Pd-  The  probability  that  a  target 
is  declared  when  a  target  is  in  fact  present. 

2.  Probability  of  False  Alarm,  Pfa-  The  probability  that  a 
target  is  declared  when  a  target  is  in  fact  not  present. 

3.  Probability  of  Miss,  Pm.  The  probability  that  a  target  is 
not  declared  when  a  target  is  in  fact  present. 

Note  that  Pm  =  1  -  Pd ,  thus,  Pd  and  Pfa  suffice  to  specify 
all  of  the  probabilities  of  interest  in  radar  system.  Therefore, 
the  above  three  probabilities  of  punctured  binary  sequence- 
pairs  in  radar  system  are  simulated  using  Matlab,  as  shown  in 
Fig. 6  and  Fig. 7.  In  addition,  the  performance  of  the  13-length 
Barker  code  and  31 -length  P4  code  are  provided  in  order  to 
compare  with  the  performance  of  punctured  binary  sequence- 
pairs  of  corresponding  lengths.  In  the  simulation  model,  we 
ran  Monte-Carlo  simulation  for  105  times  at  each  SNR  value, 
the  Doppler  shift  frequency  which  is  kept  less  than  1  unit 
(normalized  to  the  inverse  of  the  length  of  the  code,  in  units 
of  1  /NTS)  is  randomly  determined  by  Matlab,  and  the  time 
delay  is  assumed  to  be  zero.  We  use  threshold  detection  in 
coherent  system  and  the  threshold  is  adaptively  determined  in 
the  simulation. 

In  Fig. 6(a),  we  plotted  the  probabilities  of  miss  Pm  of  13- 


length  punctured  binary  sequence-pair  and  the  same  length 
Baker  code.  Observe  the  figure,  the  probability  of  miss  Pm 
of  the  system  using  13-length  punctured  binary  sequence- 
pair  is  lower  than  13-length  Barker  code  especially  when  the 
SNR  is  not  large.  It  is  in  accordance  with  the  result  shown 
in  Fig. 4  that  when  Doppler  shift  is  kept  less  than  1  and 
the  time  delay  is  zero,  the  amplitude  of  punctured  binary 
sequence-pair  falls  more  sharply  than  Barker  code.  In  Fig. 6(b), 
we  plotted  the  probabilities  of  miss  targets  detection  of  3 1  - 
length  punctured  binary  sequence-pair  and  those  of  the  same 
length  P4  code.  The  probability  of  miss  targets  detection  of  the 
system  using  31 -length  punctured  binary  sequence-pair  is  less 
than  31 -length  P4  code  especially  when  the  SNR  is  not  large. 
When  SNR  is  larger  than  17  dB,  both  probabilities  of  miss 
targets  detection  of  the  system  approach  zero.  However,  the 
probability  of  miss  targets  of  P4  code  is  lower  than  punctured 
binary  sequence-pair.  Comparing  both  Fig. 6(a)  and  Fig. 6(b), 
longer  punctured  code  performs  better  especially  when  the 
SNR  is  not  very  large. 

In  addition,  we  also  plotted  the  probability  of  detection 
versus  probability  of  false  alarm  of  the  coherent  receiver  in 
Fig-7. 

Fig. 7(a)  illustrates  performance  of  13-length  punctured  bi¬ 
nary  sequence-pair  and  the  same  length  Baker  code  when  the 
SNR  is  lOdB  and  14dB.  Having  the  same  SNR  value  such  as 
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Fig.  5.  Doppler  shift  of  3 1 -length  codcs(timc  dclay=0):  (a)  Punctured  binary  Fig.  6.  Probability  of  miss  targets  detection  (No  time  delay,  Doppler  shift  less 
sequence-pair  (b)  P4  code  than  1):  (a)  13-length  Punctured  binary  sequence-pair  VS.  13-length  Barker 

codc(b)  31 -length  Punctured  binary  sequence-pair  VS.  31 -length  P4  code 


lOdB  or  14dB  in  the  figure,  the  Pd  of  13-length  punctured  bi¬ 
nary  sequence-pair  is  larger  than  Pd  of  13-length  Barker  code 
while  the  Pfa  of  the  first  code  is  also  smaller  than  Pfa  of  the 
latter  code.  In  some  other  words,  13-length  punctured  binary 
sequence-pair  has  much  higher  target  detection  probability 
while  keeping  a  lower  false  alarm  probability.  Furthermore, 
observe  Fig.7(a),  13-length  punctured  binary  sequence-pair 
even  has  much  better  performance  at  lOdB  SNR  value  than 
13-length  Barker  code  at  14dB  SNR  value.  Similarly,  Fig. 7(b) 
shows  that  the  performance  of  31 -length  punctured  binary 
sequence-pair  is  much  better  than  the  P4  code  of  same  length. 
According  to  the  above  results,  it  is  easy  to  see  that  our  newly 
provided  punctured  binary  sequence-pair  is  very  promising  to 
be  an  alternative  pulse  compression  code  in  the  Radar  system. 

VI.  Conclusion 

A  kind  of  new  code-punctured  binary  sequence-pair  and 
its  properties  have  been  investigated  in  this  paper.  The  sig¬ 
nificant  advantage  of  the  new  code  over  conventional  phase 
compression  code  is  the  considerably  reduced  sidelobe  as 
low  as  zero  and  correspondingly  the  significantly  improved 
PSR  of  the  autocorrelation  function.  In  addition,  the  length  of 
punctured  binary  sequence-pair  known  can  be  as  long  as  31. 
The  disadvantage  of  the  new  sequence-pair  is  that  it  is  not 
so  tolerant  of  Doppler  shift  when  the  Doppler  shift  is  larger 


than  1.  We  also  apply  the  punctured  binary  sequence-pair  to 
the  target  detection  simulation  in  the  radar  system.  According 
to  the  simulation  results,  it  is  easy  to  observe  that  13-length 
punctured  binary  sequence-pair  has  better  performance  than 
13-length  Barker  code.  Similarly,  the  31 -length  punctured 
binary  sequence-pair  performs  better  than  31 -length  P4  code 
when  the  Doppler  shift  is  kept  less  than  1  unit(nomializcd 
to  the  inverse  of  the  length  of  the  code,  in  units  of  1  /NTS) 
in  the  radar  target  detection  system.  As  a  result,  the  general 
conclusion  from  the  results  presented  in  this  paper  is  that 
the  punctured  binary  sequence-pair,  which  has  much  longer 
code  length  and  better  autocorrelation  sidelobe  property  than 
the  biphase  code  such  as  Barker  code,  and  simpler  structure 
than  those  polyphase  codes  such  as  P4  code,  effectively 
increases  the  variety  of  candidates  for  pulse  compression  codes 
especially  for  long  code. 
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1.  ABSTRACT 

Based  on  zero  correlation  zone(ZCZ)  concept,  the  definition 
and  properties  of  ZCZ  sequence-pair  set  are  given  in  this 
paper.  A  method  is  presented  to  use  the  perfect  punctured 
sequence-pair  together  with  Hadamard  matrix  in  the  zero 
correlation  zone  to  construct  the  perfect  punctured  ZCZ 
sequence-pair  set.  According  to  performance  analysis ,  perfect 
punctured  ZCZ  sequence-pair  set  has  good  autocorrelation 
and  cross  correlation  properties  when  doppler  shift  is 
not  large.  The  radar  target  detection  system  simulation 
results  also  show  that  perfect  punctured  ZCZ  sequence-pair 
outperforms  other  conventional  pulse  compression  codes, 
such  as  the  well  known  polyphase  code-P4  code.  As  a 
result,  perfect  punctured  ZCZ  sequence-pair  set  can  be  good 
candidate  for  pulse  compression  code. 

II.  Introduction 

Pulse  compression,  known  as  a  technique  to  raise  the  signal 
to  maximum  sidelobe(signal-to-sidelobe)  ratio  to  improve  the 
target  detection  and  range  resolution  abilities  of  the  radar 
system,  allows  a  radar  to  simultaneously  achieve  the  energy 
of  a  long  pulse  and  the  resolution  of  a  short  pulse  without 
the  high  peak  power  which  is  required  by  a  high  energy  short 
duration  pulse  [1].  One  of  the  waveform  designs  suitable  for 
pulse  compression  is  phase-coded  waveform  design  that  a  long 
pulse  of  duration  T  is  divided  into  N  subpulses  each  of  width 
Ts.  Each  subpulse  has  a  particular  phase,  which  is  selected  in 
accordance  with  a  given  code  sequence.  The  pulse  compres¬ 
sion  ratio  equals  the  number  of  subpulses  N  =  T/Ts  %  BT , 
where  the  bandwidth  is  B  ~  1  /Ts.  In  general,  a  phase- 
coded  waveform  with  longer  code  word,  in  some  other  words, 
higher  pulse  compression  ratio,  can  have  lower  sidelobe  of 
autocorrelation,  relative  to  the  mainlobe  peak,  so  its  main 
peak  can  be  better  distinguished.  The  relative  lower  sidelobe 
of  autocorrelation  is  very  important  since  range  sidelobe  are  so 
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harmful  that  they  can  mask  main  peaks  caused  by  small  targets 
situated  near  large  targets.  In  addition,  the  cross-correlation 
property  of  the  pulse  compression  codes  should  be  considered 
in  order  to  reduce  the  interference  among  radars  when  we 
choose  a  set  of  pulse  compression  codes,  because  in  the  real 
world,  a  radar  may  not  work  alone,  such  as  in  the  RSN( Radar 
Sensor  Network). 

Much  time  and  effort  was  put  for  designing  sequences  with 
good  autocorrelation  and  cross  correlation  properties  for  radar 
target  ranging  and  target  detection.  However,  it  is  known 
that  for  most  good  binary  sequences  of  length  N(N  >  13) 
the  attainable  sidelobe  levels  are  approximately  \/N  [2]  [3] 
and  the  mutual  cross  correlation  peaks  of  sequences  of  the 
same  length  are  much  larger  and  are  usually  in  the  order  of 
2y/[N)  to  3 \/(N).  Set  of  binary  sequences  of  length  N  with 
autocorrelation  sidelobes  and  cross-correlation  peak  values 
both  of  approximately  y/( N)  are  only  achieved  in  paper  [4]. 
In  addition  to  binary  sequences,  polyphase  codes,  with  better 
Doppler  tolerance  and  lower  range  sidelobes  such  as  the  Frank 
and  PI  codes,  the  Butler-matrix  derived  P2  code  and  the 
linear-frequency  derived  P3  and  P4  codes  were  provided  and 
intensively  analyzed  in  [5]  [6]  [7].  Nevertheless,  the  range 
sidelobe  of  the  polyphase  codes  can  not  be  as  low  as  zero 
either.  The  structure  of  polyphase  codes  is  more  complicated 
and  is  not  easy  to  generate  comparing  with  binary  codes. 
Therefore,  based  on  ZCZ  (zero-correlation  zone)  [8]  concept, 
we  propose  the  ZCZ  sequence-pair  set,  which  can  reach  zero 
autocorrelation  sidelobe  during  zero  correlation  zone  and  zero 
mutual  cross  correlation  peaks  during  the  whole  period.  We 
also  propose  and  analyze  a  method  that  perfect  punctured 
sequence-pair  together  with  Hadamard  matrix  are  used  in  the 
zero  correlation  zone  to  construct  the  perfect  punctured  ZCZ 
sequence-pair  set  and  subsequently  apply  it  to  radar  detection 
system.  For  the  performance  evaluation  of  this  proposed  code, 
an  example  is  presented,  investigated  and  studied  in  the  radar 
targets  detection  simulation  system  in  the  paper.  Because 
of  the  ideal  property  performance  and  good  target  detection 
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performance  in  simulation  system,  our  proposed  new  code  can 
be  another  useful  candidate  for  pulse  compression  application, 
especially  in  the  RSN. 

The  rest  of  the  paper  is  organized  as  follows.  Section  2 
introduces  the  definition  and  properties  of  ZCZ  sequence-pair 
set.  In  Section  3,  the  perfect  punctured  ZCZ  sequence-pair  set 
is  introduced,  and  a  method  using  perfect  punctured  sequence- 
pair  and  Hadamard  matrix  to  construct  ZCZ  sequence-pair  set 
is  given  and  proved.  In  Section  4,  the  properties  and  ambiguity 
function  of  perfect  punctured  ZCZ  sequence  pair  set  are 
simulated  and  analyzed.  The  performance  of  perfect  punctured 
ZCZ  sequence-pair  set  is  investigated  in  radar  targets  detection 
simulation  system  comparing  with  P4  code.  In  Section  6, 
conclusions  are  drawn  on  perfect  punctured  ZCZ  sequence- 
pair  set. 

III.  The  Definition  and  Properties  of  ZCZ 
Sequence-pair  Set 

Zero  correlation  zone(ZCZ)  is  a  new  concept  provided 
by  Fan  [8]  [9]  [10]  in  which  the  autocorrelation  and  cross 
correlation  sidelobes  are  zero  while  the  time  delay  is  kept 
within  the  value  r  instead  of  the  whole  period  of  time  domain. 

We  consider  ZCZ  sequence-pair  set(X,  Y),  X  to  be  a  set  of 
K  sequences  of  length  N  and  Y  to  be  a  set  of  K  sequences 
of  the  same  length  N: 

*W€X  p  =  0, 1,2,  ...,K  —  1  (1) 

y{9)e  Y  q  =  0,1,2,..., K  -1  (2) 

The  autocorrelation  function  for  sequence-pair  (xv,yp) is 
defined  by: 

N-i 

(r)  =  £  0  <  r  <  N  -  1  (3) 

/  =  0 

The  cross  correlation  function  for  sequence-pair  (xp,yp) 
and  ( xq,yq),p  /  q  is  defined  by: 

N-i 

£*<->»,*>  (T)  =  Y  XlP)y(«'+r)modN-°  <T<N-\  (4) 

1=0 
N-  i 

CxMyM (r)  =  Y,  A((ri;)mo<yN,0  <  r  <  N-  1  (5) 

z=o 

For  pulse  compression  sequences,  some  properties  are  of 
particular  concern  in  the  optimization  of  any  design.  They 
are  peak  sidelobe  level,  the  energy  of  autocorrelation  side- 
lobes  and  energy  of  their  mutual  cross  correlation  [4],  So 
the  peak  sidelobe  level  which  represents  a  source  of  mutual 
interference  that  can  obscure  weaker  targets  can  be  presented 
as  maxx  |/?x/'yp(r)|,  r  E  Zo(zero  correlation  zone)  for  ZCZ 
sequence-pair  set.  The  another  optimization  criterion  for  the 
set  of  sequence-pair  is  the  energy  of  autocorrelation  sidelobes 
together  with  the  energy  of  cross  correlation.  By  minimizing 
the  energy,  it  can  be  distributed  evenly  and  the  peak  autocor¬ 
relation  level  can  be  minimized  as  well  [4].  Here,  the  energy 


of  ZCZ  sequence-pair  set  can  be  employed  as: 

K-\  Z0  K-  1K-1N-1 

£  =  ^x(P)y<P)(T)  +  ^I(,,)y(,,)(T)  (6) 

p=0  r  =  l  p=0  g=0  r=0 

(p  ^  q) 


According  to  (7),  it  is  obvious  to  see  that  the  energy  can 
be  kept  low  while  minimizing  the  autocorrelation  and  cross 
correlation  of  the  sequence-pair  set. 

Then,  the  ZCZ  sequence-pair  set  can  be  constructed  to 
minimize  the  autocorrelation  and  cross  correlation  of  the 
sequence-pair  set  and  the  definition  of  ZCZ  sequence-pair  set 
can  be  expressed: 

Definition  2-1  Assume(xjp\  y[p^)  to  be  sequence-pair  set 
of  length  N  and  number  of  sequence-pair  K ,  where  p  = 
1,2,3, N  -  l,i  =  0, 1, 2, K  -  1,  if  all  the  sequences 
in  the  set  satisfy  the  following  equation: 


N- 1 

N-\ 

^x(p)y«.)  (t)  =  Y  XiP) 

y (i+r)mod(N)  /  ^  *^'(i +r)mod(N) 

1=0 

/= 0 

f  A N, 

for  r  =  0,p  =  q 

=  \  0, 

for  r  =  0,p  ^  (7  (7) 

l  0, 

for  0  <  |r|  <  Zq 

where  0  <  A  <  1,  then  (x\p\y\p^)  is  called  ZCZ  sequence- 

pair  set,  ZCZP(N,  K ,  Zq)  is  a  symbol  in  short. 


IV.  Perfect  Punctured  ZCZ  Sequence-pair  Set 
A.  Definition  of  Perfect  Punctured  ZCZ  Sequence-pair  Set 

Matsufuji  and  Torii  have  provided  some  methods  of  con¬ 
structing  ZCZ  sequences  in  [1 1]  [12].  In  this  section,  we  apply 
perfect  punctured  sequence-pair  [13]  in  zero  correlation  zone 
to  construct  the  perfect  punctured  ZCZ  sequence-pair  set. 

Definition  3-1  [13]  Sequence  a  =  (uq,  ui,  •••:  wjv-i)  is  the 
punctured  sequence  for  v  =  (uq,  Vi, 

{0,  if  j  E  p  punctured  bits 
Vj ,  if  j  E  Non-punctured  bits 

Where  p  is  the  number  of  punctured  bits  in  sequence  vy 
suppose  Vj  =  (-1, 1),  u  is  p-punctured  binary  sequence,  (u,  v) 
is  called  a  punctured  binary  sequence-pair. 

Theorem  3-1  [13]  The  autocorrelation  of  punctured 
sequence-pair  (u,  v)  is  defined 

N- 1 

=  ^  ^  ^  ^  T  ^  N  —  1  (9) 

1=0 


If  the  punctured  sequence-pair  has  the  following  autocorre¬ 
lation  property: 


Ruv(t) 


£,  if  r  =  0modAr 
0,  others 


(10) 


the  punctured  sequence-pair  is  called  perfect  punctured 
sequence-pair  [13).  Where,  E  =  £/=o '  u,v(l+r)mmIN  = 
N  -  p,  is  the  energy  of  punctured  sequence-pair. 
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If  (x[p\  y[p^)  in  Definition  2-1  is  constructed  by  per¬ 
fect  punctured  sequence-pair  and  a  certain  matrix,  such  as 
Hadamard  matrix  or  an  orthogonal  matrix,  where 

x\p)  e  (-1,1),  i  =  0, 1,2, TV  —  1 
yt(,)  €  (-1,0,1),  1  =  0, 1,2,..., AT-  1 


N-\ 

N- 1 

II 

M 

y(i+T)modN  /  v 

i=o 

1=0 

(  an, 

for  r  =  0,p  =  q 

=  <  0, 

for  r  =  0,  p  7^  q 

l  0. 

for  0  <  |t|  <  Zo 

where  0  <  A  <  1,  then  (,x\p\y\p^)  can  be  called  perfect 
punctured  ZCZ  sequence-pair  set. 


B.  Design  for  Perfect  Punctured  ZCZ  Sequence-pair  Set 

Based  on  odd  length  perfect  punctured  binary  sequence 
pairs  and  a  Hadamard  matrix,  a  perfect  punctured  ZCZ 
sequence-pair  set  can  be  constructed  on  following  steps: 

Step  1 :  Considering  an  odd  length  perfect  punctured  binary 
sequence-pair  (u.u),  the  length  of  each  sequence  is  N \ 

u  =  u0,ui,...,uNl^i,Ui  €  (-1, 1), 

V  =  V0,  Vl,  VNi-liVi  €  (-1,0, 1), 
i  =  0,1,2 -  l,Ni  odd 


Step  2:  Consider  Walsh  sequences  set  B,  the  length  of  the 
sequence  is  N2  which  is  equal  to  the  number  of  the  sequences. 
In  some  other  words,  a  Hadamard  matrix  of  order  TV  2  is 
considered. 


B  =  (b0,b\...,bNt~l), 


E  =  %,b\,...,b'N]_i), 


Rb-b> 


f  N2,  i fi  =  j 
\  0,  if i?j 


Step  3:  Doing  bit-multiplication  on  the  perfect  punctured 
binary  sequence-pair  and  each  line  of  Walsh  sequences  set 
B(Hadamard  matrix),  then  sequence-pair  set  {X,Y)  is  ob¬ 
tained. 


b'  =  (&o,6i.->6jva-i)>*  =  0, 1,  -,N2  -  1, 

x)  =  WjmodN,tjmo<iNj,0  <  i  <  N2  -  1,0  <  j  <  N  -  1, 

X  =  (I0,!1,...,!"1'1), 

y )  =  vjmodN,bljmodNl,  0  <  i  <  N2  -  1,0  <  j  <  N  -  1, 
Y  =  (y°,y\...,  y"’-1) 

Where  GCD(N\,  N2)  =  1,  common  divisor  of  N\  and 
N2  is  1,  TV  =  N\  *  N2.  The  sequence-pair  set  (A”,y)  is 
perfect  punctured  ZCZ  sequence-pair  set  and  N\  —  1  is  the  zero 
correlation  zone  Z0.  The  length  of  each  sequence  in  perfect 
punctured  ZCZ  sequence-pair  set  \s  N  =  N\*  N2  that  depends 
on  the  product  of  length  of  perfect  punctured  sequence-pair 
and  the  length  of  Walsh  sequence  in  Hadamard  matrix.  The 


number  of  sequence-pair  in  perfect  punctured  ZCZ  sequence- 
pair  set  rests  on  the  order  of  the  Hadamard  matrix.  The 
sequence  x 1  in  sequence  set  X  and  the  corresponding  sequence 
y'  in  sequence  set  Y  construct  a  sequence-pair  (x\  y1)  that  can 
be  used  as  a  pulse  compression  code. 

The  correlation  property  of  the  sequence-pair  in  perfect 
punctured  ZCZ  sequence-pair  set  is: 

RX'y>(r )  =  RxJy'(T)  =  Ruv(TmodNl)Rb,b>(Tino(lN2) 

=  Rvv(TniodNi)Rb,b,(TTnodN2) 

(  EN2 ,  if  r  =  0,  t  =  j 
=  <  0,  if  0  <  |r|  <  TV,  -  l,i  =  j  (12) 

{  0,  'f  J 

where  N\  -  1  is  the  zero  correlation  zone  Z0. 

Proof: 

1)  When  ?  =  j, 

r  =  0, 

Ruv(0)  =  E,Rb>bj(0)  =  N2} 
Rx^(0)=Ruv(0)Rb^(0)  =  EN2] 
0<|r|<N1-l> 

Ruv(t)  =  0, 

Rxiy>(7)  =  Ruv{7^odN\)R^bJ(TinodN2)  =  0; 

2)  When  i  ±  j, 

r  =  0, 

Rt*  6>(0)  =  o, 

Rx'yA 0)  =  Rxjy*  (0) 

=  Ruv(7™odN\  )Rb*i>j  (TinodN2)  =  0; 

0  <  |r|  <  Ni  -  1, 

Ruv(7)  =  0, 

Rx'yj(r)  =  RxJy'(r) 

=  Ruv(77nodN\)Rbtbj(TmodN2)  =  0. 

According  to  Definition  2-1,  the  sequence-pair  set  con¬ 
structed  by  the  above  method  is  ZCZ  sequence-pair  set 

V.  Properties  of  Perfect  Punctured  ZCZ 
Sequence-pair  set 

Considering  the  perfect  punctured  ZCZ  sequence-pair  set 
that  is  constructed  by  the  method  mentioned  in  the  last  section, 
the  autocorrelation  and  cross  correlation  properties  can  be 
simulated  and  analyzed  with  Matlab.  For  example,  the  perfect 
punctured  ZCZ  sequence-pair  set  {X,Y)  is  constructed  by 
31 -length  perfect  punctured  binary  sequence-pair  ( u,v),u  = 
(1111-1-1-11  —  11-1111  —  1  —  1-1-11-1-11-1-1111- 
111  -  l),v  =  (1111000101011100001001001110110)  and 
Hadamard  matrix  H  —  [1111;1  —  11  —  1;  11  —  1  —  1;  1  —  1  —  11] 
of  order  4.  We  follow  the  three  steps  presented  in  Section  3  to 
construct  the  124-length  perfect  punctured  ZCZ  sequence-pair 
set.  The  number  of  sequence-pair  here  is  4  and  the  length  of 
each  sequence  is  31  *  4  =  124.  The  first  line  of  each  matrix 
X  =  [x\]x2]  x3;  £4]  and  Y  =  [yi ;  2/2;  2/3;  2/4]  constitute  a  cer¬ 
tain  perfect  punctured  ZCZ  sequence-pair  (x\,y\).  Similarly, 
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the  second  line  of  each  matrix  X  and  Y  constitute  another 
perfect  punctured  ZCZ  sequence-pair  (£2,2/2)- 

xx  =  (  1111-1-1-11-11-1111-1-1-1-1 

1-1-11-1-1111-111-11111-1-1 
-11-11-1111-1-1-1-11-1-11-1 
-mi  -  in  -  inn  - 1  - 1  -  n  -  li  -  in 
i-i-i  —  l-ii-i-ii-i-iiii-m-i 
nii-i-i-ii-ii-iiii-i-i-i-ii 
-l  -  li  -  l  -  mi  -  in  -  i  ), 
yx  =  (  11110001010111000010010011101101111000 
10101110000100100111011011110001010111 
00001001001 110110111 100010101 1 10000100 
1001110110  ); 

X2  =  (  1-11-1-11-1-1-1-1-1-11-1 

-11-1111-1-1-111-1111-1-1-1 
1-111-1111111-111-11-1-1-1111 
-1-11-1-1-1111-11-1-11-1-1 
-1  —  1  —  1-11—1-11-1111  —  1  —  1  —  11 
1-1111-1-1-11-111-1111111-111 
-11-1-1-1111-1-11-1-1-111), 
y2  =  (  1  -  11  -  1000-  10-  10-  11  -  10000100-  1 
001  -  1101  -  10  -  11  -  1100010101  -  110000 
-100100-  11  -  10  -  1101  -  11  -  1000-  10 
-10-  11  -  10000100-  1001  -  1101  -  10-1 
1  -  1100010101  -  110000-  100100-  11  -  10 
-no ). 

Perfect  punctured  ZCZ  sequence-pairs  (£1,2/1)  and  (£2,2/2) 
are  simulated  and  investigated  in  the  following  parts. 

A.  Autocorrelation  and  Cross  Correlation  Properties 

The  autocorrelation  property  and  cross  correlation  property 
of  124-length  perfect  punctured  ZCZ  sequence  pair  set  ( X ,  y) 
are  shown  in  Fig.l  and  Fig. 2  . 

From  the  Fig.l  and  Fig. 2,  the  sidelobe  of  autocorrelation 
of  ZCZ  sequence-pair  set  can  be  as  low  as  0  when  the  time 
delay  is  kept  within  Z0  =  N\  =  31  (zero  correlation  zone) 
and  the  cross  correlation  value  is  kept  as  low  as  0  during  the 
whole  time  domain. 

B.  Ambiguity’  function 

When  the  transmitted  impulse  is  reflected  by  a  moving 
target,  the  reflected  echo  signal  includes  a  linear  phase  shift 
which  corresponds  to  a  Doppler  shift  /j  [14].  As  a  result 
of  the  Doppler  shift  /(i,  the  main  peak  of  the  autocorrelation 
function  is  reduced  and  so  as  to  the  SNR  degradation  shown 
as  following: 


•  100  -so  0  SO  too  150 

Delayt/f^ 


Fig.  1 .  Periodic  autocorrelation  property  of  perfect  punctured  ZCZ  sequence- 
pair  set 
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Fig.  2.  Periodic  cross  correlation  property  of  perfect  punctured  ZCZ 
sequence-pair  set 


In  addition,  the  sidelobe  structure  is  changed  because  of  the 
Doppler  shift. 

The  ambiguity  function  which  is  usually  used  to  analyze 
the  radar  performance  within  Doppler  shift  can  be  found  in 
[14]  shown  as  following: 

y(t,FD)  =  [°°  x(sy2"FDSx-(s  -  t)ds  =  A(t,FD)  (14) 
7-00 

where  t  is  the  time  delay  and  Fp  is  the  Doppler  shift 

An  equivalent  definition  can  be  given  in  terms  of  the  signal 
spectrum  by  applying  basic  Fourier  transform  properties 

/oo 

X‘(F)X(F-  FD)ej2nF,dF  (15) 

-OO 

The  ambiguity  function  is  defined  as  the  magnitude  of 
A(t,Fp)  [H]  as  following: 

A(fFp)  =  \A(CFd)\  (16) 


\d\dD  =  10  log 


Jo  x(s)xm(s)ds 
x(s)ej2nfdTcx*  (s)ds 


(13) 


However,  we  have  ZCZ  sequence-pair  set  in  the  paper,  so 
the  transmitting  code  and  the  receiving  code  are  not  the  same. 
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The  ambiguity  function  for  ZCZ  sequence-pair  set  can  be 
rewritten  as 

/oo 

xip)(sy^o.^M-(s  _  t)ds 

■  oo 

where  p,q  =  0.1, 2...,  K  -  1 


A  pair  [t ,  t  D) 


X 


<  ,  .lr\  .  _ 

(Fjy  v9j(F 


FD)ej2~FtdF 


where  p,  q  =  0, 1, 2...,  K  —  1  (17) 

The  ambiguity'  function  is  defined  as  the  magnitude  of 
Ap(nr(t'FD)  as  following: 

Aparr(t,FD)  =  \  Apair  (t ,  Fd)\  (18) 

In  addition,  assume  the  ZCZ  sequence-pair  set  are  (X.K), 
x(p)  6  x,  a:<">  =  EnJo  -  nTs)  and  ?/<»>  e  V,  y <«>  = 

ijn’1  (t  -  nT$),  since  the  periodic  correlation  is  used 
instead  of  aperiodic  correlation  in  this  paper.  The  j4paT>(t,  Fd) 
in  one  period  of  length  yVTs  can  be  expressed  as: 

\APair(t,FD)\  (19) 

=  I  /  u(p)(s  -  nTs)v^’(s  -  riT .  -  ((VTS  -  t)) 

Vo 

AN-dt, 

eU2*lD,)ds  +  /  U('’>(s-  nTs)v{q)’ 

(s  -  tiTj  -  0e°'2,rFDs)ds|,  p,q  =  0, 1.2...,  K  -  1 

At  the  same  time,  when  p  =  q,  equation  (20)  can  be  used  to 
analyze  the  autocorrelation  performance  within  Doppler  shift; 
and  when  q  p,  equation  (20)  can  be  used  to  analyze  the 
cross  correlation  performance  within  Doppler  shift.  Equation 
(20)  is  plotted  in  Fig. 3  in  a  three-dimensional  surface  plot 
to  analyze  the  radar  performance  of  perfect  punctured  ZCZ 
sequence-pair  set  within  Doppler  shift.  Here,  maximal  time 
delay  is  1  unit  (normalized  to  length  of  the  code,  in  units  of 
NTS)  and  maximal  doppler  shift  is  5  units  for  cross  correlation 
and  3  units  for  autocorrelation  (normalized  to  the  inverse  of 
the  length  of  the  code,  in  units  of  l/NTs). 

In  Fig. 3(a),  there  is  relative  uniform  plateau  suggesting  low 
and  uninform  sidelobes,  minimizing  target  masking  effect  in 
zero  correlation  zone  of  time  domain,  where  Zo  =  31,  -31  < 
r  <  31.  From  Fig. 3(b),  considering  cross  correlation  property 
between  any  two  perfect  punctured  ZCZ  sequence-pairs  among 
the  ZCZ  sequence-pair  set,  124-length  perfect  punctured  ZCZ 
sequence-pair  set  is  tolerant  of  Doppler  shift  when  Doppler 
shift  is  not  large.  When  the  Doppler  shift  is  zero,  the  range 
sidelobe  of  cross  correlation  of  our  proposed  code  is  zero  in 
the  whole  time  domain. 

As  synchronization  technology  develops  quickly  in  the 
industrial  world,  time  delay  can,  to  some  extent,  be  well 
limited,  it  is  necessary  to  investigate  the  performance  of  our 
proposed  code  without  time  delay.  When  t  =  0,  the  ambiguity 
function  can  be  expressed  as: 

ipa,r(0,  Fd)=  [  ’  x(s  -  nTs)y"(s  -  nTJe^^ds  (20) 
Jo 
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Fig.  3.  Ambiguity  function  of  124-lcngth  ZCZ  scqucncc-pair  set:  (a) 
autocorrelation  (b)  cross  correlation 


And  the  doppler  shift  performance  without  time  delay  is 
presented  in  the  Fig. 4. 

Fig. 4(a)  illustrates  that  without  time  delay,  when  the 
Doppler  shift  is  less  than  1  unit,  autocorrelation  value  of 
perfect  punctured  ZCZ  sequence-pair  set  falls  sharply  and 
the  trend  of  the  amplitude  over  the  whole  frequency  domain 
decreases  as  well.  Fig.4(b)  shows  that  there  are  some  convex 
surfaces  in  cross  correlation  performance.  Observe  Fig. 4(a) 
and  Fig.4(b),  when  Doppler  frequencies  equal  to  multiples 
of  the  pulse  repetition  frequency  (PRF  =  1  /PRI  =  1  /Ts) 
except  Doppler  frequency  is  equal  to  2  PRF  for  cross  corre¬ 
lation,  all  the  ambiguity  value  turns  to  be  zero.  That  is  the 
same  as  many  widely  used  pulse  compression  binary  code 
such  as  Barker  code.  In  a  word,  behavior  in  time  delay  and 
in  response  to  Doppler  shift  shown  by  ambiguity  function  of 
perfect  punctured  ZCZ  sequence-pair  set  can  be  as  good  as 
conventional  pulse  compression  binary  code. 

VI.  System  Simulation  in  Radar  System 

According  to  [14],  Po(Probalfihty  of  Detection), 
PpA (Probability  of  False  Alarm)  and  Fa/ (Probability  of 
Miss)  are  three  probabilities  of  most  interest  in  the  radar 
system.  Note  that  Pm  =  1  —  Pl>,  thus,  Po  and  Pfa  suffice 


5 


Page  132  of  134 


OM 


I 

l  “j 

«  031 
Z  02; 

oi( 

0‘ 

0 


05 


1  15  2  25  3 

Dopper  shift  v  * Mtb 


(a) 


07 


(b) 


Fig.  4.  Doppler  shift  of  124-length  ZCZ  sequence-pair  sct(timc  dclay=0): 
(a)  autoeorrclation  (b)  cross  correlation 
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Fig.  5.  Probability  of  miss  targets  detection  (No  time  delay,  Doppler  shift 
less  than  1):  124-lcngth  perfect  punctured  ZCZ  sequence-pair  VS.  1 24-lcngth 
P4  code 
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Fig.  6.  Probability  of  detection  versus  probability  of  false  alarm  of  the 
coherent  receiverfNo  time  delay,  Doppler  shift  less  than  I) :  124-lcngth  perfect 
punctured  ZCZ  sequence-pair  VS.  124-lcngth  P4  code 


to  specify  all  of  the  probabilities  of  interest  in  radar  system. 
Therefore,  we  simulated  the  above  three  probabilities  of 
perfect  punctured  ZCZ  sequence-pair  in  radar  system  in 
this  section.  The  performance  of  124-length  P4  code  is 
also  provided  in  order  to  compare  with  the  performance  of 
punctured  binary  perfect  punctured  ZCZ  sequence-pair  of 
corresponding  length.  In  the  simulation  model,  105  times  of 
Monte-Carlo  simulation  has  been  run  for  each  SNR  value. 
The  Doppler  shift  frequency  that  is  kept  less  than  1  unit 
(normalized  to  the  inverse  of  the  length  of  the  code,  in 
units  of  1  /NTS)  is  a  random  variable,  and  the  time  delay  is 
assumed  zero.  Threshold  detection  is  used  in  this  coherent 
system  where  the  threshold  is  adaptively  adjusted. 

In  Fig. 5,  the  probabilities  of  miss  target  detection  Pm  of  the 
system  using  124-length  perfect  punctured  ZCZ  sequence-pair 
are  lower  than  124-length  P4  code  especially  when  the  SNR  is 
not  large.  When  SNR  is  larger  than  18  dB,  both  probabilities 
of  miss  targets  of  the  system  approaches  zero.  However,  the 
probabilities  of  miss  targets  of  P4  code  fall  more  quickly  than 
perfect  punctured  ZCZ  sequence-pair. 

We  plotted  the  probability  of  detection  Pd  versus  probabil¬ 
ity  of  false  alarm  PpA  of  the  coherent  receiver  in  Fig. 6. 


Fig. 6  shows  performance  of  124-length  perfect  punctured 
ZCZ  sequence-pair  and  the  same  length  P4  code  when  the 
SNR  is  lOdB  and  14dB.  Within  the  same  SNR  value  ei¬ 
ther  lOdB  or  14dB,  the  probabilities  of  detection  of  124- 
length  perfect  punctured  ZCZ  sequence-pair  are  much  larger 
than  probabilities  of  detection  of  124-length  P4  code,  and 
meanwhile  probabilities  of  false  alarm  of  the  first  code  arc 
also  smaller  than  PpA  of  the  latter  code.  In  some  other 
words,  124-length  perfect  punctured  ZCZ  sequence-pair  has 
higher  target  detection  probability  while  keeping  a  lower 
false  alarm  probability.  Furthermore,  observe  Fig. 6,  124-length 
perfect  punctured  ZCZ  sequence-pair  even  has  much  better 
performance  with  lOdB  SNR  than  124-length  P4  code  with 
14dB  SNR. 

Vll.  Conclusion 

The  definition  and  properties  of  ZCZ  sequence-pair  set  are 
discussed  in  this  paper.  Based  on  perfect  punctured  sequence- 
pair  and  Hadamard  matrix,  a  constructing  method  for  the 
perfect  punctured  ZCZ  sequence-pair  set  has  been  investigated 
along  with  its  properties.  The  significant  advantage  of  the 
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perfect  punctured  ZCZ  sequence-pair  set  is  a  considerably 
reduced  sidelobe  as  low  as  zero  in  the  zero  correlation  zone 
and  zero  mutual  cross  correlation  value  in  the  whole  time 
domain.  The  disadvantage  of  our  proposed  code  is  that  the 
number  of  the  sequences  in  the  set  depends  on  the  order 
of  Hadamard  matrix  that  is  limited  by  2 k,(k  =  0,1,...). 
According  to  the  radar  system  simulation  results  shown  in 
Fig. 5  and  Fig. 6,  it  is  easy  to  observe  that  124-length  perfect 
punctured  ZCZ  sequence-pair  set  has  better  performance  than 
124-length  P4  code  when  the  Doppler  shift  is  kept  less  than 
1  unit(normalized  to  the  inverse  of  the  length  of  the  code, 
in  units  of  \/NT s).  As  a  result,  the  general  conclusion  can 
be  drawn  from  the  results  presented  in  this  paper  that  the 
perfect  punctured  ZCZ  sequence-pair  set,  which  has  much 
better  autocorrelation  and  cross  correlation  properties  than  the 
optimum  biphase  codes(longer  than  13),  whose  autocorrelation 
sidelobes  and  cross  correlation  peak  value  have  been  found 
to  be  both  approximately  \/]V,  can  effectively  increase  the 
variety  of  candidates  for  pulse  compression  codes  Because 
of  the  ideal  cross  correlation  properties  of  perfect  punctured 
ZCZ  sequence-pair  set,  future  work  would  be  focused  on  the 
application  of  the  perfect  punctured  ZCZ  sequence-pair  set  in 
multiple  radar  system. 
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